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Presented  herein  is  a  user's  manual  for  the  ELSHOK  computer  code.  This 
computer  program  calculates  the  transient  response  of  a  submerged  ring-stiffened 
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fluid  is  based  upon  the  Doubly  Asymptotic  Approximation. 


SUMMARY 


Presented  herein  is  a  user's  manual  for  the  ELSHOK  computer  code. 
This  computer  program  calculates  the  transient  response  of  a  submerged 
ring-stiffened  shell  of  revolution,  with  or  without  internal  structure, 
to  an  underwater  shock  wave  having  an  arbitrary  direction  of  impingement 
Linearly  elastic  structures  are  considered,  and  the  surrounding  fluid  is 
treated  as  an  infinite  acoustic  medium.  Modal  structural  analysis  is 
employed  throughout.  When  the  shell  contains  internal  structure,  a 
substructuring  technique  is  applied  in  which  the  modes  of  the  empty  ring 
stiffened  shell  and  the  fixed-base  modes  of  each  piece  of  internal 
structure  are  coupled  through  the  use  of  dynamic  boundary  conditions. 

The  structure-fluid  interaction  is  approximated  by  means  of  the  Doubly 
Asymptotic  Approximation  expressed  in  terms  of  functions  which  are 
orthogonal  over  the  wet  surface  of  the  submerged  shell. 
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I.  INTRODUCTION 


The  development  of  analytical  methods  and  computer  codes  for  the  analysis 
of  the  response  of  submerged  structures  to  underwater  shock  waves  has  made  con¬ 
siderable  progress  In  recent  years.  In  particular,  Weldllnger  Associates,  under 
the  sponsorship  of  the  Defense  Nuclear  Agency  (DNA)  and  the  Office  of  Naval 
Research  (ONR) ,  has  developed  the  ELSHOK  (ELASTIC  SHOCK)  computer  code  for 
analyzing  such  problems  when  only  elastic  structural  deformations  are  involved. 
For  the  case  of  more  Intense  loadings  In  which  elasto-plastic  structural  defor¬ 
mations  occur,  the  EPSA  code  of  Ref.  [1]  may  be  used. 

The  present  report  describes  the  ELSHOK  computer  code.  This  computer 
program  calculates  the  transient  response  of  a  submerged  ring-stiffened  shell 
of  revolution,  with  or  without  internal  structure,  to  an  underwater  shock  wave 
emanating  from  an  explosive  source  placed  at  an  arbitrary  location  away  from 
the  shell.  Linearly  elastic  structures  are  considered,  and  the  surrounding 
fluid  is  treated  as  an  infinite  acoustic  medium.  Modal  structural  analysis 
is  employed  throughout.  When  the  shell  contains  internal  structure,  the  sub- 
structuring  technique  of  Refs.  [2]  and  [3]  is  applied  in  which  the  free-free 
modes  of  the  empty  ring-stiffened  shell  and  the  fixed-base  modes  of  each  piece 
of  internal  structure  are  coupled  through  the  use  of  dynamic  boundary  conditions. 
The  structure-fluid  interaction  is  approximated  in  ELSHOK  by  means  of  the 
Doubly  Asyiiq>totic  Approximation  (DAA)  of  Ref.  [4],  expressed  in  terms  of 
functions  which  are  orthogonal  over  the  wet  surface  of  the  submerged  shell 
(Ref.  [2]). 

The  development  of  ELSHOK  has  been  in  tandem  with  an  experimental 
program  aimed  ".t  providing  a  high  quality  data  base  for  judging  com¬ 
puterized  preu  .  on  mods.  This  combined  analytical/experimental 


7 


effort  continues  to  address  the  submerged  shock  response  of  structures  of 
increasing  complexity.  To  date,  comparisons  of  test  data  and  calculations 
performed  using  ELSHOK  have  shown  that  ELSHOK  is  a  viable  design/analysis 
tool. 

Documentation  for  the  current  operational  version  of  ELSHOK  is  provided  in 
this  manual  for  theory,  code  design  and  implementation,  input  data,  operating 
instructions,  data  management,  and  two  sample  problems.  It  is  hoped  that  the 
information  supplied  here  is  sufficiently  detailed  to  enable  a  new  user  to 
perform  an  ELSHOK  calculation.  The  authors  wish  to  broaden  the  ELSHOK  user 
community  and  to  establish  communication  links  between  users  and  developers 
so  that  the  ELSHOK  code  may  be  enhanced  to  satisfy  future  requirements. 

It  should  be  pointed  out  at  this  time  that  modifications  to  ELSHOK  have 
been  made  to  permit  the  use  of  geometrically  non-linear  members  as  shock 
mounts.  These  mounts  may  serve  as  connections  between  (1)  internal  structure 
and  the  ring-stiffened  shell  and  (2)  pieces  of  internal  structure.  Thus, 
the  next  generation  of  the  ELSHOK  code  will  allow  hard-mounted,  soft  (spring) - 
mounted,  and  shock-mounted  Internal  appendages  to  be  treated  in  a  routine 
manner.  The  Improved  version  of  ELSHOK  will  be  released  after  the  completion 
of  extensive  checking,  and  a  revised  user's  manual  will  then  be  Issued. 


II.  THEORY 


a.  General  Information 

The  ELSHOK  computer  code  is  a  suite  of  computer  programs  assembled 
and  developed  to  calculate  the  submerged  shock  response  of  a  linearly 
elastic  structure  subjected  to  an  underwater  shock  wave.  Specifically, 
a  ring-stiffened  shell  of  revolution  of  finite  length,  with  or  without 
internal  structure,  is  considered  to  be  immersed,  initially  at  rest,  in 
an  infinite  acoustic  fluid  and  to  be  excited  by  an  acoustic  spherical  wave, 
of  arbitrary  pressure  profile,  emanating  from  an  explosive  source  located  at 
an  arbitrary  point  away  from  the  structure.  It  is  required  to  compute 
the  transient  response  of  the  shell  and  any  internal  appendages  and  to 
provide  results  in  the  form  of  printed  and/or  plotted  velocity-time  his¬ 
tories  and  associated  pseudo-velocity  shock  spectra. 

In  ELSHOK,  the  structure-fluid  interaction  is  approximated  by  means 
of  the  DAA  of  Ref.  [4].  Following  Ref.  [2],  the  form  of  the  DAA  used  is 
that  obtained  when  the  normal  fluid  displacement  of  the  structure-fluid 
interface  is  expanded  in  a  series  of  functions  (surface  expansion  functions) 
which  are  orthogonal  over  the  wet  surface  of  the  structure.  The  ele¬ 
ments  of  the  matrices  in  the  resulting  DAA  are  obtained  by  matching  exact 
pressure-velocity  relations  at  zero  and  infinite  frequencies.  Thus,  in 
transient  problems,  the  DAA  yields  exact  results  at  early  and  late  times 
and,  by  the  nature  of  its  formulation,  brings  about  a  smooth  transition 

between  these  two  limits.  The  DAA  accounts  for  the  effects  of  the 
fluid  by  quantities  defined  solely  on  the  wet  surface  of  the  submerged 
body  and  may  be  thought  of  as  uncoupling  the  fluid  field  from  the 
structural  field. 
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The  ELSHOK  code  determines  the  transient  response  of  the  submerged 


structural  system  described  above  by  using  component  modal  analysis,  or  sub¬ 
structuring,  as  discussed  in  Ref.  [3].  In  this  technique,  the  complete 
structural  system  is  assumed  to  be  comprised  of  a  main  structure  S  (the  shell) 
and  any  attached  substructures  or  Internal  appendages  a.  The  vibration  modes 
of  each  component  are  determined  separately,  and  the  equations  of  motion  of 
the  system  S  +  a  are  obtained  by  enforcing  compatibility  of  deformation  at 
the  points  of  attachment.  The  implementation  of  substructuring  in  ELSHOK  also 
allows  for  the  case  of  no  attached  appendages  (i.e.,  for  an  empty  shell  S) . 
b.  Response  Equations 

As  shown  in  Refs.  [2]  and  [3],  the  use  of  the  DAA  permits  the  equations 
governing  the  shock  response  of  a  submerged  body  S,  with  internal  structure 
a,  to  be  written  as 


PjO./  Xy  (Qj  +  . 


(1) 


q^(0)  =  q.(0)  =  0 


(j  =  1.2 . Vg)  (2) 


pcpj  ‘a 


I 


Q^dt  - 


J.  V’i' 


(j  ^ 


(3) 


Qjdt 


(j  =  1,2, .. . ,Vp)  (4) 
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where  the  are  the  generalized  coordinates  of  the  main  structure  S, 

and  are,  respectively,  the  free-free  natural  frequencies  and  corresponding 

generalized  masses  of  S,  p  and  c  denote,  respectively,  the  mass  density 

of  and  speed  of  sound  in  the  fluid,  v  and  v  represent,  respectively, 

S  F 

the  number  of  shell  modes  and  surface  expansion  functions  used,  and  a  dot 
denotes  differentiation  with  respect  to  time  t.  Additional  symbols  appear¬ 
ing  in  the  above  equations  will  be  defined  below.  It  may  be  observed  that 
the  effect  of  the  fluid  in  the  interaction  problem,  as  described  by  the 
DAA  of  Eq.  (3),  is  simply  to  apply  an  additional  loading  on  the  shell  S,  result- 
ing  from  the  terms  in  Eq.  (1).  No  further  analysis  of  the  fluid  field  is 

required  during  the  structure-fluid  interaction  problem. 

The  elements  of  the  fluid-shell  transformation  matrix  appearing  in 


Eqs.  (1)  and  (3)  are  given  by 


■r  i 


4'.  dA 

J 


(1  l»2,...,Vg  ^  j  —  l,2,...,Vp) 


j  A 


where  ate  the  normal  components  of  the  free-free  in-vacuo  mode  shapes 
of  the  empty  structure  S  (i.e.,  without  internal  structure),  are  the 
surface  expansion  functions,  A  is  the  wet  surface  area  of  the  body  S, 
and  the  y^  are  normalization  constants  defined  by 


(i»j  *  1*2, ...,Vp) 


in  which  6^^  is  the  Kronecker  delta. 
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If  Pj  and  Uj.  represent,  respectively,  the  incident  pressure  and  normal 
fluid  particle  velocity  on  the  wet  surface  of  S,  the  expansion  coefficients 
associated  with  the  incident  shock  wave  [Eqs.  (1)  and  (3)]  may  be  expressed 


(3  ®  1,2 , . . . »Vp) 


U  =  —  .  dA 

IJ  ^  )  I  J 

A 


(j  “  1,2 , . . . »V_) 


The  ej£  terms  in  Eq.  (3)  are  the  elements  of  the  inverse  of  the  symmetrical 

virtual  mass  matrix  formulated  in  terms  of  the  <p . .  Moreover,  the  Q  .  in 

3  03 

Eq.  (1)  are  the  generalized  forces  exerted  on  the  body  S  by  the  internal 
structure  a.  These  terms  contain  the  effects  of  the  elastic  reactions 
of  c  on  S  and  of  the  inertial  forces  due  to  any  concentrated  masses  joined 
to  S.  General  matrix  equations  for  the  are  given  in  Ref.  [3]. 

If  a  substructure  a  is  attached  to  S  without  the  use  of  inter¬ 
vening  energy-absorbing  devices  or  shock  mounts.  Ref.  [3]  shows  that  the 
physical  degrees  of  freedom  of  o  may  be  partitioned  into  the  two  sets 


i=l 


y  C  q 
ii^i 


(j  =  1,2 . N) 


‘^aji'^ai  ®ji^' 


ji  ai 


(j  =  1,2,. ..,N)  (10) 
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where  the  are  the  N  physical  degrees  of  freedom  of  a  which  are  con¬ 

strained  to  move  with  S,  the  are  the  remaining  N  active  physical 
degrees  of  freedom  of  a,  and  the  are  the  generalized  coordinates  of 
a  .  In  Eqs.  (9)  and  (10),  the  are  coefficients  obtained  from  the 
modes  of  S  resulting  from  the  enforcement  of  compatibility  of  deformation 
at  the  points  of  attachment  of  a  to  S,  the  are  the  elements  of  the 

fixed-base  modes  of  a,  and  the  g^^  are  influence  coefficients  (sometimes 
termed  "constraint  modes")  accounting  for  the  support  or  base  motions  of  o. 

As  discussed  in  Ref.  [3],  the  partitioning  of  Eqs.  (9)  and  (10)  leads 
to  the  following  equations  for  the  generalized  coordinates  of  an  elastic 
Internal  appendage  a  connected  to  a  main  structure  S: 


N  "S 

‘  uL  4=1 


(j  "  1,2, .. .  >'^p) 


(j  -  1 >2 , . . . »^p) 


in  which  the  are  the  fixed-base  natural  frequencies  of  the  internal 
structure  o  and  the  are  terms  resulting  from  the  support  on  base 
motions  of  o  .  Reference  [3]  contains  general  matrix  equations  for  the 


G  ,  ,  the  C  ,  and  the  g,,  .  Response  equations  corresponding  to  Eqs.  (11) 

JK  JK  JK 

and  (12)  are  required  for  each  substructure  included  in  a  given  analysis. 


c.  Surface  Expansion  Functions 


Surface  expansion  functions  are  spatial  functions  having  the  property 
of  orthogonality  over  the  wet  surface  of  a  submerged  body.  As  pointed 
out  in  Ref.  [2],  when  surface  expansion  functions  are  used  in  a  series 
expansion  of  the  normal  displacement  of  the  wet  surface,  they  lead  to 
matrices  in  the  DAA  which  are  we 11 -conditioned  for  inversion.  The  sur¬ 
face  expansion  functions  used  in  any  calculation  must  be  capable  of 
representing  the  normal  displacement  of  the  structure-fluid  interface. 

In  ELSHOK,  the  determination  of  the  surface  expansion  functions 
requires  the  wet  surface  of  the  structure  under  study  to  be  partitioned 
into  a  central  region  (which  need  not  be  a  right  circular  cylinder)  and 
two  end  closures  of  arbitrary  revolved  shape,  as  shown  in  Fig.  1.  When 
symmetry  is  employed  along  the  length  of  the  structure,  only  the  left 
end  closure  must  be  specified.  Since  ELSHOK  allows  a  user  to  highlight 
a  portion  of  the  structure  when  localized  responses  are  expected,  a 
special  section  of  the  structure,  called  the  "compartment  region"  or 
"compartment  model",  may  also  be  specified.  The  use  of  the  compartment 
region  is  discussed  in  detail  in  Appendices  A  and  B. 

Points  on  the  wet  surface  are  located  by  specifying  the  meridional 
arc  length  s,  measured  along  the  reference  curve  used  to  generate  the 
wet  surface  of  revolution,  and  the  circumferential  coordinate  angle  0 
(Fig.  1).  Following  the  shell-of-revolution  analysis  employed  by  ELSHOK 
for  the  main  body  S,  the  9-dependence  of  the  surface  expansion  functions 
is  eliminated  by  means  of  separation  of  variables.  Thus,  the  circumferential 
distributions  of  these  functions  are  taken  to  be  cos  N9  or  sin  N0,  as 
appropriate,  where  N  is  the  integral  number  of  circumferential  waves  or 


circumferential  harmonics.  The  meridional  distributions  of  the  surface 


expansion  functions  may  then  be  obtained  for  one  given  value  of  N  at  a 
time. 


For  N=0  (axisymmetrical)  and  N  >  1,  the  meridional  distribution 
of  a  convenient  set  of  functions  for  generating  orthogonal  functions 
capable  of  describing  the  normal  motion  of  the  central  region  of 
the  structure  illustrated  in  Fig.  1  may  be  written  symbolically  as 
=  cos  (M  IT  z/L)  (M  =  0,1,2>...)  on  central  region 


(1 


=  0  elsewhere 


where  L  is  the  meridional  arc  length  of  the  central  region,  M  is  the 
number  of  meridional  half-waves  in  L,  and  z  is  the  meridional  coordinate 
measured  from  the  left  end  of  the  central  region  (0  <  z  <  L) . 


For  N=0,  the  meridional  distribution  of  a  convenient  set  of  functions 
for  generating  orthogonal  functions  capable  of  describing  the  normal 
motion  of  the  end  closures  may  be  written  symbolically  as 

5 ,  =  a  cos  (^)  (ii  =  0,2,4,...)  on  each  end  closure 

(1 

F  =0  elsewhere 
k 

in  which  SL  is  the  meridional  arc  length  of  the  generator  of  an  end 
closure,  M  is  the  number  of  quarter-waves  in  £,  a^  =  +  1  is  a  factor 
used  to  construct  symmetrical  and  antisymmetrical  functions,  and  s 
is  the  meridional  coordinate  measured  from  the  appropriate  pole 
(0  <  s  £  H) .  For  N  ^  1,  Eq.  (14)  must  be  replaced  by 


5 


(15) 


£  *  a  sin  (M  =  1,3,5,...)  on  each  end  closure 

k  c  Zx, 


£,  =  0  elsewhere 
k 

In  order  to  avoid  a  directional  singularity  at  the  pole  of  each  end 
closure  when  the  circumferential  distribution  is  supplied. 


For  N=0  (axisyimnetrical)  and  N  >  1,  the  meridional  distribution 
of  a  convenient  set  of  functions  for  generating  orthogonal  functions 
capable  of  describing  the  localized  normal  motion  of  the  compartment 
region  may  be  written  symbolically  as 


*  sin  (M  TT  z/L  )  (M  =  1,2,3,...)  on  compartment  region 


=  0  elsewhere 

where  L  is  the  meridional  arc  length  of  the  compartment  region,  M  is  the 
number  of  meridional  half-waves  in  L,  and  z  is  the  meridional  coordinate 


measured  from  the  left  end  of  the  compartment  region  (0  £  z  £  L) . 


In  the  usual  application  of  ELSHOK,  a  central  region  and  two  end  closures 
(or  one  if  symmetry  is  employed)  are  considered  without  the  use  of  a  compart¬ 
ment  region.  For  this  case,  the  surface  expansion  functions  are  generated 
from  the  convenient  functions  of  Eqs.  (13)-(15)  for  each  required  value  of  N. 
When  a  compartment  region  is  employed  in  addition  to  the  central  and  end- 
closure  regions,  the  representation  of  the  normal  motion  of  the  wet  surface 
of  the  submerged  structure  changes.  For  Ns0,2,3,...,  the  surface  expansion 
functions  are  generated  from  Eq.  (16)  and  are  employed  to  capture  local 
responses.  For  N*l,  the  surface  expansion  functions  are  generated  from 
Eqs.  (13) -(16)  and  are  employed  to  capture  both  global  and  local  responses. 
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Once  the  decision  has  been  made  about  which  of  the  above  sets  of 
convenient  functions  are  to  be  used  to  describe  the  normal  displace¬ 
ment  of  the  structure-fluid  interface  in  a  given  problem,  the  surface 
expansion  functions  may  then  be  generated.  The  ELSHOK  code  accomplishes 
this  task  by  means  of  the  Gram-Schmidt  orthogonalization  process  (Ref.  [5]  ). 
Normalization  of  the  surface  expansion  functions  is  done  according  to  Eq.(6), 
d.  Incident  Loading 

As  may  be  seen  from  Eqs.  (1)  -(8),  the  solution  of  the  governing 
equations  requires  a  knowledge  of  both  the  incident  pressure  and  the 
incident  fluid  particle  velocity  applied  by  the  shock  wave  to  the  wet 
surface  of  the  submerged  structure.  In  ELSHOK,  the  incident  wave  is  assumed 
to  be  a  spherical  wave  propagating  through  a  linear  acoustic  fluid. 

Thus,  the  fluid  particle  velocity  may  be  determined  from  a  knowledge  of 
the  incident  pressure.  Two  representations  of  the  incident  pressure 
are  included  in  ELSHOK  :  a  known  pressure  profile  such  as  that  deter¬ 
mined  by  experiment  and  the  empirical  decaying  exponential  form  of  Ref. 

[6]  for  spherical  charges.  Neither  loading  option  considers  the  effects 
of  cavitation. 

For  incident  loading  specified  in  terms  of  a  known  pressure  profile,  the 
incident  pressure  and  normal  fluid  particle  velocity  at  a  given  point  on  the 


wet  surface  may  be  written  as 


R-S 


R-S 


Pl(R.t)  -  ^  P(t - H(t - f-) 


(17) 


u^(R,t)  *  [^p^(R,t)  f  Pt(R,t)  dT]  cos  (t) 


(18) 


where  R  is  the  distance  from  the  origin  of  the  incident  spherical  wave 

to  the  point  of  interest,  S  is  the  stand-off  distance  (the  distance 

o 

between  the  origin  of  the  incident  wave  and  the  closest  point  on  the  body) , 

H  is  the  Heaviside  unit  function,  P  is  the  known  pressure-time  function 

at  the  point  where  the  shock  wave  first  contacts  the  submerged  body 

(R  =  S  )  ,  and  d)  is  the  angle  between  the  spherical  wave  front  and  the 
o 

inward  normal  to  the  wet  surface  of  the  submerged  body  at  the  point  of 
interest.  It  may  be  observed  that  the  time  t  is  measured  from  the  first 
arrival  of  the  shock  front  at  the  surface  of  the  submerged  structure. 

When  the  decaying  exponential  form  is  employed,  the  incident  pres¬ 
sure  is  given  by 

R-S  R-S 

Pj(R,t)  =  P^(R,W)  exp  [-  (t - /  e^(R,W)]  H(t - f-)  (19) 

in  which  P^(R,W)  is  the  initial  peak  pressure,  6^  (R,W)  is  the  time 
constant  of  exponential  decay,  and  W  is  the  weight  of  the  spherical 
charge  used  to  produce  the  shock  wave.  For  this  loading  option,  the 
incident  normal  fluid  particle  velocity  is  approximated  by  using  Eq.  (19) 
in  Eq.  (18)^ yielding 

p  0  R-S 

^  ^^o  ■  ■  ~r  ^  ^20) 

Such  a  procedure  for  determining  the  incident  fluid  particle  velocity 
is  approximate  because,  as  shown  in  Ref.  [6]*  the  spatial  decay  of  P^ 
differs  slightly  from  the  1/R  decay  in  a  true  spherical  wave. 


18 


III.  ORGANIZATION  AND  OPERATION 


a.  Outline 

The  ELSHOK  computer  code  may  be  used  to  calculate  the  linearly 
elastic  shock  response  of  a  submerged  ring-stiffened  shell  of  revolution 
containing  internal  structure.  A  typical  shell-substructure  configura¬ 
tion  is  displayed  in  Fig.  1.  The  analysis  of  such  a  problem  using  ELSHOK 
is  based  upon  a  modal  substructuring  procedure  which  does  not  involve 
the  modes  and  natural  frequencies  of  the  combined  structural  system  (shell 
and  appendages).  In  addition,  as  discussed  in  Ref.  [3],  a  system  stiffness 
matrix  is  not  required,  since  the  substructuring  technique  employs  inter¬ 
action  forces  and  moments  and  compatibility  of  deformation  at  the  shell- 
substructure  junctions  to  solve  the  dynamic  response  problem.  A  modal 
mass  matrix  of  size  v  ^  by  v  g  is  required  for  those  problems  involving 
concentrated  or  lumped  masses  attached  to  the  shell. 

Separate  modal  analyses  of  the  shell  S  and  any  internal  structure 
a  are  required  when  using  ELSHOK.  The  shell,  a  body  of  revolution,  is 
analyzed  using  the  B0S0R4  finite  difference  code  of  Ref.  [7].  The 
Internal  structure,  if  any,  is  analyzed  by  means  of  the  finite  element 
code  SAPIV  of  Ref.  [8].  User  familiarity  with  these  codes  is  assumed. 
However,  since  some  modifications  to  B0S0R4  and  SAPIV  have  been  made  to 
allow  the  use  of  these  codes  in  the  ELSHOK  program  suite,  some  infor¬ 
mation  about  them  will  be  given  as  needed. 

The  ELSHOK  code  is  written  in  FORTRAN  IV  for  use  on  Control  Data 
Corporation  (CDC)  computers.  The  code  is  modular  in  construction  and 
operates  on  CDC  SCOPE,  KRONOS,  NOS,  and  NOS/BE  operating  systems  of  various 
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levels.  Machine  dependency  is  minimal  and  the  transfer  of  ELSHOK  to  com¬ 
puters  of  other  makes  is  not  a  major  problem. 

The  source  decks  for  the  ELSHOK  code  are  maintained  in  the  format  of 
the  CDC  UPDATE  batch  editing  system  (Ref.  [9]).  Thus,  the  entire  ELSHOK 
code  is  stored  in  a  collection  of  files  knovm  as  "program  libraries". 

Every  card  image  in  a  program  library  is  assigned  a  unique  identifier. 

This  enables  a  user  to  reference  any  card  image  during  an  UPDATE  modi¬ 
fication  run.  Hence,  card  images  may  be  inserted  into  or  deleted  from 
a  program  library  to  alter  a  program  module.  The  modified  program  library 
may  then  be  processed  by  the  FORTRAN  compiler  and  the  resulting  compiled 
or  object  code  may  be  executed.  The  UPDATE  system  allows  a  user  great 
flexibility  in  the  use  of  ELSHOK. 

All  of  the  major  components  of  the  ELSHOK  code,  except  B0S0R4, 
make  use  of  a  technique  usually  referred  to  as  "dynamic  allocation  of 
central  memory"  or  "main  memory  management".  In  this  technique,  only 
those  arrays  required  during  a  program  execution  step  are  stored  in  a 
column  matrix,  called  the  "master  array",  allocated  to  blank  common.  An 
Individual  array  is  located  by  specifying  the  location  of  its  first  word 
within  the  master  array;  the  location  of  its  last  word  is  determined  such 
that  space  is  allocated  for  the  minimum  array  size  needed  for  the  problem 
under  study.  Moreover,  all  arrays  of  substantial  size  are  variably  or 
adjustably  dimensioned  in  subroutines  using  them.  Hence,  the  capacity 
of  a  code  employing  main  memory  management  is  limited  only  by  the  central 
memory  available  to  the  computer  used  and  not  by  the  number  of  modes, 
nodal  points,  etc,  as  when  arrays  are  dimensioned  to  fixed  sizes. 
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The  implementation  of  the  dynamic  allocation  in  ELSHOK  does  not  allow  the 
size  of  the  master  array  in  blank  common  to  be  changed  during  program 
execution.  It  is,  therefore,  necessary  to  dimension  the  master  array  to 
its  full  size  prior  to  compilation.  Thus,  two  steps  are  usually  required 
when  executing  a  given  component  of  ELSHOK.  First,  the  needed  size  of 
the  master  array  is  determined,  most  often  by  means  of  an  execution 
option  available  with  most  of  the  components  of  ELSHOK.  Secondly,  the 
dimension  of  the  master  array  is  set,  the  program  compiled,  and  the 
execution  begun.  Although  a  certain  amount  of  data  transfer  to  and 
from  a  disk  file  (scratch  file)  is  required  for  the  efficient  use  of  the 
available  computer  core  via  main  memory  management,  the  resulting 
Increase  in  the  size  of  a  problem  which  may  be  solved  by  ELSHOK  is  by  far 
the  more  important  consideration. 

The  major  components  of  ELSHOK  are: 

1.  B0S0R4  -  structural  analyzer  for  shell  (Ref.  [7]) 

2.  ACESNID  -  virtual  mass  processor  (Ref.  [10]) 

3.  PIFLASH  -  shell-fluid  processor 

4.  SAPIV  -  structural  analyzer  for  substructure  (Ref.  [8]) 

5.  PICRUST  -  substructure  processor 

6.  USLOB  -  time  integration  processor 

7.  PUSLOB  -  plotting  processor. 

Several  other  minor  components  will  be  alluded  to  at  the  appropriate 
times.  Each  component  of  the  ELSHOK  code  has  been  written  such  that  any 
consistent  set  of  units  may  be  employed  throughout  the  analysis  of  a 
given  problem.  Units  may  be  changed  when  preparing  plots  of  velocity¬ 
time  histories  and  shock  spectra. 


The  use  of  the  above  program  components  in  the  four  phases  of  an 
ELSHOK  calculation  is  Illustrated  in  Fig.  2.  A  brief  description  of  each 
phase  follows. 

b.  Phase  I  -  Shell  and  Fluid  Data 
1.  Structural  Analysis  of  Shell 

The  first  step  in  Phase  1  is  to  model  the  rotationally  symmetrical 
portion  of  the  structure  being  analyzed  (Fig.  1)  for  use  with  the  B0S0R4 
computer  code.  The  resulting  idealized  mathematical  model  may  contain 
bulkheads  and  rings,  but  must  exclude  internal  appendages  and  concentrated 
masses.  The  B0S0R4  code  is  quite  general  and  applies  to  segmented,  ring- 
stiffened,  branched  shells  of  revolution  having  various  meridional  geom¬ 
etries,  wall  constructions,  and  ring  reinforcements.  Circumferential 
rings  may  be  treated  as  discrete  stiffeners  or  represented  in  the 
orthotropic  approximation  (smeared)  when  closely  spaced.  Stringers  must 
be  treated  in  the  orthotropic  approximation  to  preserve  structural 
axisymmetry.  Symmetry  about  the  midpoint  of  the  length  may  be  accounted 
for  by  appropriate  boundary  conditions  at  the  plane  of  symmetry.  In 
ELSHOK,  the  mathematical  model  of  the  entire  main  body  S  (or  half,  if 
symmetry  is  employed)  is  termed  the  "full  model". 

In  problems  where  only  a  portion  of  a  structure  is  of  interest, 
as  is  the  case  when  a  localized  loading  is  applied  to  a  section  of  a 
structure  enclosed  by  stiff  bulkheads,  the  user  may  highlight  the  region 
of  Interest.  This  is  accomplished  in  ELSHOK  by  employing  two  different 
mathematical  models  of  the  structure  when  using  B0S0R4:  the  full  model 
described  above  and  the  model  of  the  region  of  interest,  called  the 


compartment  model  .  The  full  model  is  used  to  capture  the  gross  effects 
of  the  response,  such  as  the  rigid  body  translation  and  whipping  of  the 
entire  structure  (N  =  1)  and  the  rolling  and  twisting  of  the  entire 
structure  (N  =  0  torsional) .  The  compartment  model  is  used  to  capture 
the  local  details  of  the  response  in  the  region  of  interest  (N  =  0 
breathing,  N  > 2) .  The  use  of  this  technique  is  optional,  and  the  full 
model  alone  may  be  employed  for  the  analysis  of  the  submerged  shock 
response.  More  information  about  the  use  of  this  approximation  may  be 
found  in  Appendix  A. 

Once  the  mathematical  representation  of  the  main  body  S  has  been  deci¬ 
ded  upon,  the  in-vacuo  free-free  modes  and  natural  frequencies  of  S  must  be 
determined  using  the  B0S0R4  code.  For  compatibility  with  ELSHOK,  a 
separate  B0S0R4  calculation  must  be  performed  for  each  circumferential 
harmonic  N  included  in  the  analysis.  When  only  a  full  model  is  employed, 
the  breathing  modes  and  any  torsional  modes  required  (N  =  0)  may  be 
determined  in  one  calculation  or  in  separate  calculations,  as  desired. 

For  each  B0S0R4  calculation,  the  resulting  shell  geometry,  material 
properties,  generalized  masses,  modes,  and  natural  frequencies  are 
written  on  a  file  (disk  or  magnetic  tape)  called  a  "shell  mode  file". 
Operating  instructions  for  the  B0S0R4  code  may  be  found  in  Appendix  A. 

2.  Calculation  of  Virtual  Mass 

The  second  step  in  Phase  I  is  the  computation  of  the  virtual  mass 
array,  which  provides  the  late-time  contribution  of  the  DAA.  In  the 
ELSHOK  suite  of  programs,  this  is  accomplished  by  means  of  the  ACESNID 
(ACCESSION  TO  ^ERTIA  AND  DAMPING)  code  of  Ref.  [10].  In  ACESNID,  the  virtual 
mass  array  is  determined  from  the  solution,  based  on  simple  sources,  of  a 


low-frequency  steady-state  (incompressible)  problem  in  which  normal  dis¬ 
placements  corresponding  to  surface  expansion  functions  are  applied  to  the 
surface  of  a  cavity  of  revolution  in  the  infinite  fluid  having  the  same 
shape  and  size  as  the  wet  surface  of  the  structure  under  consideration. 

Program  ACESNID  obtains  the  geometry  of  the  cavity  in  the  fluid  from 
the  geometry  supplied  on  a  shell  mode  file  for  a  full  model.  The  sur¬ 
face  expansion  functions  required  are  computed  via  the  Gram-Schmidt 
orthogonallzatlon  process  (Ref.  [5])»  as  discussed  earlier.  The  symmetrical 
virtual  mass  array  for  each  circumferential  wave  number  n  may  be  determined 
in  separate  calculations,  or  all  of  the  required  values  of  N  may  be 
considered  in  one  calculation,  depending  on  the  application.  The  virtual 
mass  and  data  needed  to  re-compute  the  surface  expansion  functions  when 
required  elsewhere  are  written  on  a  file,  called  the  "virtual  mass  file", 
containing  virtual  mass  data  for  all  values  of  N  used.  This  file  may  be 
prepared  by  one  execution  of  ACESNID  for  all  values  of  N  considered  or  by 
merging  the  results  from  separate  ACESNID  calculations.  Appendix  B  con¬ 
tains  operating  instructions  for  the  ACESNID  code. 

3.  Completion  of  Data  for  Shell  and  Fluid 

The  third  and  final  step  in  Phase  I  is  to  complete  the  data  for 
describing  the  shell  and  the  fluid.  For  this  purpose,  the  PIFLASH  (PREPARE 
INPUT  FOR  FLUID  AND  SHELL)  computer  code  is  provided.  This  code  takes  the 
data  from  all  of  the  shell  mode  files  and  from  the  virtual  mass  file  and 
reorganizes  it  to  facilitate  the  solution  of  the  response  equations. 

During  the  execution  of  PIFLASH,  the  modes  obtained  from  B0S0R4  are  all 
normalized  to  the  total  mass  of  the  entire  main  body  S.  Thus,  when  symmetry 
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is  employed,  the  mass  of  the  mathematical  model  is  doubled.  User  supplied 
input  determines  how  many  of  the  shell  modes  available  are  to  be  provided 
for  the  time-history  calculation.  Any  concentrated  masses  included  in  the 
mathematical  model  of  the  structure  under  study  are  accounted  for  during 
the  execution  of  PIFLASH.  The  PIFLASH  code  constructs  a  file,  called  the 
"shell-fluid  file",  containing  all  of  the  information  needed  to  describe 
the  shell,  any  concentrated  masses  attached  to  the  shell,  and  the  fluid 
during  the  subsequent  analysis.  If  no  internal  appendages  or  substructures 
are  included  in  the  submerged  structure  being  studied,  the  shell-fluid 
file  serves  as  the"input  file"  for  the  time-history  calculation  described 
later.  Detailed  information  about  the  use  of  PIFLASH  may  be  found  in 
Appendix  C. 

c.  Phase  II  -  Substructure  Data 
1.  Structural  Analysis  of  Internal  Appendages 

The  first  stage  of  Phase  II  consists  of  the  development  of  a  finite 
element  model  of  each  piece  of  internal  structure  a.  The  SAPIV  computer 
program  is  included  in  the  ELSHOK  suite  of  computer  codes  for  this  pur¬ 
pose.  Structural  systems  comprised  of  a  number  of  different  types  of 
finite  elements  may  be  analyzed  using  SAPIV.  The  SAPIV  library  of 
finite  elements  contains  a  three-dimensional  truss  element,  a  three- 
dimensional  beam  element,  a  plane  stress  and  plane  strain  element,  a  two- 
dimensional  axisymmetric  solid,  a  three-dimensional  solid,  a  thick  shell 
element,  a  thin  plate  or  thin  shell  element,  a  boundary  element,  and  a 
pipe  element.  Any  concentrated  masses  attached  to  the  substructures  must 
be  accounted  for  in  the  mathematical  models  used  with  SAPIV.  If  no  sub¬ 
structures  are  present  in  a  given  investigation.  Phase  II  does  not  apply 


and  must  be  skipped. 


Once  the  mathematical  model  of  each  substructure  has  been  completed, 
the  fixed-base  modes  and  natural  frequencies  of  each  substructure  must  be 
determined.  A  separate  SAPIV  calculation  must  be  performed  for  each 
internal  appendage.  For  each  such  calculation,  the  resulting  geometry, 
stiffness  and  mass  information,  modes,  and  natural  frequencies  are 
written  on  a  file  called  a  "substructure  mode  file".  The  manner  in  which 
SAPIV  must  be  used  for  an  ELSHOK  calculation  is  described  in  APPENDIX  D. 


2 .  Completion  of  Data  for  Substructures 

The  second  and  final  step  in  Phase  II  is  to  complete  the  data 
required  for  each  and  every  internal  appendage.  Program  PICRUST 
(PREPARE  INPUT  FOR  CALCULATING  RESPONSE  OF  SUBSTRUCTURE)  accomplishes 
this  task.  The  PICRUST  code  takes  the  data  from  a  given  substructure 
mode  file,  rearranges  it,  and  computes  a  number  of  additional  items 
needed  to  perform  the  solution  of  the  response  equations.  A  separate 
calculation  is  required  for  each  substructure. 

For  a  given  substructure,  the  PICRUST  code  allows  the  user  to 
determine  how  many  of  the  fixed-base  modes  available  are  to  be  provided 
for  the  time-history  calculation.  The  user  also  specifies  the  con¬ 
nectivity  between  the  shell  and  the  substructure  under  consideration. 
During  an  execution  of  PICRUST,  influence  coefficients  corresponding  to 
the  base  or  support  motions  of  a  substructure  (the  of  Eq.  (10)  or 
the  constraint  modes  of  Ref.  [3])are  computed.  Modal  coefficients  for 
the  calculation  of  interaction  forces  and  moments  between  the  shell  and 
a  substructure  are  also  determined.  For  each  internal  appendage  studied. 


the  results  of  the  PICRUST  execution  are  stored  on  a  file  called  a  "sub¬ 


structure  file".  The  substructure  files  for  all  substructures  in  a  given 


investigation  must  be  merged  with  the  shell-fluid  file  for  use  in  the  sub¬ 


sequent  shock  response  problem,  as  discussed  under  Phase  III,  Operating 
instructions  for  the  PICRUST  code  are  given  in  Appendix  E. 
d.  Phase  III  -  Submerged  Shock  Response 

1.  Completion  of  Input  File 

Almost  all  of  the  input  data  required  for  the  structure-fluid 
interaction  calculation  must  be  supplied  in  a  single  file,  termed  the 
"input  file".  For  problems  without  internal  appendages,  the  shell-fluid 
file  described  under  Phase  I  is  the  input  file  for  the  shock  response 
calculation.  For  problems  involving  substructures,  the  shell-fluid  file 
must  be  merged  with  all  of  the  substructure  files  of  Phase  II.  This  may 
be  easily  accomplished  by  means  of  simple  control  card  (job  control 
language)  operations,  or  a  FORTRAN  computer  program  may  be  written  for 
this  purpose.  The  usual  ELSHOK  analysis  employs  the  control  card  operations 
as  explained  in  Appendix  F. 

2.  Time  Integration 

The  ELSHOK  suite  of  computer  programs  includes  the  computer  code  USLOB 
(UNDERWATER  SHOCK  LOADING  OF  BODIES)  for  solving  the  governing  equations 
of  the  current  formulation  of  the  structure-fluid  interaction  problem. 

The  USLOB  time  integration  processor  employs  the  Runge-Kutta  technique, 
in  the  modification  due  to  Gill  (Ref.  [11]), to  perform  the  integration  in 
time.  In  addition  to  the  input  file  described  above,  the  user  must  also 
supply  time  step  information  and  a  description  of  the  shock  wave  loading. 

At  the  discretion  of  the  user,  input  may  be  supplied  to  cause  the  USLOB 
code  to  provide  printed  velocity-time  histories  for  points  on  the  shell 
and  any  internal  appendages. 


To  permit  the  production  of  plotted  velocity-time  histories  in  Phase  IV 
of  an  ELSHOK  calculation,  program  USLOB  provides  up  to  two  post-processing 
files.  The  first  file,  called  the  "generalized  velocity  file",  contains 
the  generalized  velocities  of  the  shell,  (t)  ,  and  any  substructures, 
q^j(t),  and  is  always  produced  during  an  execution  of  USLOB.  The  second 
file,  called  the  "shell-substructure  transformation  file",  is  produced 
only  when  internal  appendages  are  attached  to  the  shell.  This  file  con¬ 
tains  the  shell-substructure  transformation  matrices  [the  ^  of  Eq.  (9)] 
needed  to  describe  the  motion  of  the  shell,  at  the  shell-substructure  con¬ 
nection  points,  with  reference  to  the  coordinate  systems  used  for  the  sub¬ 
structures.  Further  details  of  the  USLOB  code  may  be  found  in  Appendix  F. 
e.  Phase  IV  -  Plots  of  Velocity-Time  Histories 


1.  General  Information 

After  the  completion  of  Phases  I-III  of  an  ELSHOK  calculation,  it  is 
usually  necessary  to  present  the  computed  results  in  the  form  of  plotted 
velocity-time  histories.  The  PUSLOB  (PLOTS  FOR  USLOB)  computer  code  is 
provided  for  this  purpose.  The  input  file,  the  generalized  velocity  file, 
and  the  shell-substructure  transformation  file,  if  applicable,  must  be 
available  to  complete  Phase  IV.  All  plotting  using  PUSLOB  is  performed 
on  a  TEKTRONIX  graphics  terminal  (AGIO,  4051,  4052,  e.g.)in  conjunction  with 
a  TEKTRONIX  hard  copy  unit  (4631,  e.g.).  User  provided  input  selects  from 
among  the  various  output  options,  specifies  labeling,  adjusts  scales,  and 
the  like. 

2 .  Summary  of  Output  Options 

Using  PUSLOB,  plots  of  velocity-time  histories  may  be  produced  for 


any  point  on  the  shell  S  and  for  anv  nodal  point  on  any  internal  appendage  a. 
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At  the  discretion  of  the  user,  velocity  plots  for  the  shell  may  be  pro¬ 
vided  in  the  local  coordinate  system  and/or  in  the  global  coordinate  system. 
In  the  local  system,  the  inward  normal  (n) ,  meridional  (s) ,  and  circum¬ 
ferential  (9)  directions  of  Fig.  1  are  employed.  In  the  global  system, 
the  athwartship  (-  Z ) ,  fore-aft  (+X) ,  and  downward  (-Y)  directions  of  Fig. 

1  are  used.  All  plotting  for  a  substructure  is  performed  with  reference 
to  the  coordinate  system  of  the  given  substructure.  In  addition  to  plotted 
results,  coded  card-image  time  histories  may  be  prepared  for  transmittal  to 
other  computer  installations  or  for  additional  post-processing. 

If  pseudo-velocity  shock  spectra  are  required,  these  may  be  pro¬ 
duced  from  the  card-image  file  mentioned  above.  At  present,  shock  spectra 
are  not  produced  by  the  PUSLOB  code,  but  by  a  separate  small  computer 
program  of  the  type  available  at  mf'-t  scientific  establishments.  Software 
for  a  TEKTRONIX  4051  or  4052  is  avai  able  with  ELSHOK  to  produce  pseudo¬ 
velocity  shock  spectra  from  velocity-time  histories,  after  processing  the 
time  histories  on  a  digital  tablet  (digitizer).  The  inclusion  of  the 
spectra  capability  in  PUSLOB  is  a  modification  planned  for  the  future. 
Appendix  G  gives  more  information  about  the  use  of  the  PUSLOB  computer 


IV.  SAMPLE  PROBLEMS 
a.  Sample  Problem  1 

Consider  the  structure  shown  in  Fig.  3  to  be  immersed,  initially  at 

rest,  in  an  infinite  expanse  of  water  and  to  be  excited  by  side-on  loading 

in  the  form  of  a  spherical  step  wave  of  constant  strength  P=10  psi  (68.947 
2 

kU/m  ) .  The  acoustic  wave  is  assumed  to  emanate  from  a  point  located 
opposite  the  longitudinal  midpoint  of  the  structure  at  a  stand-off  distance 
of  69.9  ft  (21.305  m) .  Moveover,  the  incident  wave  is  assumed  to  impinge 
upon  the  submerged  structure  so  as  to  excite  doubly  symmetrical  motion 
(with  respect  to  the  incident  wave).  The  orientation  of  the  shock  wave 
may  be  seen  in  Fig.  3.  This  sample  problem  involves  only  a  full  model 
of  the  shell  S  and  illustrates  the  use  of  the  straightforward  modeling 
technique  discussed  in  Appendix  A.  One  substructure  a  is  included  in 
the  problem,  namely,  the  long  beam  and  its  support  brackets  (Details 
1  and  2  of  Fig.  3).  The  bars  and  mounting  plates  (Details  3  and  4, 
respectively,  of  Fig.  3)  are  regarded  as  concentrated  masses  rigidly 
attached  to  the  shell  S. 

The  symmetrical  in-vacuo  mode  shapes  and  corresponding  natui al 
frequencies  of  the  empty  shell  of  revolution  S  were  obtained  by  applying 
the  B0S0R4  computer  program  to  the  two-segment  mathematical  model 
pictured  in  Fig.  4.  For  this  computation,  the  numerous  light  stiffeners 
were  treated  in  the  orthotropic  approximation,  while  the  heavy 
stiffeners  were  treated  as  discrete  rings.  The  end  closure  of  Fig.  4 
was  modeled  as  a  flat  plate.  Only  families  of  modes  having  circum¬ 
ferential  wave  numbers  N=0,l,2,3  were  considered,  and  30  modes  were 
computed  for  each  of  these  values  of  N.  Less  than  3  min  of  computer  time 


(central  processor  time)  was  used  on  a  CDC  6600  for  these  calculations. 
Elimination  of  the  torsional  and  rolling  modes  for  N=0  reduced  the  num¬ 
ber  of  modes  of  S  employed  in  the  solution  to  115.  Appendix  A  contains 
a  listing  of  the  input  data  prepared  for  the  B0S0R4  code  for  Sample 
Problem  1. 

The  virtual  mass  applicable  to  doubly  symmetrical  response  of 
the  submerged  structure  of  Fig.  3  was  computed  using  the  ACESNID  com¬ 
puter  code.  For  the  current  problem,  the  cavity  in  the  fluid  was 
represented  by  12  annular  source  bands  on  each  end  closure  and  75  source 
bands  on  the  cylindrical  section.  In  performing  the  calculation  of  the 
virtual  mass,  the  speed  of  sound  in  the  fluid  was  taken  to  be  5000  ft /sec 
(1524  m/sec),  corresponding  to  a  diameter  transit  time  of  0.56  msec.  For 
each  circumferential  harmonic  N  considered,  the  normal  motion  of  the  wet 
surface  of  the  shell  S  was  described  by  15  surface  expansion  functions 
corresponding  to  Eq.  (13)  for  the  cylindrical  section  and  3  surface 
expansion  functions  corresponding  to  Eqs.  (14)  and  (15)  for  the  end  closure 
Thus,  72  surface  expansion  functions  in  all  were  employed  in  the  solution. 
About  8  min  of  central  processor  time  was  required  on  a  CDC  6600  to 
determine  the  virtual  mass.  The  input  data  for  the  virtual  mass  calcu¬ 
lation  of  Sample  Problem  1  may  be  found  in  Appendix  B. 

The  symmetrical  in-vacuo  mode  shapes  and  corresponding  natural  fre¬ 
quencies  of  the  substructure  a  were  obtained  by  applying  the  SAPIV 
computer  code  to  the  assemblage  of  finite  elements  displayed  in  Fig.  5. 

For  this  calculation,  24  beam  elements  were  used  to  connect  25  nodal  points 
and  the  substructure  a  was  constrained  to  move  in  the  X  -  Z  (horizontal) 
plane  because  of  the  doubly  symmetrical  nature  of  the  problem  under  study. 
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It  should  be  pointed  out  that  an  element  release  was  employed  to  pre¬ 


vent  the  development  of  a  bending  moment  in  the  support  bracket  at  the 
junction  of  the  beam  and  bracket  and  that  the  weight  of  the  support  bracket 
was  neglected.  The  first  5  symmetrical  fixed-based  modes  were  computed 
and  employed  in  the  subsequent  shock  response  calculation.  About  3  sec 
of  CDC  6600  computer  time  was  required  for  the  execution  of  the  SAPIV  code. 
The  set  of  input  data  used  for  this  calculation  is  listed  in  Appendix  D. 

The  shell-fluid  file  (Appendix  C)  was  produced  by  means  of  the  PIFLASH 
computer  program.  Less  than  30  sec  of  central  processor  time  on  a  CDC 
6600  was  needed  for  this  purpose.  The  set  of  input  data  prepared  for  PIFLASH 
for  Sample  Problem  1  may  be  found  in  Appendix  C.  Processing  to  produce 
the  substructure  file  (Appendix  E)  for  Sample  Problem  1  was  accomplished 
by  means  of  the  PICRUST  computer  code  in  less  than  15  sec  of  CDC  6600  com¬ 
puter  time.  Appendix  E  contains  the  input  data  employed  for  this  execu¬ 
tion  of  PICRUST.  The  input  file  for  the  response  calculation  of  Sample 
Problem  1  was  constructed  from  the  shell-fluid  and  substructure  files 
above  using  the  control  statement  procedure  discussed  in  Appendix  F.  This 
procedure  required  less  than  2  sec  of  central  processor  time  on  a  CDC  6600. 

Solution  of  the  transient  structure-fluid  interaction  problem 
described  earlier  was  performed  by  means  of  the  USLOB  code.  A  constant 
time  step  of  0.000025  sec  (about  1/25  of  a  diameter  transit  time)  was  used, 
and  the  calculation  was  continued  for  more  than  8  diameter  transit  times. 
Nearly  8  min  of  CDC  6600  central  processor  time  was  required  to  com¬ 
pute  the  shock  response  of  the  submerged  structure.  Appendix  E  con¬ 
tains  the  input  data  employed  for  this  calculation.  Plotting  of  the 
results  was  accomplished  by  means  of  the  PUSLOB  computer  code,  about  15 
sec  of  CDC  6600  computer  time  being  required  for  the  problem  at  hand. 


The  input  data  employed  to  produce  the  desired  velocity-time  plots  may  be 

found  in  Appendix  G.  Figures  6-13  display  the  results  of  the  calculation 

*) 

performed  as  Sample  Problem  1.  These  figures  are  embellished  versions  of 
those  produced  by  the  PUSLOB  code, 
b.  Sample  Problem  2 

Consider  the  structure  shown  in  Fig. 3  to  be  immersed,  initially  at 
rest,  in  an  infinite  expanse  of  water  and  to  be  excited  by  side-on  loading 
in  the  form  of  a  spherical  wave  of  decaying  exponential  profile  in  time. 

The  shock  wave  is  assumed  to  be  produced  by  the  underwater  explosion 
of  a  3  lb  (1.36  kg)  spherical  charge  of  TNT  at  a  point  located  opposite 
the  longitudinal  midpoint  of  the  structure  at  a  stand-off  distance  of 
2.802  ft  (0.854  m),or  one  model  diameter.  As  in  Sample  Problem  1,  the 
incident  wave  is  assumed  to  impinge  upon  the  submerged  structure  so  as 
to  excite  doubly  symmetrical  motion.  The  orientation  of  the  shock  wave 
may  be  seen  in  Fig.  3.  This  sample  problem  is  presented  to  demonstrate 
the  use  of  the  approximate  modeling  technique,  involving  a  region  of 
interest,  discussed  in  Appendix  A.  One  substructure  o  is  included  in 
the  problem,  namely,  the  long  beam  and  its  support  brackets  (Details 
1  and  2  of  Fig.  3).  The  bars  and  mounting  plates  (Details  3  and  4, 
respectively,  of  Fig.  3)  are  regarded  as  concentrated  masses  rigidly 
attached  to  the  shell  S. 

The  symmetrical  In-vacuo  mode  shapes  and  corresponding  natural 
frequencies  of  the  empty  shell  of  revolution  S  were  obtained  by  applying 
the  B0S0R4  computer  code  to  the  three-segment  full  model  and  to  the  one  - 
segment  compartment  model  pictured  in  Fig.  14.  As  in  Sample  Problem  1, 

*)  The  term  "PROBONE"  (from  PROBLEM  ONE)  appearing  in  the  alphameric  headings 
of  Figs.  6-13  is  the  name  of  the  input  file  used  for  Sample  Problem  1. 


the  numerous  light  stiffeners  were  treated  in  the  orthotropic  approxi¬ 
mation,  while  the  heavy  stiffeners  were  treated  as  discrete  rings.  The 
end  closure  of  Fig.  14  was  modeled  as  a  flat  plate.  Only  families  of 
modes  having  circumferential  wave  numbers  N=0, 1,2,3  were  considered  in 
the  modal  analysis  of  the  shell  S.  For  N=l,  30  modes  were  computed  using 
the  full  model  of  S.  For  each  value  of  N=0,  2,3,  20  modes  were  determined 
using  the  compartment  model  of  S.  Less  than  2  min  of  computer  time 
(central  processor  time)  was  used  on  a  CDC  6600  for  these  calculations. 

The  use  of  15  breathing  modes  for  N=0,  15  modes  for  N=2 ,  15  modes  for  N=3, 
and  30  modes  for  N=1  brought  the  total  number  of  modes  of  S  employed  in 
the  transient  response  problem  to  75.  Appendix  A  contains  the  input  data 
prepared  for  the  B0S0R4  code  for  Sample  Problem  2. 

The  virtual  mass  suitable  for  use  with  the  approximate  modeling 
technique  employed  for  the  shell  S  was  computed  using  the  ACESNID  com¬ 
puter  code.  The  cavity  in  the  fluid  was  represented  by  10  annular  source 
bands  on  each  end  closure  and  75  source  bands  on  the  cylindrical  section. 
In  performing  the  calculation  of  the  virtual  mass,  the  speed  of  sound 
in  the  fluid  was  taken  to  be  5000  ft/sec  (1524  m/sec) .corresponding  to  a 
diameter  transit  time  of  0.56  msec.  For  N=l,  the  normal  motion  of  the 
wet  surface  of  the  shell  S  was  described  by  8  surface  expansion  functions 
corresponding  to  Eq.  (13)  for  the  cylindrical  section,  3  surface  expansion 
functions  corresponding  to  Eqs.  (14)  and  (15)  for  the  end  closures,  and 
7  surface  expansion  functions  corresponding  to  Eq.  (16)  for  the  compart¬ 
ment  region.  For  each  value  of  N=0,2,3,  7  surface  expansion  functions 
corresponding  to  Eq.  (16)  were  used  for  the  compartment  region.  Thus,  39 
surface  expansion  functions  in  all  were  employed  in  the  solution.  About 
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7  min  of  central  processor  time  was  required  on  a  CDC  6600  to  determine 
the  virutal  mass.  The  input  data  for  the  virtual  mass  calculation  of 
Sample  Problem  2  may  be  found  in  Appendix  B. 

The  symmetrical  in-vacuo  fixed-base  modes  and  corresponding  natural 
frequencies  of  the  substructure  o  were  obtained  by  applying  the  SAPIV 
computer  code  to  the  assemblage  of  finite  elements  displayed  in  Fig.  5. 

Since  the  identical  mathematical  model  was  used  for  Sample  Problem  1, 
the  discussion  of  the  modeling  and  computational  details  will  not  be 
repeated  here.  It  should  be  pointed  out  that,  for  Sample  Problem  2,  the  shell 
cross  section  in  the  note  of  Fig.  5  (b)  must  be  changed  to  mesh  point  1 
of  segment  3.  The  set  of  input  data  used  for  the  modal  analysis  of  the 
substructure  a  for  Sample  Problem  2  is  listed  in  Appendix  D. 

The  PIFLASH  computer  program  was  utilized  to  produce  the  shell- 
fluid  file.  Less  than  15  sec  of  central  processor  time  on  a  CDC  6600 
was  needed  for  this  purpose.  The  set  of  input  data  prepared  for  PIFLASH 
for  Sample  Problem  2  may  be  found  in  Appendix  C.  Processing  to  produce 
the  substructure  file  for  Sample  Problem  2  was  accomplished  by  means  of 
the  PICRUST  computer  code  in  less  than  15  sec  of  CDC  6600  computer  time. 
Appendix  E  contains  the  input  data  employed  for  this  execution  of  PICRUST. 

The  input  file  for  the  response  calculation  of  Sample  Problem  2  was  con¬ 
structed  from  the  shell-fluid  and  substructure  files  above  using  the  con¬ 
trol  statement  procedure  discussed  in  Appendix  F.  This  procedure  required 
less  than  2  sec  of  central  processor  time  on  a  CDC  6600. 

Solution  of  the  transient  structure-fluid  interaction  problem 
described  earlier  was  performed  by  means  of  the  USLOB  code.  A  constant 
time  step  of  0.000025  sec  (about  1/25  of  a  diameter  transit  time)  was  used. 
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and  the  calculation  was  continued  for  more  than  8  diameter  transit  times. 


Nearly  8  min  of  CDC  6600  central  processor  time  was  required  to  compute 

the  shock  response  of  the  submerged  structure.  Appendix  E  contains  the 

input  data  employed  for  this  calculation.  Plotting  of  the  results  was 

accomplished  by  means  of  the  PUSLOB  computer  code,  about  10  sec  of  CDC 

6600  computer  time  being  required  for  the  problem  at  hand.  The  input 

data  employed  to  produce  the  desired  velocity-time  plots  may  be  found 

in  Appendix  G.  Figures  15-19  display  the  results  of  the  calculation 

*) 

performed  as  Sample  Problem  2.  These  figures  are  embellished  versions  of 
those  produced  by  the  PUSLOB  code. 


*)  The  term  "PROBTWO"  (from  PROBLEM  TOO)  appearing  in  the  alphameric  headings 
of  Figs.  15-19  is  the  name  of  the  input  file  used  for  Sample  Problem  2. 
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TYPICAL  SHELL-SUBSTRUCTURE  CONFIGURATION 


FIG.  2  ORGANIZATION  OF  THE  ELSHOK  CODE 
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FIG.  3  STRUCTURAL  MODEL  --  SAMPLE  PROBLEM 


(A)  GEOMETRY 


NOTE:  THE  B0S0R4  CODE  EMPLOYS  TWO  MORE  MESH  POINTS  IN  EACH 
SEGMENT  THAN  THE  NUMBER  PROVIDED  AS  INPUT. 


(B)  DISTRIBUTION  OF  MESH  POINTS  ALONG 
REFERENCE  MERIDIAN 


FIG.  4 


SHELL  OF  REVOLUTION 
SAMPLE  PROBLEM  1 


NOTE:  ACTIVE  PHYSICAL  DEGREES  OF  FREEDOM  ARE  2  -  COMPONENTS 
OF  TRANSLATION  AND  y- COMPONENTS  OF  ROTATION. 


(A)  FINITE  ELEMENT  MODEL  OF  SUBSTRUCTURE 


NOTE:  SHELL  CROSS  SECTION  IS  FOR  MESH  POINT  42  OF  SEGMENT  2. 

(B)  SHELL-SUBSTRUCTURE  CONNECTIVITY 

FIG.  6  INTERNAL  STRUCTURE 
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ATHWARTSHIP  VELOCITIES  AT  STATION  B 
(SEGMENT  2.  POINT  12) 
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ATHWARTSHIP  VELOCITIES  AT  STATION  C 
(SEGMENT  2,  POINT  22) 
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FIG.  13  HORIZONTAL  (z)  VELOCITY  AT 
BEAM  STATION  H  (NODE  25) 

FOR  SAMPLE  PROBLEM  1 
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(A)  GEOMETRY  OF  FULL  MODEL 


NOTE:  THE  B0S0R4  CODE  EMPLOYS  TWO  MORE  MESH  POINTS  IN  EACH 
SEGMENT  THAN  THE  NUMBER  PROVIDED  AS  INPUT. 

(B)  DISTRIBUTION  OF  MESH  POINTS  ALONG  REFERENCE 
MERIDIAN  OF  FULL  MODEL 
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(C)  GEOMtTRY  AND  DISTRIBUTION  OF  MESH  POINTS 
USED  FOR  COMPARTMENT  MODEL 
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FIG.  16  ATHWARTSHIP  VELOCITIES  AT  STATION  D 

(SEGMENT  3.  POINT  1  ) 
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FIG.  16  ATHWARTSHIP  VELOCITIES  AT  STATION  E 

(SEGMENT  3,  POINT  19) 
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1.  Background 

The  B0S0R4  computer  code  of  Ref.  [7]  is  included  in  the  ELSHOK  suite 
of  computer  programs  for  the  modal  analysis  of  the  main  body  S.  This 
code  is  quite  general  and  applies  to  segmented,  ring-stiffened,  branched 
shells  of  revolution  having  various  meridional  geometries,  wall  construc¬ 
tions,  and  ring  reinforcements.  Comprehensive  documentation  of  B0S0R4 
may  be  found  in  Ref.  [7].  A  knowledge  of  the  information  given  in  Sections 
II  and  III  of  the  present  report  will  aid  in  the  understanding  of  this 
appendix. 

In  the  formulation  employed  xn  B0S0R4,  the  independent  variables 
of  the  analysis  are  the  meridional  arc  length  s,  measured  along  the  shell 
reference  meridian  and  the  circumferential  coordinate  angle  G  (Fig.  1). 

The  dependent  variables  are  the  meridional  (u) ,  circumferential  (v) ,  and 
normal  (w)  components  of  displacement  of  the  shell  reference  surface.  The 
numerical  analysis  in  B0S0R4  is  based  upon  a  finite  difference  energy 
method,  so  that  in  the  computer  program,  the  dependent  variables  are  the 
displacement  components  at  discrete  points  (mesh  points )  on  the  shell 
reference  surface.  Further,  the  circumferential  coordinate  9  is  eliminated 
by  means  of  separation  of  variables.  Thus,  for  instance,  the  0-dependence  of 
the  normal  displacement  is  eliminated  by  the  assumption  w(s,0)  =  w^^  cos  N0  or 
w(s,0)  »  Wjj  sin  N0,  as  appropriate,  where  N  is  the  number  of  circumferential 


waves  or  circumferential  harmonics.  The  B0S0R4  code  allows  all  required 
values  of  N  to  be  treated  in  one  execution  of  the  program  or  in  a  series  of 
separate  calculations. 


2.  Structural  Models 


In  most  applications  of  ELSHOK,  the  modeling  of  the  shell  S  for  use 
with  B0S0R4  proceeds  in  the  straightforward  manner  described  in  Ref.  [7]. 
Thus,  the  user  identifies  each  segment  of  the  model,  includes  the 
appropriate  ring  stiffeners,  provides  proper  boundary  conditions  and  con¬ 
nectivity  information,  specifies  physical  properties,  and  so  forth.  The 
resulting  model,  referred  to  as  the  "full  model"  in  what  follows,  is  then 
utilized  in  the  modal  analysis  of  S  for  all  of  the  circumferential  harmonics 
N  Included  in  the  analysis.  In  some  cases,  however,  the  user  may  choose 
to  highlight  a  portion  of  the  shell  S  in  order  to  reduce  the  size  of  the 
computational  effort  at  the  expense  of  some  details  of  deformation  beyond 
the  region  of  Interest.  This  type  of  analysis  may  also  be  performed  using 
ELSHOK. 

When  a  localized  loading  is  applied  to  a  portion  of  the  shell 
enclosed  by  stiff  bulkheads  or  isolated  by  heavy  rings,  the  loading  is 
primarily  going  to  cause  deformation  near  the  point  of  impingenent  of  the 
load  on  the  shell.  Any  internal  structure  in  the  loaded  region  bounded 
by  the  stiff  bulkheads  or  heavy  rings  will  also  be  affected  strongly. 

This  pattern  of  deformation  allows  simplifications  to  be  made  regarding 
the  portion  of  the  structure  studied  carefully.  Specifically,  the  iso¬ 
lated  region  under  localized  loading  may  be  selected  as  the  region  of 
Interest  within  the  full  model.  This  is  accomplished  in  ELSHOK  by  employing 
two  different  mathematical  models  of  the  shell  S  when  using  B0S0R4:  the 
full  model  described  above  and  the  region  of  Interest,  called  the  "com¬ 
partment  model".  A  description  of  the  use  of  these  models  follows. 
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Modes  having  a  circumferential  wave  number  N=1  correspond  to  rigid  body 
displacements  and  whipping  (beam  bending)  of  the  entire  structure.  In  order 
to  permit  the  development  of  these  important  motions,  the  N=1  shell  modes  are 
determined  using  the  full  model.  The  N=0  rolling  and  torsional  modes,  those 
modes  involving  only  circumferential  motion  of  the  structure,  are  also 
determined  using  the  full  model. 

Deformations  or  motions  beyond  the  bulkheads  or  rings  isolating  the 
region  of  Interest  should  be  unimportant  for  the  N=0  axisymmetrical  or 
breathing  modes  (normal  and  meridional  displacements  only)  and  for  the 
modes  corresponding  to  N>2  (normal,  meridional,  and  circumferential  dis¬ 
placements)  .  This  is  so  because  each  bounding  bulkhead  or  heavy  ring  is 
very  stiff  in  its  own  plane  and  moves  primarily  in  that  plane  for  localized 
loading.  Thus,  these  modes  are  determined  using  the  compartment  model, 
allowing  unrestrained  longitudinal  displacements,  when  a  simplified  analysis 
is  being  performed. 

In  summary,  the  ELSHOK  code  allows  a  user  to  perform  the  modal  analy¬ 
sis  of  the  main  body  S  by  employing  B0S0R4  in  conjunction  with  either  of 
two  different  modeling  techniques.  The  first  is  the  usual  straightforward 
manner  in  which  a  mathematical  model  is  developed  and  is  employed  for  all 
values  of  N  Included  in  a  given  analysis.  The  second  is  an  approximate 
technique,  allowing  a  user  to  highlight  a  portion  of  the  body  S.  This 
approach  uses  (1)  the  full  model  for  the  N=1  modes  and  the  N=0  rolling  and 
torsional  modes  and  (2)  the  compartment  model  for  the  N=0  axisymmetrical  or 
breathing  modes  and  the  modes  having  N>2  .  Once  the  mathematical  model  of  the 
shell  S  has  been  decided  upon,  the  user  determines  the  in-vacuo  free-free 
modes  and  natural  frequencies  routinely  using  the  B0S0R4  code. 
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3.  Rules  of  Operation 


For  compatibility  with  the  current  version  of  the  ELSHOK  code,  the  fol¬ 
lowing  rules  must  be  observed  when  executing  the  B0S0R4  code: 

(1)  Symmetry  about  the  longitudinal  midpoint  of  the  shell  may  not  be 
employed  in  the  mathematical  model  of  the  shell  unless  the  entire 
structure  under  consideration  (including  the  shell  S,  any  internal 
appendages  a,  and  any  concentrated  masses)  is  symmetrical  about  the 
midpoint . 

(2)  When  symmetry  is  employed,  the  left-hand  portion  of  the  shell 
(as  depicted  in  Fig.  1)  must  be  modeled. 

(3)  When  the  approximate  modeling  technique  is  employed,  the  compartment 
model  must  be  identical  to  the  corresponding  portion  of  the  shell  in 
the  full  model. 

(4)  The  wet  segments  of  the  mathematical  model  (i.e.,  those  segments 
representing  the  portion  of  the  submerged  structure  in  contact  with 
the  fluid)  must  be  numbered  first,  starting  with  unity  at  the  initial 
segment,  at  the  left-hand  end  of  the  model.  Any  internal  segments 
may  be  numbered  in  any  manner  consistent  with  the  restrictions  of 
Ref.  [7]. 

(5)  A  separate  execution  of  B0S0R4  must  be  made  for  each  circumferential 
harmonic  N  included  in  an  analysis.  If  only  a  full  model  representa¬ 
tion  of  S  is  employed,  the  N=0  breathing  modes  and  the  N=0  rolling 
and  torsional  modes  may  be  determined  in  the  same  execution  of  B0S0R4, 
if  desired. 

(6)  Only  non-negative  values  of  N  (0,  +1,  +2,  ...)  must  be  used.  The 
resulting  modes  produced  by  B0S0R4  are  used,  with  appropriate  sign 
changes,  in  the  subsequent  shock  response  calculation  as  modes 
having  circumferential  symmetry  with  respect  to  the  origin  of  the 
circumferential  coordinate  angle  9  of  Fig.  1.  When  antisymmetrical 
modes  are  required,  they  are  also  accounted  for  by  means  of  appropriate 
sign  changes.  The  sign  convention  used  in  ELSHOK,  to  which  the  com¬ 
ponents  of  the  modes  obtained  from  B0S0R4  are  made  to  conform,  is 
shown  in  Fig.  1. 
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(7)  In  any  given  shock  response  problem,  modes  must  be  computed  and 
saved  for  at  least  N=0  and  N*l. 

4.  Modifications  to  BOSOR4 

A  number  of  minor  changes  have  been  made  to  B0S0R4  to  facilitate  its 
use  in  the  ELSHOK  suite  of  computer  programs.  A  brief  description  of  each 
of  these  modifications  follows: 

(1)  The  capacity  of  B0S0R4  has  been  Increased  from  50  modes  to  100  modes 
for  each  circumferential  wave  number  N  considered. 

(2)  An  additional  file  has  been  supplied  and  assigned  to  logical  unit  3. 

This  file  is  referred  to  by  the  local  or  logical  file  name  TAPES  and 
corresponds  to  the  shell  mode  file  for  a  given  N  described  in  Section 
III.  File  TAPES  contains  shell  geometry,  physical  properties,  generalized 
masses,  natural  frequencies,  and  mode  shapes,  in  coded  format,  for  a 
given  circumferential  harmonic  N.  The  shell  mode  file  must  be  saved 
after  each  execution  of  B0S0R4  to  provide  data  needed  subsequently. 

(3)  The  plotting  routines  supplied  with  the  standard  B0S0R4,  and  all  sub¬ 
routine  calls  Involving  them,  have  been  eliminated.  Plots  of  the 
mode  shapes  obtained  from  B0S0R4  are  produced  from  the  data  on  TAPES 
by  means  of  a  small  plotting  processor  (BOSRMP)  supplied  with  ELSHOK. 

This  processor  prepares  plots  on  TEKTRONIX  graphics  terminals. 

(4)  One  additional  input  item  has  been  added  after  the  usual  input  data. 

The  purpose  of  this  additional  data  is  to  provide  a  clean  or  errror-free 
termination  of  a  B0S0R4  execution.  Appropriate  documentation  may  be 
found  in  the  next  subsection. 

The  above  changes  have  been  made  through  the  use  of  the  UPDATE  batch 
editing  system  (Ref.  [9])  and  may  be  easily  removed  to  return  B0S0R4  to 
its  standard  configuration. 


5.  Additional  Input 


When  using  B0S0R4  with  the  ELSHOK  code,  all  of  the  input  data  described 
in  Ref.  [7]  must  be  prepared  for  the  problem  of  interest.  Moreover,  one 
additional  input  item  should  be  placed  after  the  usual  set  of  data  to 
terminate  the  execution  of  B0S0R4  without  error.  A  description  of  this 
input  item  follows: 

ITEM  NO.  VARIABLE  FORMAT 

last  TITLE(l)  AlO 


where  TITLE(l)  =  first  word  of  the  alphameric  title  array  employed  by 
B0S0R4.  The  user  should  supply  the  value  ENDbDATAbb, 
in  which  the  symbol  "b"  denotes  a  blank  space. 


Note:  If  TITLE (l)=ENDbDATAbb,  an  end  of  file  mark  is  written  on 

logical  file  TAPE3,  and  the  execution  of  B0S0R4  is  terminated 
by  a  call  to  the  system  routine  EXIT.  If  this  value  is  not 
supplied,  the  execution  of  B0S0R4  proceeds  normally  but  ter¬ 
minates  with  a  system  error  indicating  an  attempt  to  read  more 
input  data  than  supplied.  The  use  of  this  additional  input 
item  is  optional. 


6.  Sample  Control  Cards 

A  typical  set  of  control  cards  or  job  control  language  (JCL)  statements 
for  an  execution  of  B0S0R4  under  a  NOS  1  operating  system  Is  shown  below.  A 
brief  outline  of  these  control  cards  Is  presented  In  the  next  subsection. 
Various  other  arrangements  of  the  job  steps  are  possible. 


B0S0R4(T277)  RANLET 
USER<RANLETOtPASWORO) 

CHARGE (60220 IRtPROJECT) 
SETBSL(1777) 

COPYBR ( iNPUTt INFILE) 
REWIND(INFILE) 

PURGE (BOSSYOO/NA) 

DEFINE (TAPE3>BOSSYOO/CT»StHsR) 
ATTACH (LGOsBOSFTN) 

RFL(125000) 

LOSET (PRESETsINOEFfHAPsN) 
LGO(INFILE) 

GOTO (CATALOG) 

EXIT. 

CATALOG(TAPE3tR) 

REWIN0(TAPE3*INFILE) 

C0PYBF(TAPE3.C0PY) 

COPYSBFdNFILEfCOPY) 

ROUTE (COPYtDC=PR*TID»WEIDLNG) 


BOSOR4(T277)  RANLET 

Job  card  or  job  statement  containing  an  alphameric  job  name,  a  time 
limit  In  seconds  (octal),  and  an  optional  identifier. 

USER (RANLETD , PASWORD) 

CHARGE (60220 IR , PROJECT) 

SETBSL(1777) 

System  dependent  accounting  information. 

COPYBRdNPUT,  INFILE) 

REWIND (INFILE) 

Provide  copy  of  input  data  card  Images  on  local  file  INFILE. 

ATTACH (LGO-BOSFTN) 

Assign  permanent  file  BOSFTN  to  current  job  with  local  file  name  LGO. 

File  BOSFTN  contains  the  compiled  code  or  object  module  for  B0S0R4. 

PURGE (BOSSYOO/NA) 

DEFINE (TAPE3=B0SSY00/CT=S ,M=R) 

Assign  space  on  permanent  file  device  for  the  shell  mode  file  (TAPE3) 
under  the  permanent  file  name  BOSSYOO.  Information  written  on  TAPE3 
during  the  execution  of  B0S0R4  is  written  directly  on  the  permanent  file. 
RFL( 125000) 

Request  field  length  for  execution.  This  is  optional  on  many  systems. 
LDSET (PRESET-INDEF , MAP=N) 

Select  load  options. 

LGO (INFILE) 

Load  and  execute  the  B0S0R4  object  module.  The  required  set  of  input 
data  is  supplied  on  local  file  INFILE. 

GOTO (CATALOG) 

Unconditional  transfer  to  control  statement  CATALOG  if  the  execution 
terminates  normally. 

EXIT. 

Provide  transfer  of  control  when  an  execution  error  occurs.  If  this 
statement  is  encountered  after  an  error-free  program  execution,  the 
system  terminates  the  job.  If  an  error  condition  occurs  during  job 
execution,  the  system  processes  the  control  stat  ments  that  follow 
the  EXIT  statement. 

CATALOG (TAPE3,R) 

List  information  about  each  logical  record  of  file  TAPE3  on  file  OUTPUT. 
REWIND (TAPE3 , INFILE) 

C0PYBF(TAPE3,C0PY) 

COPYSBF ( INFILE , COPY) 

Copy  the  information  on  files  TAPE3  and  INFILE  to  the  local  file  COPY. 
ROUTE (COPY , DC=PR , TID-WEIDLNG) 

Dispose  the  file  COPY  to  the  remote  site  WEIDLNG.  At  job  termination, 
the  system  changes  the  local  file  OUTPUT  into  a  print  file  which  is 
disposed  to  the  site  corresponding  to  the  job  origin. 
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8.  Input  Data  for  Sample  Problem  1 


A  set  of  input  data  applicable  to  Sample  Problem  1  is  shown  below. 

These  input  cards  or  card  images,  depending  on  the  mode  of  operation  selected 
by  the  user,  direct  the  B0S0R4  code  to  perform  the  modal  analysis  of  the 
shell  S  of  Fig.  4  for  the  circumferential  harmonic  N=0.  For  each  of  the 
other  circumferential  harmonics  employed  (N=l,2,3),  the  parameters  NMINB  and 
NMAXB  on  the  fifth  card  must  be  set  to  the  value  of  N  being  considered  and 
the  alphameric  job  title  must  be  modified  accordingly. 
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9.  Input  Data  for  Sample  Problem  2 

The  set  of  input  data  listed  below  directs  the  B0S0R4  code  to  perform 
the  modal  analysis  of  the  full  model  of  the  shell  S,  shown  schematically  in 
Fig.  14,  for  the  circumferential  harmonic  N=l. 
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For  the  circumferential  harmonic  N=0,  the  set  of  input  data  shown  below 
directs  the  BOSOR4  code  to  perform  the  modal  analysis  of  the  compartment  model 
of  the  shell  S  (Fig.  14).  For  N=2  and  N=3,  the  parameters  NMINB  and  NMAXB  on 
the  fifth  card  must  be  set  to  the  value  of  N  being  considered  and  the  alpha¬ 
meric  job  title  must  be  modified  accordingly. 
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APPENDIX  B  -  THE  USE  OF  ACESNID 


1.  Background 

The  ACESNID  computer  code  of  Ref.  [10]  is  included  in  the  ELSHOK  suite 
of  computer  programs  for  determining  the  virtual  mass  array,  a  quantity 
producing  the  late-time  contribution  of  the  DAA.  Spatial  functions  having 
the  property  of  orthogonality  over  the  wet  surface  of  a  submerged  body  are 
required  in  the  formulation  upon  which  ACESNID  is  based.  These  surface 
expansion  functions  are  obtained  during  the  calculation  of  the  virtual  or 


entrained  mass.  The  ACESNID  code  is  quite  general  and  considers  a  cavity 


in  an  infinite  acoustic  fluid  corresponding  to  the  wet  surface  of  a  fully 
submerged  shell  of  revolution  having  an  arbitrary  meridional  geometry.  A 
knowledge  of  the  information  given  in  Refs.  [2]  and  [10]  and  Sections  II 
and  III  of  the  present  report  will  aid  in  the  understanding  of  this 
appendix. 

In  ACESNID,  the  motion  of  the  fluid  induced  by  prescribed  steady- 
state  vibrations  is  governed  by  a  velocity  potential  which  satisfies 
the  Helmholtz  equation.  The  interior  boundary  of  the  fluid,  a  cavity 
of  revolution,  is  divided  into  bands.  A  line  distribution  of  simple 
sources  is  placed  on  the  central  line  of  each  band,  and  acoustic 
influence  coefficients  for  pressure  and  velocity  are  determined  for  a 
given  frequency  of  excitation.  The  source  strengths  are  then  evaluated 
by  equating  the  normal  fluid  velocity  on  the  surface  of  the  cavity  to  pre¬ 
scribed  values  corresponding  to  steady-state  motions  with  spatial 
amplitude  distributions  given  by  surface  expansion  functions.  From  the 
source  strengths  follow  the  fluid  pressures  developed  on  the  surface 
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of  the  cavity  and  the  generalized  forces  required  to  determine  the  accession 
to  inertia  and  damping  associated  with  the  chosen  surface  expansion  functions 
When  applying  ACESNID  with  ELSHOK,  a  steady-state  problem  of  sufficiently 
low  frequency  must  be  considered  to  approximate  the  incompressible  behavior 
of  the  fluid.  In  such  a  case,  the  solution  to  the  Helmholtz  equation 
approaches  the  solution  to  the  Laplace  equation.  The  resulting  accession 
to  damping  becomes  negligible,  and  the  accession  to  inertia  yields  the 
virtual  mass  array,  as  discussed  in  Ref.  [2].  The  formulation  employed 
in  ACESNID  is  based  upon  a  separation  of  variables,  and  allows  the  user 
to  treat  all  of  the  required  values  of  the  circumferential  harmonic  N  in 
one  execution  or  in  a  series  of  separate  calculations. 

2.  Partitioning  of  Cavity 

In  all  applications  of  ELSHOK,  the  use  of  ACESNID  must  parallel  the 
use  of  B0S0R4.  If  the  structural  modeling  proceeded  in  the  straight¬ 
forward  manner  of  using  only  the  full  model  representation  of  the  shell 
S,  consistent  partitioning  of  the  cavity  in  the  fluid  must  be  provided 
for  the  generation  of  the  surface  expansion  functions.  Thus,  the  user 
may  not  highlight  a  region  of  interest.  On  the  other  hand,  if  the  user 
chose  to  highlight  a  portion  of  the  structural  model  and  to  apply  the 
full  and  compartment  models  in  the  fashion  described  in  Section  II  and 
Appendix  A,  the  partitioning  of  the  cavity  must  reflect  the  highlighted 
region  of  interest  to  allow  the  evaluation  of  the  appropriate  surface 
expansion  functions. 

When  generating  surface  expansion  functions  for  use  with  the 
straightforward  full  model  representation  of  the  shell  (no  region  of 
interest),  the  cavity  in  the  fluid,  or  the  wet  surface  of  the  submerged 
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structure,  must  in  general  be  partitioned  into  three  distinct  regions: 

(1)  the  left  region  corresponding  to  the  left  end  closure;  (2)  the  right 
region  corresponding  to  the  right  end  closure;  and  (3)  the  central  region, 
which  need  not  be  a  right  circular  cylinder,  lying  between  the  left  and 
right  regions.  When  symmetry  about  the  longitudinal  midpoint  of  the 
structure  is  employed  in  the  full  model,  specification  of  the  right  region 
is  not  required.  In  the  straightforward  mode  of  operation,  all  cir¬ 
cumferential  harmonics  considered  in  a  problem  under  study  may  be  treated 
in  one  execution  of  ACESNID.  Equations  (13)-(15)  govern  the  generation 
of  the  surface  expansion  functions  for  this  case. 

When  generating  surface  expansion  functions  for  use  with  the  simplified 
representation  of  the  shell,  the  cavity  in  the  fluid  must  in  general  be 
partitioned  into  four  distinct  regions.  The  first  three  are  identical 
to  those  described  above.  The  fourth  region  corresponds  to  the  compartment 
model  or  region  of  interest  within  the  full  model  of  the  shell  S.  Because 
different  requirements  must  be  met  for  N=1  and  for  N^l,  as  discussed  in 
Section  II  and  Appendix  A,  a  separate  execution  of  ACESNID  is  required 
for  each  value  of  N  Included  in  a  given  analysis.  Equations  (13)-(16)  govern 
the  generation  of  the  surface  expansion  functions  for  this  case. 

3.  Rules  of  Operation 

For  compatibility  with  the  current  version  of  the  ELSHOK  code,  the 

following  rules  must  be  observed  when  executing  the  ACESNID  code: 

(1)  The  use  of  the  B0S0R4  code  must  precede  the  use  of  ACESNID.  This 
is  so  because  the  ACESNID  code  obtains  the  geometry  of  the  cavity 
in  the  fluid  from  a  shell  mode  file  for  a  full  model. 


(2)  The  values  of  the  circumferential  harmonic  N  employed  with 
ACESNID  must  be  identical  to  those  used  with  B0S0R4.  If  N=0, 

1,2,3  were  used  for  the  modal  structural  analysis  of  the  main 

body  S,  then  N=0,l,2,3  must  be  used  for  the  calculation  of  the 
virtual  mass. 

(3)  For  the  generation  of  the  surface  expansion  functions,  the  parti¬ 
tioning  of  the  cavity  in  the  fluid  must  be  described  with 
reference  to  the  full  model  of  the  shell  S. 

(4)  The  frequency  of  excitation  specified  by  the  user  must  be  suf¬ 
ficiently  small  that  the  behavior  of  the  fluid  in  the  ACESNID 
calculation  approximates  that  of  an  incompressible  fluid. 

(5)  When  a  compartment  region  or  region  of  interest  is  considered  , 
a  separate  execution  of  ACESNID  must  be  performed  for  each 
circumferential  wave  number  N  included  in  an  analysis.  When 
straightforward  modeling,  involving  only  a  full  model,  is 
employed,  all  required  values  of  N  may  be  treated  in  one  execution 
of  ACESNID. 

(6)  The  virtual  mass  file  must  contain  the  results  computed  for  each 
circumferential  harmonic  N  in  the  order  of  increasing  number  of 
circumferential  waves  (i.e.,  in  the  order  N=0,  N=l,  N=2,...). 

4.  Modifications  to  ACESNID 

Reference  [10]  documents  a  version  of  ACESNID  prepared  for  use 
with  the  EPSA  code  of  Ref.  [Ij.  The  version  of  ACESNID  supplied  in 
the  ELSHOK  suite  of  computer  programs  differs  somewhat  from  the  ver¬ 
sion  of  Ref.  [1].  A  brief  description  of  each  of  the  major  differences 
between  these  two  versions  of  the  ACESNID  code  follows: 

(1)  Any  finite  cavity  of  revolution  may  be  considered  with  ELSHOK. 

(2)  The  input  data  specifications  have  been  changed.  Increased 
flexibility  has  been  provided  to  allow  description  of  more  general 
cavities  of  revolution  for  both  the  generation  of  the  surface 


expansion  functions  and  the  calculation  of  the  acoustic  Influence 
coefficients.  The  source  regions*  which  pertain  to  the  calculation  of 
the  acoustic  Influence  coefficients,  are  described  or  partitioned 
separately  from  the  regions  associated  with  the  surface  expansion 
functions. 

(3)  Numerical  integrations  are  performed  using  the  1/3  rule  due  to  Simpson 
(Ref.  [12]).  The  trapezoidal  integration  rule  employed  in  Ref.  [10] 
has  been  replaced.  Some  reduction  in  execution  time  has  been  observed 
as  a  result  of  this  change. 

(4)  A  file  has  been  supplied  and  assigned  to  logical  unit  4.  This  file  Is 
referred  to  by  the  local  or  logical- file  name  TAPE4.  File  TAPE4  is  a 
shell  mode  file  (TAPE3  from  a  B0S0R4  calculation)  for  a  full  model  of 
the  shell  S.  It  is  used  in  ACESNID  to  describe  the  geometry  of  the 
cavity  In  the  fluid. 

(5)  File  TAPE?  has  been  added  and  assigned  to  logical  unit  7.  This  file 
receives  the  virtual  mass  array  in  coded  card-image  format  for  each 
circumferential  harmonic  N  considered  in  an  execution  of  ACESNID. 

File  TAPE?  must  be  saved  to  provide  data  needed  subsequently. 


.  Description  of  Input  Data 


A  description  of  the  input  data  required  for  an  execution  of  the 
ACESNID  code  is  given  below.  The  input  data  specifications  are  written 
in  a  style  similar  to  FORTRAN.  It  is  assumed  that  the  user  is  familiar 
with  this  programming  language. 

Input  Set  1  FORMAT  (8A10) 

TITLE 

where  TITLE  =  alphameric  job  title  or  job  heading. 

Input  Set  2  FORMAT  (515) 

NSTART,  NFINIS,  NFREQ,  NVMASS,  NCHECK 

where  NSTART  =  initial  circumferential  harmonic  N  to  be  used  in  the 
calculation. 

NFINIS  =  final  circumferential  harmonic  N  to  be  used  in  the 

calculation.  Calculations  are  done  for  an  increment 
of  unity  in  the  circumferential  harmonic. 

NFREQ  =  number  of  excitation  frequencies.  For  ELSHOK 
applications,  use  NFREQ  =  1. 

NVMASS  =  integer  flag  indicating  whether  or  not  a  calculation  of 
virtual  mass  is  to  be  performed  (l=yes;  0=no) .  For 
ELSHOK  applications,  use  NVMASS  =  1. 

NCHECK  =  integer  flag  indicating  whether  or  not  the  computed 
source  strengths  are  to  be  checked  by  back  substitu¬ 
tion  into  the  linear  algebraic  equations  governing 
their  solution  (l=yes;  0=no) .  For  most  ELSHOK 
applications,  use  NCHECK  =  0. 

NOTE;  Since  ACESNID  utilizes  main  memory  management,  it  is  quite  easy 
to  reduce  the  core  requirement  or  central  memory  to  the  minimum 
size  required  for  a  given  execution.  This  may  be  accomplished 
by  specifying  a  negative  value  of  NSTART  in  an  otherwise  com¬ 
plete  set  of  input  data.  An  execution  of  ACESNID  with  this  set 
of  input  data  will  then  supply  information  from  which  the  desired 
size  of  the  dynamic  array  may  be  determined.  When  determining 
the  minimum  allowable  memory  size  for  a  given  problem,  the 
dimension  of  the  dynamic  array  in  blank  common  must  be  set  to  at 
least  NSEG+NFREQ+1  words,  where  NSEC  is  the  number  of  segments 
in  the  shell  mode  file  supplied  as  TAPE4. 


INPUT  SET  3 


FORMAT  (515) 


NLEFT,  NCYL,  NRITE,  NWBOSG,  NSYMF 

where  NLEFT  =  number  of  bands  on  the  left  end  closure  of  the  cavity 
In  the  fluid. 

NCYL  =  number  of  bands  on  the  cylindrical  portion  of  the 

cavity  in  the  fluid.  If  no  cylindrical  portion  exists 
for  the  problem  under  study,  use  NCYL  =  0. 

NRITE  =  number  of  bands  on  the  right  end  closure  of  the  cavity 
in  the  fluid.  This  value  must  be  supplied,  since 
ACESNID  requires  the  entire  cavity  in  the  fluid,  even 
if  a  model  of  half  of  the  submerged  shell  S  is  pro¬ 
vided  on  TAPE4  because  of  symmetry. 

NWBOSG  =  number  of  wet  segments  in  the  B0S0R4  model  contained 
on  the  shell  mode  file  supplied  as  TAPE4. 

NSYMF  =  Integer  flag  indicating  whether  or  not  the  cavity  in 
the  fluid  must  be  completed  using  symmetry  (l=yes; 

0=no) .  If  the  geometry  of  the  full  model  on  TAPE4 
describes  half  of  an  entire  structural  model,  then 
use  NSYMF  =  1. 

Note  1  ;  The  input  parameters  NLEFT,  NCYL,  NRITE  are  used  to  partition 
the  entire  Interior  boundary  of  the  fluid  into,  in  general  , 
three  distinct  regions  which  may  be  termed  "source  regions". 
The  cylindrical  portion  of  the  cavity  in  the  fluid  may  not 
be  present  In  a  given  problem.  In  essence,  these  input  param¬ 
eters  specify  the  number  of  lines  of  simple  sources 
employed  to  describe  the  behavior  of  the  fluid.  The  use  of 
source  regions  allows  the  acoustic  influence  coefficients  to 
be  computed  in  an  efficient  manner. 

Note  2  ;  Although  the  source  band  parameters  do  not  apply  to  the 

generation  of  surface  expansion  functions,  their  values  are 
governed  by  the  number  of  surface  expansion  functions 
selected.  This  is  so  because  compatibility  of  normal  motion 
is  required  at  the  structure-fluid  interface.  Thus,  a  suf¬ 
ficient  number  of  source  bands  must  be  provided  to  enable 
the  fluid  response  to  match  the  normal  motion  specified  by 
the  surface  expansion  functions.  For  most  applications  of 
ACESNID,  about  four  or  five  bands  are  needed  in  the  shortest 
half-wavelength  or  shortest  local  undulation  found  in  the 
set  of  surface  expansion  functions. 

Note  3;  In  many  applications  of  ACESNID,  the  end  closures  of  the 
fluid  cavity  correspond  to  the  end  closures  of  the 
structural  model,  and  the  central  portion  of  the  shell  of 
revolution  is  a  circular  cylinder.  However,  this  is  not  a 
requirement.  The  user  is  advised  to  refer  to  Ref.  [10]  for 
more  detailed  information. 


FORMAT  (515) 


Input  Set  4 


NORDER,  NFENDS, 
where  NORDER 


NFENDS 

NFCENT 

NFCMPT 

NOMIT 


NFCENT,  NFCMPT,  NOMIT 

=  control  integer  for  generation  of  surface  expansion 
functions . 

NORDER  =  1  :  surface  expansion  functions  are  alter¬ 
nately  symmetrical  and  antisymmetrical  about  the  mid¬ 
point  of  the  structure.  Use  this  option  for  a  general 
unsymmetrical  problem. 

NORDER  =  2  :  surface  expansion  functions  are  symmetrical 
about  the  midpoint  of  the  structure.  Use  this  option 
for  a  symmetrical  problem. 

NORDER  =  3  :  surface  expansion  functions  are  antisym¬ 
metrical  about  the  midpoint  of  the  structure.  Use 
this  option  for  an  antisymmetrical  problem. 

=  number  of  surface  expansion  functions  to  be  generated 
on  the  end  closures  of  the  structure.  An  equal  number 
of  surface  expansions  is  generated  on  each  end  closure 
unless  indicated  otherwise  by  Input  parameter  NOMIT. 

=  number  of  surface  expansion  functions  to  be  generated 
on  the  central  portion  of  the  structure. 

=  number  of  surface  expansion  functions  to  be  generated 
on  the  compartment  portion  of  the  structure.  If  a 
compartment  model  is  not  employed,  used  NFCMPT  =  0. 

=  control  parameter  for  omitting  surface  expansion 
functions  on  end  closures  of  structure. 

NOMIT  =  0  :  no  surface  expansion  functions  omitted. 

NOMIT  =  1  :  no  surface  expansion  functions  generated  on 
left  end  closure  of  structure. 

NOMIT  =  2  :  no  surface  expansion  functions  generated  on 
right  end  closure  of  structure. 


Input  Set  5 


FORMAT  (3E10.0) 


RHOFL,  VSOUND,  ERR 


where  RHOFL  =  mass  density  of  fluid. 

VSOUND  =  speed  of  sound  in  fluid. 

ERR  =  relative  accuracy  accepted  for  numerical  integration. 
This  parameter  terminates  the  numerical  integration 
procedure  used  in  ACESNID.  In  general,  use 
0.005  <  ERR<  0.0075. 


Input  Set  6 


FORMAT  (1015) 


LGEONP,  LFCNNP.  LSEFNP,  LORTH,  LGEOSP, 
LSEFSP,  LINF,  LRHS,  LSOL,  LACID 


where  LGEONP  =  output  flag  for  structural  nodal  point  geometry. 

LFCNNP  =  output  flag  for  convenient  functions  used  to 
generate  surface  expansion  functions. 

LSEFNP  =  output  flag  for  surface  expansion  functions 
evaluated  at  structural  nodal  points. 

LORTH  =  output  flag  for  orthogonality  check  of  surface 
expansion  functions. 

LGEOSP  =  output  flag  for  geometry  of  source  bands. 

LSEFSP  =  output  flag  for  surface  expansion  functions  evaluated 
at  central  lines  of  source  bands. 

LINF  =  output  flag  for  acoustic  influence  coefficients. 

LRHS  =  output  flag  for  right-hand  sides  of  equations 
for  determining  source  strengths. 

LS0L'=  output  flag  for  source  strengths. 

LACID  =  output  flag  for  accession  to  inertia  and  damping. 

Note  1  ;  In  an  output  flag  is  set  to  1,  output  is  supplied.  If  an 
output  flag  is  set  to  0,  no  output  is  supplied. 

Note  2  ;  In  almost  all  cases,  use  LINF  =  0,  since  the  acoustic 

Influence  coefficients  are  very  large  complex  arrays  for 
most  problems  of  interest. 

Input  Set  7  FORMAT  (8E10.0) 


(CPS(J),  J  =  1,  NFREQ) 

where  CPS(J)  =  the  J-th  frequency  of  excitation.  For  ELSHOK 

applications,  supply  only  CPS(l),  since  NFREQ  =  1, 
as  stipulated  above. 


Note :  The  frequency  of  excitation  supplied  as  CPS(l)  must  be  suf¬ 
ficiently  low  to  produce  steady-state  incompressible  behavior 
of  the  fluid.  For  ACESNID,  it  has  been  found  that  an 
excitation  frequency  is  sufficiently  low  when  the  quantity 
Wga^/c  is  of  the  order  of  0.0001  or  smaller.  Here,  is 
the  circular  frequency  of  excitation,  a^  is  the  maximum  radius 
of  the  central  portion  of  the  shell  S,  and  c  is  the  speed  of 
sound  in  the  fluid. 
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FORMAT  (415) 


Input  Set  8 

DO  10  K  =  1,  4 

10  READ  JSGBEG,  JPTBEG,  JSGEND,  JPTEND 

where  JSGBEG  =  segment  number  defining  the  exordium  of  the  K-th  region 
for  generating  surface  expansion  functions. 

JPTBEG  =  nodal  point,  of  segment  JSGBEG,  defining  the  exordium 
of  the  K-th  region  for  generating  surface  expansion 
functions. 

JSGEND  =  segment  number  defining  the  terminus  of  the  K-th 

region  for  generating  surface  expansion  functions. 
JPTEND  =  nodal  point,  of  segment  JSGEND,  defining  the  terminus 
of  the  K-th  region  for  generating  surface  expansion 
functions. 

K  =  1  :  supply  segment  and  nodal  point  numbers  for  wet 
surface  of  left  end  closure  of  structure  S. 

K  =  2  :  supply  segment  and  nodal  point  numbers  for  wet 
surface  of  central  region  of  structure  S. 

K  =  3  :  supply  segment  and  nodal  point  numbers  for 
wet  surface  of  right  end  closure  of  structure  S. 

K  =  4  :  supply  segment  and  nodal  point  numbers  for 
wet  surface  of  compartment  region  of  structure  S. 

Note  1  ;  The  segment  and  nodal  point  numbers  mentioned  above  refer 
to  the  numeration  supplied  by  the  B0S0R4  code  on  a  shell 
mode  file  for  a  full  model  representation  of  the  main  body 
S.  The  nodal  points  are  actually  points  on  the  reference 
meridian  used  to  generate  the  shell  of  revolution  S.  They 
are  numbered,  starting  from  unity,  within  each  segment. 

If  the  number  of  surface  expansion  functions  associated 
with  a  given  K-th  region  above  is  zero ,  omit  the  input 
parameters  for  that  value  of  K.  If  NSYMF  =  1,  omit  the 
input  parameters  for  the  right  end  closure  (K  *  3). 


Note  2  s 


Input  Set  9 

DO  20  K  =  1,  3 

20  READ  JSGBEG,  JPTBEG,  JSGEND,  JPTEND 


FORMAT  (415) 


where  JSGBEG  = 


JPTBEG  = 


JSGEND  = 


JPTEND  = 


segment  number  defining  the  exordium  of  the  K-th  source 
region  for  evaluating  acoustic  influence  coefficients, 
nodal  point,  of  segment  JSGBEG,  defining  the  exordium 
of  the  K-th  source  region  for  evaluating  acoustic 
influence  coefficients. 

segment  number  defining  the  terminus  of  the  K-th  source 
region  for  evaulating  acoustic  influence  coefficients, 
nodal  point,  of  segment  JSGEND,  defining  the  terminus 
of  the  K-th  source  region  for  evaluating  acoustic 
influence  coefficients. 

K  =  1  :  supply  segment  and  nodal  point  numbers  for 
left  end  closure  of  cavity  in  fluid. 

K  =  2  :  supply  segment  and  nodal  point  numbers  for 
cylindrical  portion  of  cavity  in  fluid. 

K  =  3  :  supply  segment  and  nodal  point  numbers  for 
right  end  closure  of  cavity  in  fluid. 


Note  1  ;  The  segment  and  nodal  point  numbers  mentioned  above  refer 

to  the  numeration  supplied  by  the  B0S0R4  code  on  a  shell  mode 
file  for  a  full  model  representation  of  the  main  body  S.  The 
nodal  points  are  actually  points  on  the  reference  meridan  used 
to  generate  the  shell  of  revolution  S.  They  are  numbered, 
starting  from  unity,  within  each  segment. 

Note  2  ;  If  NSYMF  =  1,  omit  the  input  parameters  for  the  right  end 
closure  (K  =  3) . 
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6.  Sample  Control  Cards 


A  sample  set  of  control  cards  or  job  control  language  statements 
for  an  execution  of  ACESNID  under  a  NOS  1  operating  system  is  shown  below. 
The  next  subsection  contains  a  brief  description  of  these  control  state¬ 
ments.  Various  other  arrangements  of  the  job  steps  are  possible. 


ACESNI0(T2777)  RANLET 
USER (RANLETOtPASWORD) 

CHARGE (602201RtPROJECT) 
SETBSL<7777) 

ATTACH (OLOPL=ACESNIO) 
UPDATE(F*LaA124) 

RETURN (OLDPL) 

FTN(IaCOMPILE*L*0#OPTa2*A) 

RETURN(COMPILE) 

ATTACH (TAPE4=BOSSYOO) 
RFL(165000) 

LDSET (PRESET»IN0EF,MAP=SB) 
LGO. 

GOTO (PURGE) 

EXIT. 

PURGE (VIRMAS/NA) 

SAVE(TAPE7aVIRMAS/CTaS.M*R) 

REWIND(TAPE7) 

C0PYSBF(TAPE7f OUTPUT) 


If  control  statements  of  the  above  type  are  used  to  perform  a  virtual 
mass  calculation  for  a  given  circumferential  harmonic  N,  the  files  saved 
from  the  local  file  TAPE?  for  the  individual  values  of  N  may  be  merged 
to  produce  the  virtual  mass  file  required  by  ELSHOK  in  the  following  manner: 


VMFILE(T7)  RANLET 
USER (RANLETOf PASWORO) 

CHARGE (602201RtPROJECT) 

GET(VIRMAS0.VIRMASlfVIRMAS2fVIRMAS3) 

COPYBR(VIRMASOfVMFILE) 

COPYBR (VIRMASl fVMFILE) 
C0PYBR<VIRMAS2.VMFILE) 

COPYBR (VIRMAS3tVMFILE) 

PURGE (VIRMAS/NA) 
SAVE(VMFILE*«VIRMAS/CTaStM=R) 

EXIT. 
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7.  Description  of  Control  Cards 

ACESNID(T2777)  RANLET 

Job  card  or  job  statement. 

USER (RANLETD , PASWORD) 

CHARGE(602201R, PROJECT) 

SETBSL(7777) 

System  dependent  accounting  information. 

ATTACH (OLDPL=ACESNID) 

Assign  permanent  file  ACESNID  to  current  job  with  local  file  name  OLDPL. 
File  ACESNID  contains  the  UPDATE  program  library  for  the  ACESNID  code. 
UPDATE(F,L=A124) 

RETURN (OLDPL) 

Process  UPDATE  directives  supplied  iu  the  job  stream  (on  file  INPUT) 
and  prepare  a  source  file  suitable  for  compilation. 

FTN ( I=COMP ILE , L=0 , 0PT=2 , A) 

RETURN (COMPILE) 

Produce  compiled  code  or  object  module  for  the  ACESNID  code. 

ATTACH (TAPE4=BOSSYOO ) 

Assign  permanent  file  BOSSYOO  to  current  job  with  local  file  name  TAPE4. 
File  BOSSYOO  contains  a  shell  mode  file  for  the  full  model  representation 
of  the  shell  S. 

RFL( 165000) 

Request  field  length  for  execution.  This  is  optional  on  many  systems. 
LDSET (PRESET=INDEF , MAP=SB) 

LGO. 

Select  load  options,  load,  and  execute  the  compiled  version  of  ACESNID. 

The  required  set  of  input  data  is  supplied  in  the  job  stream  (on 
file  INPUT). 

GOTO (PURGE) 

Unconditional  transfer  to  control  statement  PURGE  if  the  execution 
terminates  normally. 

EXIT. 

Terminate  execution  or  provide  transfer  of  control  when  an  execution 
error  occurs. 

PURGE (VIRMAS/NA) 

SAVE (TAPE7=VIRMAS/CT=S ,M=R) 

Save  local  file  TAPE7  on  permanent  file  device  under  the  permanent 
file  name  VIRMAS.  File  TAPE7  contains  the  virtual  mass  information 
from  the  current  execution  of  ACESNID. 

REWIND (TAPE7) 

COPYSBF (TAPE7 , OUTPUT) 

Provide  copy  of  local  file  TAPE7  on  the  OUTPUT  file. 

Note;  The  job  control  statements  for  merging  the  results  of  several  ACESNID 

calculations  into  a  single  file  (virtual  mass  file)  are  self-explanatory. 
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UPDATE  Directives  and  Input  Data  for  Sample  Problem  1 


The  UPDATE  directives  employed  to  calculate  the  virtual  mass  for 


Sample  Problem  1  are  listed  below. 


*IDENT, STORAGE 

*DELETE , ACESNID . 26 .ACESNID. 27 
COMMON  A(40000) 

NWORDA  =  40000 

The  size  of  the  dynamic  array  A  shown  above  is  slightly  in  excess  of  the  size 
determined  by  a  restricted  execution  of  the  ACESNID  code.  In  this  restricted 
execution,  the  dimension  of  the  array  A  and  the  variable  NWORDA  were  set  to  50 
words  and  the  input  parameter  NSTART  was  set  to  -1.  All  other  input  parameters 
were  identical  to  those  listed  below.  When  the  dynamic  array  is  dimensioned 
to  40000  words,  as  above,  a  field  length  of  almost  165000  words  (octal)  is 
required  for  an  execution  of  the  ACESNID  code. 

The  set  of  input  data  shown  below  directs  the  ACESNID  code  to  perform 
the  virtual  mass  calculation  for  Sample  Problem  1  for  the  circumferential 
harmonics  N=0,l,2,3.  All  references  to  shell  geometry  correspond  to  the 
structural  data  presented  in  Appendix  A  for  Sample  Problem  1. 


VIRTUAL  HASS  FOR  RIN6-STIFFENE0  CYLINDER  MITH  FLAT  END  CAPS  —  N  >  0tl*2*3 


0 

3 

1 

1 

0  nstart.nfinis.nfreo.nvnass.ncheck 

12 

T5 

12 

2 

1  NLCFT.NCYLtNRITEtNMBOSG.NSYHF 

2 

3 

15 

0 

0  NOROER.NFENDStNFCENT.NFCHPTtNOHIT 

9.45216-5 

60. 

0*3 

0.005 

RHOFLtVSOUNOtERR 

1 

1 

1 

1 

0  0 

0011  OUTPUT  FLAOS 

0.10 

CPSd) 

1 

1 

1 

12 

SEFS 

JSOBEO  *  JPTBES  *  JS6EN0  *  JPTENO 

LEFT 

2 

1 

2 

53 

SEFS 

OSBBEB.JPTBES* JSBENO.JPTENO 

CENTRAL 

1 

1 

1 

12 

SOURCES 

JSBBEOf  JPTBES • JS6EN0 • jPTENO 

LEFT 

2 

1 

2 

53 

SOURCES 

JSSBEOt JPTBES tJSSEND* JPTENO 

CYLINDER 

The  UPDATE  directives  employed  to  calculate  the  virtual  mass  for 


Sample  Problem  2  are  listed  below. 


*IDENT, STORAGE 

*DELETE , ACESNID . 26 , ACESNID . 2  7 
COMMON  A(AOOOO) 

NWORDA  =  40000 


The  size  of  the  dynamic  array  A  shown  above  is  slightly  in  excess  of  the  size 
determined  by  a  restricted  execution  of  the  ACESNID  code.  In  this  restricted 
execution,  the  dimension  of  the  array  A  and  the  variable  NWORDA  were  set  to  50 
words  and  the  input  parameter  NSTART  was  set  to  -1.  All  other  input  parameters 
were  identical  to  those  listed  below.  When  the  dynamic  array  is  dimensioned 
to  40000  words,  as  above,  a  field  length  of  almost  165000  words  (octal)  is 
required  for  an  execution  of  the  ACESNID  code. 

The  set  of  input  data  shown  below  directs  the  ACESNID  code  to  perform 
the  virtual  mass  calculation  for  Sample  Problem  2  for  the  circumferential 
harmonic  N=l.  The  virtual  mass  so  produced  is  intended  for  use  with  the  full 
model  representation  of  the  shell  S.  All  references  to  shell  geometry 
correspond  to  the  full  structural  model  presented  in  Appendix  A  for  Sample 
Problem  2. 


VlATUAL  MASS  FOR  CVLlNOCR  «ITH  FLAT  ENOS  <FUtL  MOOED  —  N  s  1 


1 

10  T< 
2 

9.AS21A-< 

1 

O.li 

1 

2 

3 

1 

2 


1 


lO.ot: 


3 


0  NSTARTtNFINlS»NFREe«NVMAS5*NCHECK 
1  NLEFT,NCVL«NRITE*Ntf»OSO«NSYNF 
0  NOROERfNFENOS»NFCENT*NFCMFT*NOMIT 
O.eOTS  RHOFL*VSOIINO»ERR 

0  0  0  0  1  I  OUTPUT  FLASS 

CRSin 

SEFS  JS0BEe»JRTBEB*dSeEN0tJRTEN0  LEFT 

SEFS  JSBBEB«JRTBE8f3S6END*JRTEN0  CENTRAL 

SEFS  JS0BE6*JRTBEB»JSBEN0tJRTEN0  CMRT 

SOURCES  JSBBE6»JRTBE8*US0EN0«URTEN0  LEFT 

SOURCES  JSBBEBvJRTBEBvdSOENOtJRTENO  CVLINOER 


For  the  circumferential  harmonic  N=0,  the  virtual  mass  associated  with 
the  compartment  model  of  the  shell  S  may  be  determined  by  executing  the 
ACESNID  code  with  the  set  of  input  data  listed  below.  For  N=2  and  N=3, 
which  must  be  treated  separately  for  Sample  Problem  2,  the  parameters 
NSTART  and  NFINIS  must  be  set  to  the  value  of  N  being  considered  and  the 
alphameric  job  title  must  be  modified  accordingly. 


VIRTUAL  MASS  FOR  CYLINDER  MITH  FLAT  ENOS  ICMPT  MOOED  ••  N  ■  0 
00110  NSTARTvNFlNIStNFREOfNVHASStNCHECK 

10  TS  10  3  1  NLEFT*NCVL»NRITE«NMB0S6»NSYNF 

2  0  0  T  0  noroer*nfends*nfcent«nfcmpt»nomit 

9.4521A-5  60.0«3  0.00T5  RHOFLtVSOUNOtERR 

1111000011  OUTPUT  FLAOS 
0.10  CPSil) 

3  1  3  19  SEFS  JSOBES.JPTBEStJSOENOtJPTENO  CMPT 

1  1  1  12  SOURCES  JSBBEB.JPTBEB.JSBENO.JPTENO  LEFT 

2  1  3  19  SOURCES  JSBBEB.JPTBEB.^ISBENDt  JPTENO  CYLINDER 


APPENDIX  C  -  THE  USE  OF  PIFLASH 


1.  General  Information 

The  PIFLASH  computer  code  is  a  post-processor  for  the  B0S0R4  and 

ACESNID  codes.  As  discussed  in  Section  III,  the  PIFLASH  code  takes  data 

from  all  of  the  shell  mode  files  (produced  by  B0S0R4)  and  from  the  virtual 

mass  f^'le  (produced  by  ACESNID)  and  reorganizes  it  to  facilitate  the 

solution  of  the  equations  governing  the  transient  response  problem  under 
*) 

investigation.  Any  concentrated  or  point  masses  which  are  attached 
to  the  shell  S  are  accounted  for  during  the  execution  of  PIFLASH. 
Familiarity  with  Appendices  A  and  B  will  aid  in  the  understanding  of  this 
appendix. 

The  labels  or  headings  employed  to  identify  the  various  portions  of 
the  output  supplied  by  PIFLASH  correspond  to  the  terminology  used  with 
the  approximate  modeling  technique  available  for  the  shell  S,  even  if 
the  straightforward  modeling  technique  (with  no  region  of  interest)  is 
being  used  in  a  particular  problem.  Thus,  geometry,  mode  shapes,  sur¬ 
face  expansion  functions,  etc  are  labeled  as  applying  to  the  full  model 
when  they  are  associated  with  files  produced  for  N=0  (torsional)  and 
N=l.  Items  associated  with  N=0  (breathing)  and  N  >^2  are  labeled  as 
applying  to  the  compartment  model. 


*) 


The  current  version  of  ELSHOK  constructs  the  equations  governing  the  shock 
response  problem  for  the  entire  structure  and  for  the  entire  cavity  in  the 
fluid,  even  if  symmetry  is  employed  in  a  given  analysis.  Future  versions 
of  ELSHOK  may  be  based  upon  a  different  approach. 
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2.  Rules  of  Operation 

(1)  Data  from  the  B0S0R4  and  ACESNID  codes  must  be  available  prior  to  the 
execution  of  the  PIFLASH  code. 

(2)  Data  must  be  supplied  to  PIFLASH  for  at  least  N=0  (breathing)  and  N=l. 

(3)  The  virtual  mass  file  must  be  supplied  to  PIFLASH  as  a  local  file 
assigned  to  logical  unit  7  and  referred  to  as  TAPE?. 

(4)  The  shell  mode  files  must  be  supplied  to  the  PIFLASH  code  as  local  files 
assigned  to  logical  units  selected  from  the  set  10,11,12,13,14,15,20. 
These  files  must  be  referred  to  as  TAPEIO,  TAPEll,  etc.  Additional 
shell  mode  files  may  be  accomodated  by  changing  the  PROGRAM  card  of 

the  PIFLASH  code  through  the  use  of  UPDATE. 

(5)  When  the  problem  under  consideration  includes  Internal  appendages,  the 
user  should  supply  control  cards  or  control  statements  to  ensure  that 
the  shell-fluid  file  (TAPE99  on  logical  unit  99)  is  comprised  of 
system-logical  records.  This  allows  simple  control  card  operations 

to  be  used  to  append  data  for  internal  appendages  to  the  shell-fluid 
file.  If  system-logical  records  are  not  used,  an  additional  processor 
must  be  furnished  to  read  and  transfer  shell,  fluid,  and  substructure 
data  to  a  single  composite  file. 


(6)  File  TAFE99  must  be  saved  to  provide  data  needed  subsequently. 


Description  of  Input  Data 


The  set  of  Input  data  required  for  an  execution  of  the  PIFLASH  code  Is 
described  below.  Following  Appendix  B,  the  Input  data  specifications  are 


written  In  a  style  similar  to  FORTRAN, 


Input  Set  1 


NUMBER,  NTORSN,  NPTM,  NSYMS,  NSYMP 


FORMAT  (515) 


where  NUMBER  = 


NTORSN  = 


NSYMS  = 


NSYMP 


total  number  of  circumferential  harmonics  to  be  processed, 
exclusive  of  the  circumferential  harmonic  associated  with 
the  torsional  ana  rolling  modes  of  the  shell  S. 

Integer  flag  Indicating  whether  or  not  the  torsional  and 
rolling  modes  of  the  shell  S  are  to  be  processed 
(l=yes;  0=no) . 

number  of  concentrated  or  point  masses  attached  ,to  the 
shell  S. 

Integer  flag  Indicating  whether  or  not  symmetry  Is 
employed  about  the  longitudinal  midpoint  of  the  shell  S 
(l=yes;  0=no) . 

Integer  flag  Indicating  whether  or  not  symmetry  Is  to  be 
taken  Into  account  for  the  point  masses  (l=yes;  0=no) . 

If  NSYMP  =  1,  all  weights  of  point  masses  supplied  to 
PIFLASH  will  be  doubled  prior  to  the  formation  of  the 
equations  of  motion.  For  the  current  version  of  PIFLASH, 
use  NSYMP  =  NSYMS. 


FORMAT  (415) 


Input  Set  2 


(NWETSG(K),  NUSESG(K),  K  =  1,  2) 


where  NWETSG(K)  = 
NUSESG(K)  = 


number  of  wet  segments  in  the  type  K  model  of  the  shell 
of  revolution  S- 

number  of  segments  used  from  the  total  number  available 
in  the  type  K  model  of  the  shell  of  revolution  S 
[NUSESG(K)  >  NWETSG(K)]. 

K  =  1  :  input  corresponds  to  the  full  model  of  the  shell  S 
K  =  2  :  input  corresponds  to  the  compartment  model  of  the 
shell  S  if  the  approximate  modeling  technique  was  used. 

If  the  straightforward  modeling  technique  was  used,  the 
input  corresponds  to  the  full  model  of  the  shell  S  (i.e., 
same  values  supplied  for  K=1  and  K=2) . 


Note ;  The  components  of  the  modes  of  the  shell  S  are  written  on  file  TAPE99 
only  for  the  first  NUSESG(K)  segments  for  each  K  [i.e.,  for  segments 
1,2,...,  NUSESG(K)].  If  a  substructure  is  attached  to  a  given  part 
of  the  shell  S,  the  appropriate  segments  must  be  included  for  both 
K=1  and  K=2.  If  a  point  mass  is  attached  to  a  given  part  of  the 
shell  S,  the  appropriate  segment  must  be  included  for  K=1  and,  if  it 
exists,  for  K=2.  A  segment  may  not  exist  for  K=2  because  it  falls 
outside  of  the  region  of  interest  when  the  approximate  modeling 
technique  is  used.  If  a  segment  is  included  for  K=l,  the  corresponding 
segment  must  be  included  for  K=2,  if  it  exits,  to  allow  subsequent 
production  of  velocity-time  histories. 


FORMAT  (ElO.O) 


Input  Set  4  FORMAT  (1515) 

(NTAPE(J),  J  =  1,  NITEMS) 

where  NITEMS  =  total  number  of  circumferential  harmonics  to  be  processed, 
inclusive  of  the  circumferential  harmonic  associated  with 
the  torsional  and  rolling  modes  if  such  modes  are  active 
in  the  problem  under  study. 

NTORSN  =  0  :  NITEMS  =  NUMBER. 

NTORSN  =  1  :  NITEMS  =  NUMBER  +  1. 

NTAPE  =  array  containing  the  logical  unit  numbers  on  which  the 
shell  mode  files  are  supplied  to  PIFLASH.  If  torsional 
and  rolling  modes  are  employed,  the  logical  unit  number 
furnished  for  these  modes  must  be  the  last  element  of  the 
NTAPE  array. 

Note ;  The  logical  unit  numbers  in  the  NTAPE  array  must  be  supplied  in 

the  order  corresponding  to  N=0  (breathing),  N=l,  N=2 . 

N=NUMBER-1,  N=0  (torsional  and  rolling).  If  the  torsional  and 
rolling  modes  are  not  active,  only  the  first  NUMBER  elements  of 
the  NTAPE  array  must  be  supplied. 


Input  Set  5 


FORMAT  (1515) 


(NJUSE(J),  J  -  1,  NITEMS) 

where  NJUSE  =  array  in  which  each  element  indicates  the  number  of  modes 
of  the  shell  S  to  be  used  for  a  given  circumferential 
harmonic  N.  In  the  NJUSE  array,  the  number  of  modes  used 
for  each  value  of  N  must  correspond  in  order  to  the 
logical  unit  number  supplied  in  the  NTAPE  array  for  the 
same  value  of  N. 
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Input  Set  6 


FORMAT  (1515) 


(KORSG(K),  K  =  1,  NKORSG) 

where  NKORSG  =  size  of  segment  correspondence  array  KORSG.  If  the 

approximate  modeling  technique  was  used  for  the  shell  S,  then 
NKORSG  is  the  number  of  segments  in  the  compartment 
model.  If  the  straightforward  modeling  technique  was  used 
for  the  shell  S,  then  NKORSG  is  the  number  of  segments 
in  the  full  mathematical  model. 

KORSG(K)  =  segment  in  full  model  of  the  shell  S  which  corresponds 

to  segment  K  in  the  compartment  model.  If  straightforward 
modeling  was  used  for  the  shell  S,  supply  the  values 
1,2, . . . ,NSEG,  where  NSEG  is  the  total  number  of  segments 
in  the  full  mathematical  model. 

Input  Set  7  FORMAT  (1515) 

DO  30  K  =  1,  NITEMS 

30  READ  (JUSE(J,K),  J  =  1,  NJUSE(K)) 

where  JUSE  =  array  containing  the  identification  numbers  of  the  modes 
of  the  shell  S  selected  for  use  by  the  user.  The  number 
of  modes  used  and  the  circumferential  harmonic  N  associated 
with  a  given  column  K  of  the  JUSE  array  must  correspond  to 
element  K  of  the  NJUSE  array.  Within  each  column  K  of  the 
JUSE  array,  the  identification  numbers  must  appear  in 
increasing  order. 
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Input  Set  8  FORMAT  (E10.0,2I5,E10.0) 

(WEIGHT(J),  LSEGM(J),  LMESH(J) .  ANGPTM(J) ,  J  =  1,  NPTM) 

where  WEIGHT(J)  =  weight  of  the  J-th  point  or  concentrated  mass 

attached  to  the  shell  of  revolution  S. 

LSEGM(J)  =  segment  of  the  shell  S  in  which  the  J-th  point 
mass  is  attached. 

LMESH(J)  =  mesh  point  in  segment  LSEGM(J)  of  the  shell  S 
at  which  the  J-th  point  mass  is  attached. 

ANGPTM(J)  =  circumferential  coordinate  angle  0  (in  degrees) 
at  which  the  J-th  point  mass  is  attached. 

Note  1;  If  no  concentrated  masses  are  attached  to  the  shell  S, 
skip  Input  Set  8. 

Note  2;  Physical  points  on  the  shell  S  are  located  with  reference  to  a 
full  mathematical  model  in  accordance  with  the  conventions 
established  for  the  B0S0R4  code.  Thus,  to  locate  a  given 
physical  point  on  the  shell  S,  the  user  must  furnish  (!)■  the 
number  of  a  segment,  (2)  the  number  of  a  mesh  or  nodal  point 
within  the  given  segment,  and  (3)  a  circumferential  coordinate 
angle  (Fig.  1).  The  nodal  points  are  actually  points  on  the 
reference  meridian  used  to  generate  the  shell  of  revolution  S. 

They  are  numbered,  starting  from  unity,  within  each  segment,  and 
may  be  thought  of  as  meridional  stations. 

Note  3;  For  slde-on  or  oblique  loading,  the  origin  of  the  circumferential 
coordinate  angle  9  (or  equivalently,  the  direction  of  the 
positive  Z-axis)  is  determined  by  the  point  at  which  the  shock 
wave  first  contacts  the  shell  S,  as  indicated  in  Fig.l.  For 
end-on  loading,  in  which  the  shock  wave  first  contacts  the 
pole  of  an  end  closure,  the  origin  of  the  circumferential 
coordinate  angle  9  is  selected  by  the  user.  In  this  case,  a 
horizontal  Z-axis  is  the  recommended  configuration. 

Note  4;  The  number  of  logical  records  written  on  the  shell-fluid  file 
by  the  PIFLASH  code  (NSF)  is  given  by  the  expression 

NSF  =  29  +  5*NITEMS  +  2*NUM3ER  +  NEXTRA 


in  which  NEXTRA  *  1  if  NPTM  >  0 
NEXTRA  -  0  if  NPTM  <  0. 


A  sample  set  of  job  control  language  statements  or  control  cards  for 
an  execution  of  PIFLASH  under  a  NOS  1  operating  system  is  shown  below.  A 


brief  description  of  these  control  statements  is  presented  in  the  next 
subsection.  Various  other  arrangements  of  the  job  steps  are  possible. 


PIFLASH<T177)  RANLET 
USER (RANLETDfPASMORO) 

CHARGE (602201RtPROJECT) 
SETBSL<777) 

ATTACH (OLOPLaPIFLASH) 
UPDATE(F.LaA124) 

RETURN <OLDPL) 

FTN(I*COMPILEfL*OtOPT«l#A) 

RETURN(COHPILE) 

6ET(TAPE7sVIRMAS) 

ATTACH (TAPElQsSOSSYOO) 

ATTACH (TAPE11=BOSSYOI) 

ATTACH (TAPEl2aBOSSY02) 

ATTACH (TAPEl3aBOSSY03) 
FILE(TAPE99*RTaS*BT*C) 
PUR6E(SHLFLU/NA> 

DEFINE (TAPE99=SHLFLU/CT»S.M»R> 
RFL(70000) 

LDSET(PRESETslNOEF,MAP=SB) 
LOSET  <FILESaTAPE99) 

L60« 

GOTO (CATALOG) 

EXIT, 

CATAL0G(TAPE99tR) 


5.  Description  of  Control  Cards 


PIFLASH(T177)  RANLET 

Job  card  or  job  statement. 

USER (RANLETD , PASWORD) 

CHARGE (60220 IR , PROJECT) 

SETBSL(777) 

System  dependent  accounting  information. 

ATTACH (OLDPL=PIFLASH) 

Assign  permanent  file  PIFLASH  to  current  job  with  local  file  name  OLDPL. 

File  PIFLASH  contains  the  UPDATE  program  library  for  the  PIFLASH  code. 
UPDATE(F,L=A124) 

RETURN (OLDPL) 

Process  UPDATE  directives  supplied  in  the  job  stream  (on  file  INPUT) 
and  prepare  a  source  file  suitable  for  compilation. 

FTN (I=COMPILE , L=0 , 0PT= 1 , A) 

RETURN (COMPILE) 

Produce  compiled  code  or  object  module  for  the  PIFLASH  code. 

GET (TAPE 7=VIRMAS) 

Retrieve-  copy  of  Indirect  access  permanent  file  VIRMAS  for  use  as 
local  file  TAPE7.  File  VIRMAS  is  the  virtual  mass  file. 
ATTACH(TAPE10=BOSSY00) 

ATTACH (TAPE 1 1 =B0SSY0 1 ) 

ATTACH(TAPE12=BOSSY02) 

ATTACH (TAPE 13=BOSSY03) 

Assign  permanent  files  BOSSYOO,  BOSSYOl,  BOSSY02,  BOSSY03  to  current 
job  with  local  file  names  TAPE 10,  TAPEll,  TAPE 12,  TAPE13,  respectively. 

These  files  are  the  shell  mode  files. 

FILE(TAPE99 ,RT=S ,BT=C) 

Specify  the  record  type  of  local  file  TAPE99.  This  assures  that  the 
appropriate  system  routines  are  loaded  for  the  processing  of  local 
file  TAPE99.  For  the  system-logical  records  indicated  (RT=S) ,  each 
record  occupies  an  Integral  number  of  central  memory  words  and  is 
terminated  by  a  system-supplied  terminating  marker.  The  block  type 
indicated  (BT=C)  is  not  applicable  to  S-type  records,  but  is  provided 
for  compatibility  with  SC0PE2. 

PURGE (SHLFLU/NA) 

DEFINE (TAPE99=SHLFLU/CT=S ,M=R) 

Assign  space  on  permanent  file  device  for  the  shell-fluid  file  (TAPE99) 
under  the  permanent  file  name  SHLFLU.  Informaf.on  written  on  TAPE99 
during  the  execution  of  PIFLASH  is  written  directly  on  the  permanent  file. 
RFL( 70000) 

Request  field  length  for  execution.  This  is  optional  on  many  systems. 
LDSET(PRESET=INDEF,MAP=SB) 

LDSET(FILES-TAPE99) 

LGO. 

Select  load  options,  load,  and  execute  the  compiled  version  of  PIFLASH. 

The  required  set  of  input  data  is  supplied  in  the  job  stream  (on  file  INPUT). 
The  FILES  parameter  appearing  on  the  LDSET  control  statement  above  must 
be  supplied  for  each  file  name  referenced  on  a  FILE  control  statement. 


GOTO (CATALOG) 

Unconditional  transfer  to  control  statement  CATALOG  if  the  execution 
terminates  normally. 


EXIT. 

Terminate  execution  or 
error  occurs. 

CATALOG (TAPE99,R) 

List  information  about 


provide  transfer  of  control  when  an  execution 
each  logical  record  of  file  TAPE99  on  file  OUTPUT. 
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UPDATE  Directives  and  Input  Data  for  Sample  Problem  1 


The  UPDATE  directives  employed  to  prepare  the  shell-fluid  file  for 


Sample  Problem  1  are  listed  below. 


*IDENT, STORAGE 

*DELETE .PIFLASH . 22 .PIFLASH . 23 
COMMON  A (10000) 

NWORDA  =  10000 


The  size  of  the  dynamic  array  A  shown  above  is  in  excess  of  the  6128  words 
actually  required  for  Sample  Problem  1.  In  general,  a  dynamic  array  dimen¬ 
sioned  to  25000  words,  corresponding  to  a  central  memory  requirement  of 
about  125000  words  (octal),  is  sufficient  for  most  problems.  A  field  length 
of  almost  70000  words  (octal)  is  required  for  an  execution  of  PIFLASH  when 
the  dynamic  array  is  dimensioned  as  above. 

The  set  of  input  data  shown  below  directs  the  PIFLASH  code  to  prepare 

% 

'  the  shell-fluid  file  for  Sample  Problem  1. 


4 

0 

4 

1 

1 

NUMBER.NT08SN t NPTM 

I.NSYMStNSYNR 

2 

2 

2 

2 

(NMETS8(K>»NUSES0(K)»K»1 

1*2) 

386.4 

eoavTY 

10 

11 

12 

13 

INT4PE(J>*U«l»NirENS> 

25 

30 

30 

30 

(N JUSE IJ) • 1 • NI TEMS ) 

1 

2 

1 K0RS8 ( K ) tK> 1 t NKORSe ) 

2 

4 

6 

T 

8 

9 

11 

13 

14  15 

16 

IT 

18 

19 

20 

JUSE 

21 

22 

23 

24 

25 

26 

2T 

28 

29  30 

JUSE 
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2 
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4 

5 

6 

T 

8 

9  10 

11 

12 

13 

14 

15 

JUSE 

16 

IT 

18 

19 

20 

21 

22 

23 

24  25 

26 

2T 

28 

29 

30 

JUSE 

1 

2 

3 

4 

5 

6 

T 

8 

9  10 

11 

12 

13 

14 

15 

JUSE 

16 

IT 

18 

19 

20 

21 

22 

23 

24  25 

26 

2T 

28 

29 

30 

JUSE 

1 

2 

3 

4 

5 

6 

T 

8 
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11 

12 

13 

14 

15 

JUSE 

16 

IT 

18 

19 

20 

21 

22 

23 

24  25 

26 

2T 

28 

29 

30 

JUSE 

32.0 

2 

12 

0.0 

MEI6NT<ll*LSEeNll) 

•LMESHIll 

«AN0PTM(1> 

32.0 

2 

12 

100.0 

«EieHT<2)*t.SE8M(2l  .LMESH(2) 

*AH6PTM(2> 

30.5 

2 

53 

90.0 

MeieHT<3) 

.I.SE8NI3I 

*LMESH(3) 

.ANOPTMO) 

30.5 

2 

53 

2T0.0 

«EI6HT(4).LSE6M(4) 

•LMESHU) 

tAN6PTM(4) 
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7.  UPDATE  Directives  and  Input  Data  for  Sample  Problem  2 


The  UPDATE  directives  employed  to  prepare  the  shell-fluid  file  for 
Sample  Problem  2  are  listed  below. 


*IDENT, STORAGE 

*DELETE .PIFLASH . 22 .PIFLASH . 23 
COMMON  A (10000) 

NWORDA  =  10000 


The  size  of  the  dynamic  array  A  shown  above  is  somewhat  in  excess  of  the 


size  actually  required  for  Sample  Problem  2 .  When  the  dynamic  array  is 
dimensioned  to  10000  words,  as  above,  a  field  length  of  almost  70000 
words  (octal)  is  required  for  an  execution  of  the  PIFLASH  code. 


The  set  of  input  data  shown  below  directs  the  PIFLASH  code  to  prepare 
the  shell-fluid  file  for  Sample  Problem  2. 


A 

0 

4 

1 

1 

3 

3 

1 

1 

3SA.4 

10 

11 

12 

13 

IS 

30 

15 

IS 

3 

2 

3 

5 

S 

7 

1 

2 

3 

4 

5 

IS 

17 

15 

19 

20 

21 
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2 

3 

4 

5 

1 

2 

3 

4 

5 

32.0 

2 

11 

0.0 

32.0 

2 

11 

100.0 

30.5 

3 

19 

90.0 

30. S 

3 

19 

270.0 

NUMBERtHTORSN«NPTM»NSYMStNSYNP 

fNHETSeiK)*NUSESe(K)»KBl«2) 

eRAVTY 

<NTAPEfJ)»J«l*NirrNS) 

<NJUSEIJ) lENS) 

(KORSe(K) tKslfNKORSe) 


9 

10 

11 

12 

14 

IS 

IS 

17 

10 

JUSE 

7 

0 

9 

10 

11 

12 

13 

14 

15 

JUSE 

22 

23 

24 

25 

2S 

27 

20 

29 

30 

JUSE 

7 

0 

9 

10 

11 

12 

13 

14 

15 

JUSE 

7 

0 

9 

10 

11 

12 

13 

14 

15 

JUSE 

MEISHT ( 1 I tLSESR I 1 ) »LMESH ( J ) , AN6RTM ( 1 ) 
«EieHTl2)tLSE6i«l2)*LMESM(2)«ANSRTM(2> 
KEISHTISI  »LSE6RO>  tLMESHO)  «  ANOPTMCS) 
VCIBHTU)  tLSEBRU)  tLMESHCAl  tANePTMCA) 
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APPENDIX  D  -  THE  USE  OF  SAP IV 


1.  Background 

The  SAPIV  finite-element  code  of  Ref.  [8]  is  included  in  the  ELSHOK 
suite  of  computer  programs  for  the  modal  analysis  of  the  internal 
appendages  or  substructures  a.  The  general-purpose  SAPIV  code  is  based 
upon  a  lumped  mass  formulation  and  applies  to  linearly  elastic  structures 
Comprehensive  documentation  of  SAPIV  may  be  found  in  Ref.  [8],  A 
knowledge  of  the  Information  given  in  Ref.  [3]  and  Sections  II  and  III  of 
the  present  report  will  aid  in  the  understanding  of  this  appendix. 

A  number  of  modifications  have  been  made  to  SAPIV  to  enable  it  to 
function  within  the  scope  of  an  ELSHOK  calculation.  Thus,  in  addition 
to  determining  the  fixed-base  modes  and  corresponding  natural  frequencies 
of  a  substructure,  the  unconstrained  or  free-free  physical  mass  and 
stiffness  matrices  (M  and  K,  respectively)  are  also  determined.  These 
matrices  are  required  for  the  calculation  of  the  constraint  modes  gj^ 
of  Eq.(lO)  and  of  the  coefficients  needed  to  evaluate  the  interaction 

forces  and  moments  developed  between  each  substructure  a  and  the  main 
body  S.  It  should  be  pointed  out  that  the  entire  physical  mass  and 
stiffness  matrices  of  a  substructure  are  not  required  for  the  subsequent 
submerged  shock  response  problem.  In  fact,  as  discussed  in  Ref.  [3], 
the  solution  of  this  response  problem  requires  only  those  rows  of  K  and 
and  only  those  main  diagonal  terms  of  M  which  correspond  to  the  physical 


degrees  of  freedom  of  a  constrained  to  move  with  S. 


2.  Rules  of  Operation 


For  compatibility  with  the  current  version  of  the  ELSHOK  code,  the 

following  rules  must  be  observed  when  executing  the  SAPIV  code: 

(1)  The  rectangular  Cartesian  coordinate  system  (x,y,z)  used  to  describe 

a  substructure  a  must  be  a  right-handed  system,  and  the  x-axis  of  this 
system  of  coordinates  must  coincide  in  direction  and  sense  with  the 
longitudinal  or  X-axis  of  the  main  body  S. 

(2)  If  a  substructure  is  symmetrical  about  the  longitudinal  midpoint  of  S, 
in  a  problem  for  which  symmetry  is  being  used,  the  user  must  model  only 
the  left  half  of  the  substructure.  If  symmetry  is  not  applicable,  the 
entire  model  of  such  a  substructure  must  be  provided. 

(3)  The  input  data  supplied  for  SAPIV  in  an  ELSHOK  application  must  be 
provided  in  two  parts.  The  first  part  is  standard  input  prepared  in 
accordance  with  the  requirements  of  Ref.  [8],  except  that  (1)  additional 
input  items  must  be  supplied  on  the  master  control  and  dynamic  analysis 
input  cards  and  (2)  the  physical  degrees  of  freedom  associated  with  the 
supports  or  bases  of  the  substructure  under  study  must  not  be  constrained 
The  second  part,  appended  to  the  first  part  without  a  separator  of  any 
kind,  supplies  the  fixity  allowing  the  calculation  of  the  fixed-base 
modes  and  corresponding  natural  frequencies. 

(4)  A  separate  execution  of  SAPIV  must  be  made  for  each  substructure 
Included  in  an  analysis. 


3.  Modifications  to  SAPIV 

A  number  of  changes  have  been  made  to  SAPIV  to  facilitate  Its  use 
in  the  ELSHOK  suite  of  computer  programs.  A  brief  description  of  each 
of  these  modifications  follows: 

(1)  A  file  has  been  supplied  and  assigned  to  logical  unit  98.  This  file 
is  referred  to  by  the  local  file  name  TAPE98  and  serves  as  a  scratch 
file. 

(2)  A  file  has  been  supplied  and  assigned  to  logical  unit  99.  This  file 
is  referred  to  by  the  local  file  name  TAPE99  and  serves  as  the  sub¬ 
structure  mode  file  of  Section  III.  File  TAPE99  contains  geometry, 
natural  frequencies,  fixed-base  modes,  and  the  unconstrained  physical 
mass  and  stiffness  matrices,  in  unformatted  binary  records,  for  a 
given  substructure.  The  substructure  mode  file  must  be  saved  after 
each  execution  of  SAPIV  to  provide  data  needed  subsequently. 

(3)  Although  not  a  modification  to  SAPIV,  it  should  be  pointed  out  here 
that  a  small  plotting  processor  (PLOTSAM)  is  supplied  with  ELSHOK. 
This  processor  prepares  plots  on  a  TEKTRONIX  graphics  terminal  and 
may  be  used  to  plot  the  finite-element  model  and  the  fixed-base 
modes  of  a  given  substructure. 

(4)  Additional  input  items  have  been  added  to  the  usual  input  data. 

These  items  were  mentioned  in  the  preceding  subsection  and  will 
be  described  in  detail  in  the  next  subsection. 

The  above  modifications  have  been  made  through  the  use  of  the  UPDATE 
batch  editing  system  (Ref.  [9])  and  may  be  easily  removed  to  return  SAPIV 
to  its  standard  configuration. 
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4.  Additional  Input 

When  using  SAPIV  with  the  ELSHOK  code,  all  of  the  applicable  input 
data  described  in  Ref.  [8]  must  be  prepared  for  the  problem  of  interest 
in  accordance  with  the  above  rules  of  operation.  Thus,  two  additional 
input  parameters  (NRIGID  and  SHIFT)  must  be  supplied  within  the 
standard  input  data.  Descriptions  of  these  parameters  are  given  below. 
Following  Appendix  B,  the  input  data  specifications  are  written  in  a 
style  similar  to  FORTRAN. 

For  ELSHOK  applications,  the  master  control  card  in  the  standard 
SAPIV  input  data  (Input  Section  II  of  Ref.  [8])  has  been  changed  to 

Master  Control  Card  FORMAT  (1015) 

NUMNP,  NELTYP,  LL,  NF,  NDYN, 

MODEX,  NAD,  KEQB,  NIOSV,  NRIGID 

where  NUMNP  »  total  number  of  nodal  points  in  the  model. 

NELTYP  =  number  of  element  groups  or  types. 

LL  =  number  of  loading  conditions.  For  ELSHOK  applications, 
use  LL  =  0. 

NF  =  number  of  frequencies  to  be  found. 

NDYN  =  analysis  type  code.  For  ELSHOK  applications,  use  NDYN  =  1. 

MODEX  =  program  execution  mode.  For  ELSHOK  applications,  use 
use  MODEX  =  0. 

NAD  =  total  number  of  vectors  to  be  used  in  a  subspace 

iteration  solution  for  eigenvalues  and  eigenvectors.  For 
ELSHOK  applications,  set  NAD  =  0  and  accept  the 
default  value. 

KEQB  =  number  of  degrees  of  freedom  (equations)  per  block  of 
storage.  For  ELSHOK  applications,  set  KEQB  =  0  and 
accept  the  value  calculated  automatically  by  SAPIV. 

NIOSV  =  integer  flag  indicating  whether  or  not  stress  information 
is  to  be  written  on  local  file  TAPEIO  for  post-processing. 

For  ELSHOK  applications,  use  NIOSV  =  0. 

NRIGID  =  integer  flag  indicating  whether  or  not  free-free  modes, 

including  rigid  body  modes,  are  to  be  determined  (l=yes;  0=no) 
For  the  current  version  of  ELSHOK,  which  employs  fixed-base 
modes,  use  NRIGID  =  0. 
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The  dynamic  analysis  card  used  to  set  parameters  for  the  calculation  of 
mode  shapes  and  frequencies  In  the  standard  SAFIV  Input  data  (Input  Section 


VII  A  of  Ref.  [8])  has  been  changed  to 


lamlc  Analysis  Card 


FORMAT  (3I5,2E10.0,I5,E10.0) 


IFPR,  IFSS,  NITEM,  RTOL,  COFQ,  NFO,  SHIFT 


where  IFPR  =  flag  for  printing  Intermediate  output  (l=yes;  0=no) . 

For  ELSHOK  applications,  use  IFPR  =®  0. 

IFSS  =  flag  for  by-passing  the  Sturm  sequence  check  in  an  eigenvalue 
calculation  (l*yes;  0=no) . 

NITEM  =  maximum  number  of  iterations  allowed  for  reaching  the  con¬ 
vergence  tolerance.  For  ELSHOK  applications,  set  NITEM  =  0 
and  accept  the  default  value. 

RTOL  =  convergence  tolerance  for  the  highest  frequency  computed. 

For  ELSHOK  applications,  set  RTOL  =0.0  and  accept  the 
default  value. 

COFQ  =  cutoff  frequency  in,  cycles  per  unit  time.  For  ELSHOK 

applications,  set  COFQ  =  0.0  to  allow  the  computation  of 
NF  frequencies. 

NFO  =  number  of  starting  iteration  vectors.  For  ELSHOK  applications, 
use  NFO  =  0. 

SHIFT  *  frequency  shift  used  in  the  calculation  of  free-free  modes 
(NRIGID  =  1) .  For  the  current  version  of  ELSHOK,  use 
SHIFT  =0.0. 


The  set  of  input  data  supplied  to  provide  the  fixity  needed  for  the 


calculation  of  the  fixed-base  modes  and  corresponding  natural  frequencies 


must  be  placed  immediately  after  the  dynamic  analysis  card  described  above. 


A  description  of  this  additional  input  follows: 
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FORMAT  (215) 


Input  Set  1 

NFACE,  NSYMTY 

where  NFACE  =  number  of  Interface  or  attachment  points  between  the  sub¬ 
structure  O  under  study  and  the  shell  S. 

NSYMTY  =  Integer  flag  Indicating  whether  or  not  symmetry  is  employed 
in  the  finite  element  model  of  the  substructure  (l=yes; 
0=no) .  For  the  current  version  of  ELSHOK,  only  symmetry 
associated  with  the  x-direction  (i.e.,  about  a  y-z  plane) 
may  be  Incorporated  into  the  model. 

Note  1  :  For  the  current  version  of  ELSHOK,  symmetry  may  not  be  employed 
unless  the  entire  structure  under  consideration  (including  the 
shell  S,  any  internal  appendages  a,  and  any  concentrated  masses) 
is  symmetrical  about  the  longitudinal  midpoint  on  the  main  body  S 

Note  2  :  When  NSYMTY  =  1,  all  effects  of  a  substructure  on  the  shell  S  are 
doubled  for  use  in  the  equations  governing  the  submerged  shock 
response  problem.  Thus,  in  a  problem  for  which  symmetry  is 
employed,  identical  substructures,  one  in  the  left  half  of  S  and 
one  in  the  right,  may  be  accounted  for  by  modeling  the  left-hand 
structure  without  the  use  of  symmetry  and  by  setting  NSYMTY  =  1. 

Input  Set  2  FORMAT  (715) 

DO  40  K  -  1,  NFACE 

40  READ  NODE,  (NFIXTY(J),  J  =  1,  6) 

where  NODE  ■  identification  number  of  a  nodal  point  serving  as  an 

Interface  or  attachment  point  between  the  substructure  a 
under  study  and  the  shell  S. 

NFIXTY  ■  array  containing  the  boundary  condition  codes  (1  or  0)  for 
nodal  point  NODE.  Elements  1-6  of  this  array  correspond 
to  the  x-translation,  y-translation,  z-translation, 
x-rotation,  y-rotation,  and  z-rotation  boundary  condition 
codes,  respectively,  for  the  indicated  nodal  point. 

NFIXTY (J)  ■  1;  physical  degree  of  freedom  J  is  constrained 

to  move  with  S  or  is  not  active  in  the  problem. 

NFIXTY (J)  ■  0;  physical  degree  of  freedom  J  is  not  con¬ 

strained  to  move  with  S  but  is  active. 

Note  ;  In  the  terminology  of  Ref. [8],  slave  degrees  of  freedom  may 
not  be  constrained  to  move  with  the  main  body  S. 

The  two  sets  of  input  data  described  above  must  be  followed  by  the 

usual  blank  or  dummy  title  and  master  control  cards  used  to  terminate  an 


execution  of  SAPIV. 


SAPIV(T177l  RANLET 
USER (RANLE70tPASW0R0) 

CHARGE (602201R«PROJECT) 

SETBSL(777) 

ATTACH<LGOaSAPFTN) 

FlLE(TAPEl*BT*C*R7aS) 

FILE<TAPE2*BTaC*R7aS) 

FILE (TAPE3*BT*C*RT*Sl 
FlLE(TAPE4*BTaC*R7aS) 

FILEe7APE7fB7aC#R7aS) 

FILE(TAPEa#BTaC#RTaS) 

FILE{TAPE9#BTaC#R7aS) 

FILE(TAPE10tBTaC*RTaS) 

FILE(7APE98«BTaCtR7aS) 

FlLE<TAPE99,B7*C.R7aS) 

PURGE (SUBMODE/NA) 

define (7APE99=SUBM00E/CTaStMsR) 

RFL(170000) 

L0SET<PRESE7aIN0EF*MAPsN) 

LDSE7(FILES*TAPE1/TAPE2/TAPE3/7APE4/TAPE7J 

LDSET(FILESaTAPE8/TAPE9/TAPE10/TAPE98/TAPE99) 

LGO, 

GOTO(CATALOG) 

EXIT. 

CATAL0G(TAPE99fR> 
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6.  Description  of  Control  Cards 

SAPIV(T177)  RANLET 

Job  card  or  job  statement. 

USER (RANLETD , PASWORD) 

CHARGE (60220 IR, PROJECT) 

SETBSL(777) 

System  dependent  accounting  information. 

ATTACH (LGO*SAPFTN) 

Assign  permanent  file  SAPFTN  to  current  job  with  local  file  name  LGO. 

File  SAPFTN  contains  the  compiled  code  or  object  module  for  SAPIV. 

FILE (TAPE 1 ,BT=C,RT=S) 

FILE (TAPE2 ,BT=C ,RT=S) 

FILE (TAPES , BT=C , RT=S ) 

FILE(TAPE4 ,BT=C ,RT=S) 

FILE (TAPE  7 , BT-C , RT=S ) 

FILE (TAPES , BT-C , RT-S) 

FILE (TAPES , BT=C , RT=S) 

FILE (TAPE 10 , BT=C , RT=S) 

FILE (TAPE98 , BT=C , RT=S) 

FILE (TAPE99 , BT=C , RT=S) 

Specify  the  record  type  of  each  of  the  indicated  local  files.  This 
assures  that  appropriate  system  routines  are  loaded  for  processing. 

PURGE (SUBMODE/NA) 

DEFINE(TAPE99=SUBM0DE/CT=S ,M=R) 

Assign  space  on  permanent  file  device  for  the  substructure  mode  file 
(TAPE99)  under  the  permanent  file  name  SUBMODE.  Information  written  on 
TAPE99  during  the  execution  of  SAPIV  is  written  directly  on  the 
permanent  file. 

RFL( 170000) 

Request  field  length  for  execution.  This  is  optional  on  many  systems. 
LDSET (PRESET=INDEF ,MAP=N) 

LDSET (FILES=TAPE 1 /TAPE2 /TAPES /TAPE4 /TAPE 7 ) 
LDSET(FILES=TAPE8/TAPE9/TAPE10/TAPE98/TAPE99) 

LGO. 

Select  load  options,  load,  and  execute  the  compiled  version  of  SAPIV.  The 
required  set  of  input  data  is  supplied  in  the  job  stream  (on  file  INPUT). 
GOTO (CATALOG) 

Unconditional  transfer  to  control  statement  CATALOG  if  the  execution 
terminates  normally. 

EXIT. 

Terminate  execution  or  provide  transfer  of  control  when  an  execution 
error  occurs. 

CATALOG (TAPE99,R) 

List  information  about  each  logical  record  of  file  TAPE99  on  file  OUTPUT. 
The  use  of  this  control  statement  is  not  recommended  for  finite  element 
calculations  involving  more  than  about  100  physical  degrees  of  freedom. 


.V.v^V. 
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7 •  Input  Data  for  Sample  Problem  1 

The  set  of  input  data  shovm  below  directs  the  SAPIV  code,  as  modified 
for  use  with  ELSHOK,  to  perform  the  modal  analysis  of  the  substructure 
included  in  Sample  Problem  1 . 


OSVtNRieiO 


SYMMETRICAL 

SUBSTRUCTURE  - 

-  05/15/02 

25 

1 

0 

S 

1 

0 

0 

0 

0 

0  NUMNFtNELTYRt. 

1 

1 

1 

0 

1 

1 

1 

20.0 

0.0 

•16.71675 

s 

1 

I 

0 

1 

1 

1 

20.0 

0.0 

-3.343750 

6 

1 

1 

0 

1 

0 

1 

0.0 

0.0 

0.0 

IS 

1 

1 

0 

1 

0 

1 

20.0 

0.0 

0.0 

16 

1 

1 

0 

1 

0 

1 

22. 

4S125 

0.0 

0.0 

24 

1 

1 

0 

1 

0 

1 

42. 

33125 

0.0 

0.0 

2S 

1 

1 

0 

1 

1 

1 

44 

.6125 

0.0 

0.0 

2 

24 

2 

0 

2 

1 

29.  OE 

06 

o.a 

7.29C«04 

2 

29.  OE 

06 

0.3 

0.0 

1 

0.0 

Q.O 

0.0 

1.0 

2.6666667 

10.666667 

2 

6.0 

0.0 

0.0 

1.0 

16.0 

O.S 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1 

1 

2 

25 

2 

2 

4 

4 

5 

25 

2 

2 

S 

S 

15 

25 

2 

2 

000001 

6 

6 

7 

1 

1 

1 

24 

24 

2S 

1 

1 

1 

MULT  e  X 
MULT  6  Y 
MULT  e  Z 


1 

1 


0 

.0 

0 

1 

1 


CMO  OF  DATA 
0  0 


0.0  o.o 

0.0  0.0 


<OUMMY  TITLE) 
0  0 


0.0  CONC  MASS 


0.0  0.0  0.0 

0.0  CLEMENT  LOAD  MULTIFLtCMS 
•••  •  •••  IFFR.IFSSt...*NrO.SHlFT 

NFACCvNSYMTY 

1  NOOC*(NFlXTV(J).Jal.O) 
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The  set  of  input  data  shovm  below  directs  the  SAPIV  code,  as  modified 


for  use  with  ELSHOK,  to  perform  the  modal  analysis  of  the  substructure 
included  in  Sample  Problem  2. 


SYMMETRICAL 

SUBSTRUCTURE  - 

-  05/15/82 

25 

1 

0 

5 

1 

0 

0 

0 

0 

0  NUMNP»NELTYPt...*N10SV«NRI6IO 

1 

1 

1 

0 

1 

1 

1 

20.0 

0.0  -16.71875  0 

5 

1 

1 

0 

1 

1 

1 

20.0 

0.0  -3.343750  1 

6 

1 

I 

0 

1 

0 

1 

0.0 

0.0  0.0  0 

15 

1 

1 

0 

1 

0 

1 

20.0 

0.0  0.0  1 

16 

1 

1 

0 

1 

0 

1 

22. 

48125 

0.0  0.0  0 

24 

1 

1 

0 

1 

0 

1 

42. 

33125 

0.0  0.0  1 

25 

1 

1 

0 

1 

1 

1 

44 

.8125 

0.0  0.0  0 

2 

24 

2 

0 

2 

1 

29.  OE 

06 

0.3 

7.29E-04 

2 

29.  OE 

06 

0.3 

0.0 

1 

8.0 

0.0 

0.0 

1.0 

2.6666667  10.666667 

2 

6.0 

0.0 

0.0 

1.0 

18.0  0.5 

0.0 

0.0 

0.0 

0.0 

MULT  6  X 

0.0 

0.0 

0.0 

0.0 

MULT  8  Y 

0.0 

0.0 

0.0 

0.0 

MULT  8  Z 

1 

1 

2 

25 

2 

2 

4 

4 

5 

25 

2 

2 

5 

5 

15 

25 

2 

2 

000001 

6 

6 

7 

1 

1 

1 

24 

24 

25 

1 

1 

1 

0 

0 

0.0 

0.0 

0.0 

0.0  0.0  0.0 

CONC  MASS 

0.0 

0.0 

0.0 

0.0 

ELEMENT  LOAD  MULTIPLIERS 

0 

0 

0 

0.0 

0.0 

0 

0.0  IFPR.IFSSf ...» 

NFO.SHIFT 

1 

1 

NFACEtNSYMTY 

1 

1 

1 

1 

1 

1 

1 

N00E»(NPIXTY(J)tJBlt6) 

ENO  OF  DATA 

(DUMMY  TITLEI 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

The  above  set  of  input  data  is  identical  to  that  employed  for  Sample 
Problem  1.  Thus,  the  execution  of  SAPIV  for  Sample  Problem  2  may  be 
skipped  if  the  substructure  mode  file  for  Sample  Problem  1  is  available 
for  the  processing  required  in  the  solution  of  Sample  Problem  2. 


APPENDIX  E  -  THE  USE  OF  PICRUST 


1.  General  Information 

The  PICRUST  computer  code  is  a  post-processor  for  the  SAPIV  code.  As 
discussed  in  Section  III,  the  PICRUST  code  takes  data  from  a  substructure 
mode  file  (produced  by  SAPIV),  adds  to  it,  and  reorganizes  it  to  facilitate 
the  solution  of  the  equations  governing  the  transient  response  problem  under 
study.  During  an  execution  of  PICRUST,  connectivity  between  a  substructure  cr 
and  the  shell  S  is  accounted  for,  constraint  modes  are  computed,  and  coefficients 
for  evaluating  the  forces  developed  at  the  points  of  attachment  between  C 
and  S  are  determined.  Familiarity  with  Appendix  D  and  Ref.  [3]  will  aid 
in  the  understanding  of  this  appendix. 

2.  Rules  of  Operation 

(1)  A  substructure  mode  file  from  an  execution  of  SAPIV  must  be  available 
prior  to  the  execution  of  the  PICRUST  code. 

(2)  A  separate  execution  of  PICRUST  is  required  for  each  substructure 
included  in  the  problem  under  consideration. 

(3)  The  substructure  mode  file  must  be  supplied  to  PICRUST  as  a  local  file 
assigned  to  logical  unit  4  and  referred  to  as  TAPE4 . 

(4)  For  ease  of  merging  files  to  complete  the  input  file  required  for  the 
submerged  shock  response  problem,  the  user  should  supply  control  cards 
or  control  statements  to  ensure  that  the  substructure  file  produced 

by  PICRUST  (TAPEIO  on  logical  unit  10)  is  comprised  of  system-logical 
records. 

(5)  File  TAPEIO  must  be  saved  to  provide  data  needed  subsequently. 


3.  Description  of  Input  Data 


The  input  data  specifications  required  for  an  execution  of  the  PICRUST 
code  are  described  below.  Following  Appendix  B,  these  specifications  are 
written  in  a  style  similar  to  FORTRAN. 

Input  Set  1  FORMAT  (1015) 

(LISTIT(J),  J  =  1,  20) 

where  LISTIT(l)  =  output  flag  for  boundary  condition  codes  of  free-free 

or  unconstrained  substructure . 

LISTIT(2)  =  output  flag  for  equation  numbers  or  problem  unknowns 
associated  with  free-free  substructure. 

LISTIT(3)  =  output  flag  for  nodal  point  geometry. 

LISTIT(4)  =  output  flag  for  main  diagonal  terms  of  diagonal  mass 
matrix  of  free-free  substructure. 

LISTIT(5)  =  output  flag  for  stiffness  matrix  of  free-free  sub¬ 
structure. 

LISTIT(6)  =  output  flag  for  connectivity  data  for  interface  or  attach¬ 
ment  points  between  current  substructure  and  main  body  S. 

LISTIT(7)  =  output  flag  for  equation  numbers  of  free-free  substructure 

after  partitioning  of  physical  degrees  of  freedom  to  corres¬ 
pond  to  Eqs.  (9)  and  (10). 

LISTIT(8)  =  output  flag  for  main  diagonal  terms  of  diagonal  mass  matrix 
of  free-free  substructure  after  partitioning. 

LISTIT(9)  =  output  flag  for  stiffness  matrix  of  free-free  substructure 
after  partitioning. 

LISTIT(IO)  =  output  flag  for  constraint  modes  in  blocking  (by  rows) 
used  for  their  solution. 

LISTIT(ll)  =  output  flag  for  constraint  modes  in  blocking  (by  columns) 
used  to  display  all  physical  degrees  of  freedom  in  a  given 
constraint  mode. 

LISTIT(12)  =  output  flag  for  fixed-base  natural  frequencies. 

LISTIT(13)  =  output  flag  for  fixed-base  modes  displayed  with  status 
of  six  degrees  of  freedom  at  each  node  shown. 

LISTIT(14)  =  output  flag  for  fixed-base  modes  displayed  in  compacted 
form  (i.e.,  with  only  active  degrees  of  freedom  shown). 

LISTIT(15)  =  output  flag  for  expansion  coefficients  of  a  modal  series 
representation,  in  terms  of  fixed-base  modes  of  substruc¬ 
ture,  of  constraint  modes. 

LISTIT(16)  =  output  flag  for  skyline  bandwidth  of  those  rows  of  parti¬ 
tioned  free-free  stiffness  matrix  corresponding  to  physi¬ 
cal  degrees  of  freedom  of  substructure  constrained  to 
move  with  shell  S. 
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LISTIT(17)  =  output  flag  for  portion  of  above  rows  of  partitioned 
free-free  stiffness  matrix  [LISTIT(16)]  corresponding 
to  physical  degrees  of  freedom  of  substructure  not 
constrained  to  move  with  shell  S. 

LISTIT(18)  =  output  flag  for  interface  force  coefficients  associated 
with  fixed-base  modes  of  substructure. 

LISTIT(19)  =  output  flag  for  portion  of  above  rows  of  partitioned 
free-free  stiffness  matrix  tLISTIT(16)]  corresponding 
to  physical  degrees  of  freedom  of  substructure  constrained 
to  move  with  shell  S. 

LISTIT(20)  =  output  flag  for  Interface  force  coefficients  associated 
with  constraint  modes  of  substructure. 

Note  1;  If  an  output  flag  is  set  to  1,  output  is  supplied.  If  an  output 
flag  is  set  to  0,  no  output  is  supplied. 

Note  2:  In  almost  all  cases,  use  LISTIT(5)  =  LISTIT(9)  =  LISTIT(17)  = 

LISTIT(19)  =  0,  since  the  structural  stiffness  matrix  is  a  very 
large  array  for  most  problems  of  interest. 


Input  Set  2  FORMAT  (315) 

NJUSE,  NHWSOB,  NHWBAR 

where  NJUSE  =  number  of  fixed-base  modes  of  substructure  to  be  used 

for  subsequent  response  calculation.  The  value  of  NJUSE 
must  not  exceed  the  value  of  NF  used  in  the  SAPIV  calculation. 

NHWSOB  =  number  of  hectowords  defining  size  of  storage  blocks  nsed 
in  caxculation  of  constraint  modes  of  substructure. 

NHWSOB  <  0:  constraint  modes  computed  in  one  block  (entire 
array  in  core),  if  possible. 

NHWSOB  >0;  constraint  modes  computed  in  blocks  of  specified 
or  smaller  size,  depending  on  available  storage. 

NHWBAR  =  number  of  hectowords  defining  size  of  storage  blocks  used 

for  interface  force  coefficients  associated  with  constraint 
modes  of  substructure.  For  the  current  version  of  ELSHOK, 
use  NHWBAR  =  0.  This  specification  directs  the  PICRUST 
code  to  write  the  entire  array  of  coefficients  on  the  sub¬ 
structure  file  in  one  block. 

Note ;  If  the  size  of  the  constraint  mode  array  is  smaller  than  about 

10000  words,  the  use  of  NHWSOB  =  0  is  recommended.  The  constraint 
mode  array  has  NUDOF  rows  and  NCDOF  columns,  where  NUDOF  is  the 
number  of  physical  degrees  of  freedom  of  o  not  constrained  to 
move  with  S  and  NCDOF  is  the  number  of  physical  degrees  of  free¬ 
dom  of  a  constrained  to  move  with  S. 


108 


INPUT  Set  3 


FORMAT  (3I5,5X,E10.0) 


DO  50  K  =  1,  NFACE 

50  READ  NIPSUB,  LBOSEG,  LBOSPT,  ANGDEG 

where  NIPSUB  =  node  number  of  substructure  a  which  is  connected  to  the 

shell  S  at  the  K-th  interface  or  attachment  point  between 
a  and  S. 

LBOSEG  =  segment  of  the  shell  S  in  which  the  K-th  interface  point 
is  located. 

LBOSPT  =  mesh  point  in  segment  LBOSEG  of  the  shell  S  at  which 
the  K-th  connection  between  a  and  S  is  made. 

ANGDEG  =  circumferential  coordinate  angle  9  (in  degrees)  at 
which  the  K-th  interface  point  is  located  (Fig.l). 

Note  1  :  The  number  of  interface  points  (NFACE)  is  retrieved  from  the 
substructure  mode  file  and  need  not  be  supplied  by  the  user 
during  an  execution  of  PICRUST. 

Note  2;  The  individual  sets  of  the  above  four  connection  parameters 
may  be  supplied  to  the  PICRUST  code  in  any  convenient  order. 

Note  3;  Physical  points  on  the  shell  S  are  located  with  reference  to 
a  full  mathematical  model  in  accordance  with  the  conventions 
established  for  the  B0S0R4  code. 

Input  Set  4  FORMAT  (ElO.O) 

DEGROT 

where  DEGROT  =  orientation  angle  (in  degrees)  of  coordinate 
axes  of  substructure  a  with  respect  to  global 
coordinate  axes  of  shell  S.  The  angle  DEGROT  (a^  of  Fig.  1) 
defines  the  angle  between  the  positive  z-axis  of  a  and 
the  positive  Z-axis  of  S.  Positive  angles  are  measured 
counterclockwise  (when  looking  in  the  positive  X-direction) 
from  the  Z-axis  of  S  to  the  z-axis  of  o. 


FORMAT  (1515) 


Input  Set  5 

(JUSE(J),  J  =  1,  NJUSE) 

where  JUSE  =  array  containing  the  identification  numbers  of  the  fixed- 
base  modes  of  substructure  o  selected  for  use  by  the  user. 

Note:  The  number  of  logical  records  written  on  the  substructure  file 
by  the  PICRUST  code  (NSB)  is  given  by  the  expression 

NSB  =  15  +  NBLOCK  +  NBLFOR 


in  which  NBLOCK  is  the  number  of  storage  blocks  needed  for  the 
constraint  modes  and  NBLFOR  =  1  for  the  current  version  of  ELSHOK 


4.  Sample  Control  Cards 


A  sample  set  of  job  control  cards  or  job  control  language  statements 
for  an  execution  of  PICRUST  under  a  NOS  1  operating  system  is  shown  below. 
The  next  subsection  contains  a  brief  description  of  these  control  statements 
Various  other  arrangements  of  the  job  steps  are  possible. 


PICRUST (T277)  RANLET 
USER (RANLETOvPASMORD) 

CHARGE (602201RtPROJECT) 

SETBSL(777) 

ATTACH ( 0LDPL*P ICRUST ) 

UPDATE (F»L»A124) 

RETURN<OLOPL) 

FTN(I=COHPILE#L«0#OPT*l*PLs25000#A) 
RETURN (COMPILE) 

FILE(TAPE4.8TaC.RT«S) 

FILE<TAPE6#0TaC*RTaS) 

FILE(TAPE7.BT=C.RT*S) 

FlLE(TAPEa.BTsCtRT*S) 

FILE(TAPE9*BT*C.RT*S) 

FILE(TAPE10.BTaCfRT*S) 

ATTACH <TAPE4*SUBM00E) 

PURGE <SYMSUB/NA) 

DEFINE (TAPE10sSYMSU8/CT»StM*R) 
RFL(llOOOO) 

LDSET (PRESET* INDEFtMAPsSB) 
LDSET(FILES*TAPE4/TAPE10) 

LDSET (FILES*TAPE6/TAPE7/TAPE8/TAPE9) 
LGO. 

GOTO (CATALOG) 

EXIT. 

CATALOG(TAPEIO*R) 


5.  Description  of  Control  Cards 


PICRUST(T277)  RANLET 

Job  card  or  job  statement. 

USER (RANLETD , P ASWORD) 

CHARGE (60220 IR , PROJECT) 

SETBSL(777) 

System  dependent  accounting  information. 

ATTACH (OLDPL=PICRUST) 

Assign  permanent  file  PICRUST  to  current  job  with  local  file  name  OLDPL. 
File  PICRUST  contains  the  UPDATE  program  library  for  the  PICRUST  code. 
UPDATE(F,L=A124) 

RETURN (OLDPL) 

Process  UPDATE  directives  supplied  in  the  job  stream  (on  file  INPUT) 
and  prepare  a  source  file  suitable  for  compilation. 

FTN ( I=C0MP ILE,L=0,0PT=1,PL=25000,A) 

RETURN (COMPILE) 

Produce  compile  code  or  object  module  for  the  PICRUST  code. 

FILE (TAPE4 , BT=C , RT=S) 

FILE (TAPE6 ,BT=C ,RT=S) 

FILE (TAPE7 , BT-C , RT=S ) 

FILE (TAPES . BT»C , RT=S) 

FILE (TAPES , BT=C . RT=S ) 

FILE (TAPE 10 , BT=C , RT=S) 

Specify  the  record  type  of  each  of  the  indicated  local  files.  This 
assures  that  appropriate  system  routines  are  loaded  for  processing. 

ATTACH (TAPE4-SUBM0DE) 

Assign  permanent  file  SUBMODE  to  current  job  with  local  file  name  TAPE4. 
File  SUBMODE  is  the  substructure  mode  file. 

PURGE (SYMSUB/NA) 

DEFINE (TAPE10=SYMSUB/CT-S ,M=R) 

Assign  space  on  permanent  file  device  for  the  substructure  file  (TAPEIO) 
under  the  permanent  file  name  SYMSUB.  Information  written  on  TAPEIO 
during  the  execution  of  PICRUST  is  written  directly  on  the  permanent  file. 
RFL( 110000) 

Request  field  length  for  execution.  This  is  optional  on  many  systems. 
LDSET (PRESET=INDEF ,MAP»SB) 

LDSET (FILES=TAPE4 /TAPE 10) 

LDSET (FILES=TAPE6 /TAPE7 /TAPES /TAPE9 ) 

LGO. 

Select  load  options,  load, and  execute  the  compiled  version  of  PICRUST. 

The  required  set  of  input  data  is  supplied  in  the  job  stream  (on  file 
INPUT) .  The  FILES  parameter  must  be  supplied  on  a  LDSET  control  statement 
for  each  local  file  name  referenced  on  a  FILE  control  statement. 

GOTO (CATALOG) 

Unconditional  transfer  to  control  statement  CATALOG  if  the  execution 
terminates  normally. 

EXIT. 

Terminate  execution  or  provide  transfer  of  control  when  an  execution 
error  occurs. 

CATALOG (TAPE 10, R) 

List  information  about  each  logical  record  of  file  TAPEIO  on  file  OUTPUT. 
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6.  UPDATE  Directives  and  Input  Data  for  Sample  Problem  I 


The  UPDATE  directives  employed  to  prepare  the  substructure  file  for 
Sample  Problem  1  are  listed  below. 

*IDENT, STORAGE 

*DELETE , P ICRUST . 1 A , P ICRUST . 1 5 
COMMON  A (10000) 

NWORDA  =  10000 

The  size  of  the  dynamic  array  A  shown  above  is  in  excess  of  the  2221  words 
actually  required  for  Sample  Problem  1.  In  general,  a  dynamic  array  dimen¬ 
sioned  to  50000  words,  corresponding  to  a  central  memory  requirement  of 
about  210000  words  (octal),  is  sufficient  for  most  problems.  When  the 
dynamic  array  is  dimensioned  to  10000  words,  as  above,  a  field  length  of 
almost  110000  words  (octal)  is  required  for  an  execution  of  the  PICRUST  code 
The  set  of  input  data  shown  below  directs  the  PICRUST  code  to  prepare 
the  substructure  file  for  Sample  Problem  1. 
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.  UPDATE  Directives  and  Input  Data  for  Sample  Problem  2 


The  UPDATE  directives  employed  to  prepare  the  substructure  file  for 
Sample  Problem  2  are  listed  below. 


*IDENT, STORAGE 

*DELETE,PICRUST. 14.PICRUST. 15 
COMMON  A( 10000) 

NWORDA  =  10000 


The  size  of  the  dynamic  array  A  shown  above  is  somewhat  in  excess  of  the  size 
actually  required  for  Sample  Problem  2.  A  field  length  of  almost  110000  words 
(octal)  is  required  for  an  execution  of  PICRUST  when  the  dynamic  array  is 
dimensioned  as  above. 


The  set  of  input  data  shown  below  directs  the  PICRUST  code  to  prepare 


the  substructure  file  for  Sample  Problem  2. 
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APPENDIX  F  -  THE  USE  OF  USLOB 


1,  General  Information 

The  USLOB  computer  code  is  included  in  the  ELSHOK  suite  of  computer 
programs  to  solve  the  governing  equations  of  the  structure-fluid  interaction 
problem  outlined  in  Section  II.  The  USLOB  time-integration  processor  employs 
the  Runge-Kutta  technique,  in  the  modification  due  to  Gill  (Ref. [11]  ),  to 
perform  the  integration  in  time.  As  discussed  in  Section  III,  the  major 
portion  of  the  input  data  required  for  the  structure-fluid  interaction 
calculation  must  be  supplied  in  a  file  termed  the  "input  file".  In 
addition  to  this  file,  the  user  must  also  supply  time  step  information 
and  a  description  of  the  shock  wave  loading.  Printed  velocity-time 
histories  may  also  be  obtained. 

2.  Completion  of  Input  File 

Before  the  USLOB  code  may  be  applied  to  a  submerged  shock  response 
problem,  the  input  file  of  Section  III  must  be  completed.  For  problems 
without  internal  appendages,  the  shell-fluid  file  (Appendix  C)  serves 
as  the  input  file.  For  problems  involving  substructures,  all  sub- 

*) 

structure  files  (Appendix  E)  must  be  appended  to  the  shell-fluid  file. 

If  these  files  are  comprised  of  system -logical  records,  the  merging  of 
the  various  files  may  be  easily  accomplished  by  means  of  simple  control 
card  or  job  control  language  operations,  as  shown  below.  If  the  shell- 
fluid  and  substructure  files  are  not  comprised  of  system -logical 
records,  a  FORTRAN  computer  program  may  be  written  to  complete  the 
Input  file.  The  use  of  system-logical  records  is  the  recommended  procedure. 

*)  For  the  current  version  of  USLOB,  the  substructure  files  may  be  appended 
to  the  shell-fluid  file  in  any  convenient  order. 


The  following  set  of  control  cards  or  control  statements  may  be  used 


to  complete  the  Input  file  for  Sample  Problem  1: 


MERGE (T77)  RANLET 
USER (RANLETOtPASWORO) 

CHARGE (602201RtPROJECT) 

ATTACH (PART1«SHLFLU) 

ATTACH (PART2»SYMSUB) 

PURGE (PROBONE/NA} 

OEF I NE ( NEWF I LEaPROBONE/C  T^S • MaR ) 
COPYBR (PARTI •NEMFlLEtSBI 
COPYBR (PART2tNEWFlLEf 17) 
WRITEF(NEWFILEtl) 

CATALOG (NEWFlLEtR) 

EXIT. 


In  the  above  job  stream,  the  entire  shell-fluid  file  (SHLFLU)  and  the 
entire  substructure  file  (SYMSUB)  are  copied  to  the  input  file  (PROBONE) 
on  a  record-by-record  basis.  This  is  accomplished  by  using  the  system 
routine  COPYBR  to  copy  the  individual  binary  records  to  the  permanent 
file  PROBONE.  The  number  of  records  to  be  copied  may  be  obtained  by 
means  of  the  control  statement  CATALOG  or  may  be  computed  using  the 
expressions  given  in  Appendices  C  and  E.  An  end-of-flle  mark  (EOF)  is  pro 
vlded  on  the  input  file  for  the  structure-fluid  interaction  problem  by 
means  of  the  system  routine  WRITEF.  A  similar  set  of  control  cards  may 
be  used  to  complete  the  input  file  for  Sample  Problem  2. 

3.  Rules  of  Operation 

(1)  The  input  file  must  be  supplied  to  the  USLOB  code  as  a  local  file 
assigned  to  logical  unit  4  and  referred  to  as  TAPE4. 

(2)  The  generalized  velocity  file  (TAPElO  on  logical  unit  10)  must  be 
saved  if  plotted  velocity-time  histories  are  to  be  made  subsequently. 

(3)  If  the  problem  under  study  Involves  a  substructure  or  more  than  one 
substructure,  the  shell-substructure  transformation  file  (TAPES  on 
logical  unit  5)  must  be  saved  if  plots  of  substructure  responses  are 
to  be  made  subsequently. 
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.  Description  of  Input  Data 


The  Input  data  specifications  required  for  an  execution  of  the  USLOB 
code  are  described  below.  Following  Appendix  B, these  specifications  are 


written  in  a  style  similar  to  FORTRAN. 


Input  Set  1 


NTIME,  NSKIP,  NCHRG,  NQUAD,  NFINE,  KOUPLE,  NSUBS 


FORMAT  (715) 


where  NTIME  =  number  of  time  points  used  in  the  numerical  integration 
of  the  governing  equations.  The  corresponding  number  of 
time  steps  is  NTIME-1. 

NSKIP  =  Integer  parameter  selecting  the  time  interval  at  which 

the  generalized  velocities  are  saved  for  subsequent  pro¬ 
cessing.  The  time-history  information  is  saved  at  time 
points  starting  with  the  initial  time  point  (t  =  0)  and 
Incremented  by  NSKIP  thereafter.  For  example,  if 
NSKIP  =  2,  time-history  information  is  saved  at  time 
points  1,3, 5, etc. 

NCHRG  =  Integer  parameter  for  selecting  the  representation  of 
the  incident  spherical  shock  wave. 

NCHRG  =  1;  empirical  decaying  exponential  pressure-time 
history  of  Ref  [6]. 

NCHRG  =  2:  discretized  pressure-time  history  containing 
up  to  25  pressure-time  points. 

NQUAD  =  number  of  rays  dividing  a  quadrant  of  any  transverse 

cross  section  of  the  shell  S  into  a  coarse  (fixed)  grid 
used  to  locate  the  shock  front  and  to  evaluate  the  con¬ 
tributions  of  the  pressure-  and  velocity-time  histories 
away  from  the  front.  For  most  ELSHOK  applications,  use 
NQUAD  =  9. 

NFINE  *  number  of  rays  dividing  a  sector  defined  by  NQUAD  into 
a  fine  grid  used  to  evaluate  the  contributions  of  the 
pressure-  and  velocity-time  histories  in  the  vicinity  of 
the  shock  front.  The  grid  defined  by  NFINE  is  established 
at  the  shock  front  and  in  essence  moves  with  the  front. 

For  most  ELSHOK  applications,  use  NFINE  =  11. 

KOUPLE  =  integer  parameter  for  selecting  the  manner  of  introducing 
the  inertia  coupling,  produced  by  concentrated  masses, 
into  the  modal  equations  of  the  shell  S.  If  no  con¬ 
centrated  masses  are  attached  to  S  (NPTM  =  0  in  PIFLASH) , 
supply  any  valid  integer  to  satisfy  the  input  format. 
KOUPLE  *  1  ;  full  (exact)  coupling  of  the  modes  of  S. 
KOUPLE  =•  2  :  partial  (approximate)  coupling  of  the  modes 
of  S. 

NSUBS  *  number  of  substructures  included  in  the  analysis. 
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Note  1  : 


Note  2 


Note  3 


To  avoid  loss  of  the  generalized  velocities  at  the  final  time 
point,  the  parameter  NTIME  must  be  selected  in  conjunction  with 
the  parameter  NSKIP  such  that  the  calculation  end£  on  a  time 
point  for  which  the  generalized  velocities  are  saved. 

:  When  KOUPLE  =  2,  simplifications  based  upon  engineering  judgment 
are  employed  in  the  representation  of  the  inertia  coupling,  due 
to  concentrated  masses,  of  the  modes  of  the  shell  S,  leading  to 
reduced  central  memory  requirement  and  computational  time  for  a 
given  problem.  Specifically,  when  KOUPLE  =  2,  full  inertia 
coupling  between  the  N=0  torsional  and  rolling  modes  and  the  N=1 
rigid  body  and  whipping  modes  of  the  shell  S  is  included  in  the 
analysis.  However,  the  inertia  coupling  between  the  N=0  breathing 
modes  and  the  modes  having  N^2  with  the  above  sets  of  modes  is 
neglected.  Inertia  coupling  is  included  between  the  modes  in  each 
of  the  groups  of  modes  N=0  (breathing),  2,  3,  ...,  but  no  inertia 
coupling  between  these  individual  groups  is  considered.  This 
partial  coupling  has  been  found  to  yield  acceptable  results  when 
used  in  conjunction  with  the  simplified  modeling  technique  of 
Appendix  A  (full  and  compartment  models  of  S)  for  problems  in 
which  localized  slde-on  loading  acts  upon  the  portion  of  S  cor¬ 
responding  to  the  compartment  model. 

:  Since  USLOB  utilizes  main  memory  management,  it  is  quite  easy  to 
reduce  the  core  requirement  or  central  memory  to  the  minimum 
size  required  for  a  given  execution.  This  may  be  accomplished  by 
specifying  a  zero  or  negative  value  of  NTIME  in  an  otherwise  complete 
set  of  input  data.  An  execution  of  USLOB  with  this  set  of  input 
data  will  then  supply  information  from  which  the  desired  size  of 
the  dynamic  array  may  be  determined.  When  determining  the  minimum 
allowable  memory  size  for  a  given  problem,  the  dimension  of  the 
dynamic  array  in  blank  common  and  the  variable  NWORDA  which  supplies 
this  dimension  to  the  various  subprograms  in  the  USLOB  code  must  be 
set  to  at  least 

NWORDA  =  NUMBER  *  (12  +  7  *  NTORSN)  +.  15*NSUBS  +  6 

where  the  parameters  NUMBER  and  NTORSN  are  identical  to  the  input 
items  of  the  same  names  defined  in  Input  Set  1  of  Appendix  C. 
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Input  Set  2 

DELT,  XLOAD,  RLOAD,  SURCUT 


FORMAT  (4E10.0) 


where  DELT  =  size  of  time  step. 

XLOAD  =  longitudinal  coordinate,  referred  to  the  X-axis  of  the 

shell  S  (Fig.  1) ,  of  the  point  from  which  the  shock  wave 
emanates. 

RLOAD  =  radial  polar  coordinate,  measured  from  the  X-axis  of  the  shell 
S,  of  the  point  from  which  the  shock  wave  emanates. 

SURCUT  =  surface  cutoff  time. 

SURCUT  <0.0  :  no  surface  cutoff  effects  taken  into  account. 
SURCUT  >0.0  :  the  incident  wave  is  set  to  zero  at  the 
surface  cutoff  time  specified.  This  time  should  correspond 
to  the  time  at  which  the  surface  cutoff  wave  first  reaches 
the  submerged  structure. 

Note:  The  user  should  refer  to  Appendix  C  for  a  discussion  of  the 

location  of  the  origin  of  the  circumferential  coordinate  angle  6. 


IF  (NCHRG  .EQ.  2)  GO  TO  60 
WCHRG,  PZMLT,  PZEXP,  THMLT,  THEXP 
GO  TO  65 


FORMAT  (5E10.0) 


where 


Note: 


WCHRG  =  weight  of  spherical  charge  (W,  below). 

PZMLT  =  multiplicative  constant  (k^,  below)  for  Incident  pressure. 

PZEXP  =  spatial  decay  constant  below)  for  Incident  pressure. 

THMLT  =  multiplicative  constant  (k^.  below)  for  time  constant  of 
exponential  decay. 

THEXP  =  spatial  decay  constant  (k^,  below)  for  time  constant  of 
exponential  decay. 

The  above  Input  parameters  apply  to  the  empirical  pressure-time 
history  described  by  (Ref.  [6]) 

Pj(R,t)  =  kj  /R)^2  gjjp 

6^(R,t)  =  k^  W^^^  (W^^^  /R)  ^ 

where  p^  Is  the  Incident  pressure,  0^  Is  the  time  constant  of  expo¬ 
nential  decay,  R  Is  the  distance  from  the  origin  of  the  spherical 
wave  to  the  point  of  Interest,  t  represents  time,  and  t  denotes  the 
time  after  arrival  of  the  shock  wave  at  the  point  of  Interest.  The 
following  diagram  Illustrates  the  pressure-time  history  for  this  option 


TIME  AFTER  ARRIVAL 


Input  Set  3  for  Discretized  Pressure-Time  History  (NCHRG  =2) 


60  NSHAPE  FORMAT  (15) 

(TC(J),  PC(J),  J  =  1,  NSHAPE)  FORMAT  (8E10.0) 

DECAY  FORMAT  (ElO.O) 

where  NSHAPE  =  number  of  pressure-time  points  to  be  furnished  for 

description  of  incident  pressure. 

TC  =  array  containing  values  of  time  after  arrival  at 
which  incident  pressures  are  supplied. 

PC  =  array  containing  values  of  pressure  corresponding 
to  times  in  array  TC. 

DECAY  =  time  constant  of  exponential  decay  for  the  decaying 
exponential  tail  ending  the  incident  pressure-time 
history.  If  a  decaying  exponential  tail  is  not 
required,  use  DECAY  _<  0.0. 

Note  1 ;  A  typical  pressure-time  history  for  this  option  is  shown  below. 

NCHAR6  •  2 


In  the  above  diagram,  th.  symbols  P_  and  0^  correspond  to  the  input 
parameters  PC (NSHAPE)  and  DECAY,  respectively. 

Note  2;  For  the  current  version  of  USLOB,  the  input  parameter  NSHAPE  must 
not  exceed  25.  This  limitation  is  easily  removed  by  means  of 
the  UPDATE  batch  editing  system. 


FORMAT  (1015) 


where  LISTIT(l)  =  output  flag  for  geometry  of  shell  S. 

LISTIT(2)  =  output  flag  for  natural  frequencies  and  cir¬ 
cumferential  wave  numbers  of  modes  of  S. 

LISTIT(3)  =  output  flag  for  normal  component  of  modes  of  S. 

LISTIT(4)  =  output  flag  for  circumferential  wave  numbers  of 

surface  expansion  functions. 

LISTIT(5)  =  output  flag  for  surface  expansion  functions. 

LISTIT(6)  =  output  flag  for  fluid-shell  transformation  matrix 
of  Eq.  (5) . 

LISTIT(7)  =  output  flag  for  virtual  mass  array. 

LISTIT(8)  =  output  flag  for  Inverse  of  virtual  mass  array. 

LISTIT(9)  =  output  flag  for  coordinates  of  nodal  points  of 

every  substructure. 

LISTIT(IO)  =  output  flag  for  shell-substructure  transformation 
matrix  of  Eq.  (9)  of  every  substructure. 

LISTIT(ll)  =  output  flag  for  constraint  modes  of  Eq.  (10)  of 
every  substructure. 

LISTIT(12)  =  output  flag  for  equation  numbers  or  problem  unknowns 
associated  with  every  substructure. 

LISTIT(13)  =  output  flag  for  fixed-base  natural  frequencies  of 
every  substructure. 

LISTIT(14)  =  output  flag  for  fixed-base  modes  of  every 
substructure. 

LISTIT(15)  =  output  flag  for  portion  of  modal  system  mass  matrix. 

associated  with  the  shell  S  after  inclusion  of  effects 
of  shell  and  concentrated  masses. 

LISTIT(16)  “  output  flag  for  expansion  coefficients  of  constraint 
inodes  of  every  substructure. 

LISTIT(17)  =  output  flag  for  portion  of  modal  system  mass  matrix 
associated  with  coupling  between  the  shell  S  and 
every  substructure. 

LISTIT(18)  =  output  flag  for  interface  force  coefficients  associated 
with  connections  between  the  shell  S  and  every 
substructure . 

LISTIT(19)  *  output  flag  for  portion  of  modal  system  mass  matrix 

associated  with  the  shell  S  after  inclusion  of  effects 
of  shell,  concentrated  masses,  and  every 
substructure. 

LISTIT(20)  =  output  flag  for  generalized  velocities  of  the  shell 
S  and  every  substructure. 

Note;  If  an  output  flag  is  set  to  1,  output  is  supplied.  If  an  output 

flag  is  set  to  0,  no  output  is  supplied. 


NPTSHL  FORMAT  (15) 

(LBOSEG(J),  LBOSPT(J),  J  =  1,  NPTSHL)  FORMAT  (1215) 


where  NPTSHL  =  number  of  shell  stations  at  which  normal  velocity-time 
histories  are  to  be  printed.  If  no  shell  responses  are 
desired  from  USLOB,  use  NPTSHL  <  0  and  supply  no  further 
data  for  this  set  of  input. 

LBOSEG(J)  =  segment  of  the  shell  S  in  which  the  J-th  velocity 
station  is  located. 

LBOSPT(J)  =  mesh  point  in  segment  LBOSEG(J)  of  the  shell  S  at 

which  a  normal  velocity-time  history  is  to  be  printed. 


Note  1;  The  normal  velocity-time  histories  at  the  above  shell  stations 
are  computed  at  circumferential  coordinate  angles  6=0,  90, 
180,  270  deg. 

Note  2;  Physical  points  on  the  shell  S  are  located  with  reference  to 
a  full  mathematical  model  in  accordance  with  the  conventions 
established  for  the  B0S0R4  code. 


Input  Set  6 


DO  70  K  »  1,  NSUBS 

READ  NPTSUB  FORMAT  (15) 

70  READ  (NSTAT(J),  NQUAN(J),  J  =  1,  NPTSUB)  FORMAT  (1215) 

where  NPTSUB  =  number  of  velocity-time  histories  which  are  to  be  printed 
for  the  K-th  substructure.  If  no  responses  are  desired 
from  USLOB  for  the  K-th  substructure,  use  NPTSUB  <  0  and 
supply  no  further  data  for  the  K-th  substructure . 

NSTAT(J)  =  nodal  point  of  the  K-th  substructure  at  which  a  velocity¬ 
time  history  is  to  be  printed. 

NQUAN(J)  =  Identification  number  of  physical  degree  of  freedom  at 

nodal  point  NSTAT(J)  for  which  a  velocity-time  history  is 
to  be  printed. 

NQUAN(J)  =  1;  local  x-component  of  translation. 

NQUAN(J)  *  2:  local  y-component  of  translation. 

NQUAN(J)  =  3:  local  z-component  of  translation. 

NQUAN(J)  =  4:  local  x-component  of  rotation. 

NQUAN(J)  =5;  local  y-component  of  rotation. 

NQUAN(J)  =  6:  local  z-component  of  rotation. 


Note;  The  substructures  are  processed  sequentially  in  the  order  that  the 
substructure  files  were  appended  to  the  shell-fluid  file. 


USL0B(T1777)  RANLET 
USER (RANLETOvPASMORO) 

CHARGE <602201R«PROUECT> 

SETBSL<7777) 

ATTACH ( OLOPL»USLOB ) 

UPDATE (F«L>A124) 

RETURN (OLDPL) 

FTN(IsCOHPILE«L«0tOPT»2fPL«30000tA> 
RETURN (COMPILE) 
FILE(TAPE4#BT*C#RT*S) 

ATTACH (TAPE4-PR0B0NE) 

PURGE (PROB105/NA) 

DEFINE (TAPE5*PROB105/CT*S»M*R) 

PURGE (PROB 11 0/NA) 

DEFINE ( TAPE 1 OapROB 1 1 0/CT«S • M«R ) 
RFL(145000) 

LDSET (PRESET*INOEF,HAP*SB) 
LDSET<FILES*TAPE4) 

LGO. 

GOTO (CATALOG) 

EXIT. 

CATALOG(TAPE5fR) 

CATALOG(TAPElOfR) 
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6.  Description  of  Control  Cards 


USLOB(T1777)  RANLET 

Job  card  or  job  statement. 

USER (RANLETD , PASWORD) 

CHARGE (60220 IR , PROJECT) 

SETBSL(7777) 

System  dependent  accounting  information. 

ATTACH (0LDPL=USL0B) 

Assign  permanent  file  USLOB  to  current  job  with  local  file  name  OLDPL. 

File  USLOB  contains  the  UPDATE  program  library  for  the  USLOB  code. 

UPDATE (F,L=A124) 

RETURN (OLDPL) 

Process  UPDATE  directives  supplied  in  the  job  stream  (on  file  INPUT) 
and  prepare  a  source  file  suitable  for  compilation. 

FTN ( I-COMPILE , L=0 , 0PT=2 , PL=30000 , A) 

RETURN(COMPILE) 

Produce  compiled  code  or  object  module  for  the  USLOB  code. 

FILE (TAPE4 , BT=C , RT=S ) 

Specify  the  record  type  of  the  Indicated  local  file.  This  assures 
that  appropriate  system  routines  are  loaded  for  processing. 

ATTACH (TAPE4-PR0B0NE) 

Assign  permanent  file  PROBONE  to  current  job  with  local  file  name  TAPE4. 

File  PROBONE  is  the  input  file  for  Sample  Problem  1. 

PURGE (PROB105/NA) 

DEFINE (TAPE5-PROB 1 05 /CT»S , M»R) 

PURGE (PROBllO/NA) 

DEFINE (TAPEIO-PROBI 10/CT-S ,M»R) 

Assign  space  on  permanent  file  device  for  the  shell-substructure 
transformation  file  (TAPES)  and  the  generalized  velocity  file  (TAPEIO) 
under  the  permanent  file  names  PROB105  and  PROBllO,  respectively. 

Information  written  on  TAPES  and  TAPEIO  during  the  execution  of  USLOB 
is  written  directly  on  the  permanent  files. 
ilFL(  145000) 

Request  field  length  for  execution.  This  is  optional  on  many  systems. 

LDSET (PRESET-INDEF ,MAP=SB) 

LDSET(FILES-TAPE4) 

LGO. 

Select  load  options,  load,  and  execute  the  compiled  version  of  USLOB. 

The  required  set  of  input  data  is  supplied  in  the  job  stream  (on  file  INPUT), 
The  FILES  parameter  must  be  supplied  on  a  LDSET  control  statement  for 
each  local  file  name  referenced  on  a  FILE  statement. 

GOTO (CATALOG) 

Unconditional  transfer  to  first  occurrence  of  control  statement  CATALOG 
if  the  execution  terminates  normally. 

EXIT. 

Terminate  execution  or  provide  transfer  of  control  when  an  execution 
error  occurs. 

CATALOG (TAPES, R) 

CATALOG (TAPE 10, R) 

List  Information  about  each  logical  record  of  files  TAPES  and  TAPEIO  on 
file  OUTPUT. 


*IDENT, STORAGE 
*DELETE , USLOB . 24 . USLOB . 25 
COMMON  A (25000) 
NWORDA  *  25000 


The  size  of  the  dynamic  array  A  shown  above  is  somewhat  in  excess  of  the 


size  determined  by  a  restricted  execution  of  the  USLOB  code.  In  this 


restricted  execution,  the  dimension  of  the  array  A  and  the  variable  NWORDA 


were  set  to  75  words  and  the  input  parameter  NTIME  was  set  to  0  (zero) . 


All  other  input  parameters  were  identical  to  those  listed  below.  When  the 


dynamic  array  is  dimensioned  as  above,  a  field  length  of  almost  145000 


words  (octal)  is  required  for  an  execution  of  the  USLOB  code. 


The  set  of  input  data  shown  below  directs  the  USLOB  code  to  perform 


the  shock  response  calculation  for  Sample  Problem  1. 


191  2  2  9  11  1  1  NTINEtNSKiPtNCHRetNOUAOtNFINEtKOUPLEtNSUBS 

2.50E-0S  127.S12S  aSS.6129  0.0  DELTtXLOADtRLOAOtSURCUT 

2  NSMAPE 

0.0  10.0  1.0  10.0  (TCUltPCCJUJaltNSHAPE) 

0.0  DECAY 

0000000000  OUTPUT  FLASS  1-10 

0000000000  OUTPUT  FLASS  11-20 

S  NPTSHL 

2  1  2  12  2  22  2  92  2  53  LBOSEStLBOSPT 

A  NPTSUB 

1  3  *  3  19  3  29  3  NSTATtNOUAN 
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8.  UPDATE  Directives  and  Input  Data  for  Sample  Problem  2 


The  UPDATE  directives  employed  to  compute  the  submerged  shock  response 
for  Sample  Problem  2  are  listed  below. 

*IDENT, STORAGE 
*DELETE , USLOB . 24 , USLOB . 2  5 
COMMON  A( 15000) 

NWORDA  =  15000 

The  size  of  the  dynamic  array  A  shown  above  is  somewhat  in  excess  of  the 
size  determined  by  a  restricted  execution  of  the  USLOB  code.  In  this 
restricted  execution,  the  dimension  of  the  array  A  and  the  variable  NWORDA 
were  set  to  75  words  and  the  input  parameter  NTIME  was  set  to  0  (zero). 

All  other  input  parameters  were  identical  to  those  listed  below.  When  the 
dynamic  array  is  dimensioned  as  above,  a  field  length  of  almost  125000 
words  (octal)  is  required  for  an  execution  of  the  USLOB  code. 

The  set  of  input  data  shown  below  directs  the  USLOB  code  to  perform 
the  shock  response  calculation  for  Sample  Problem  2. 


191  2  1  9 

2.50E>0S  127.8125 

3.8  3.S80375*5 

8  8  8  8 

8  8  8  8 

4 

1  1  2  19 

4 

13  8  3 


11  2  1  NTINE»NSKI9tNCN8etNOUAO*NFINC*KOUPLC*N$UBS 

58.4375  8.9  DCLT.XLOAO.KLOADf SIMCUT 

1.13  3.357439-5  -0.22  «CHR8»PZMLTtPZEXPtTNHL7tTNCXP 


8 

8 


8 

8 


8 

9 

19 


8  8 

8  9 

NPTSHL 

NP75U8 


OUTPUT  FLA8S  1-18 
OUTPUT  FLA8S  11-20 

LBOSEOtLBOSPT 


15  3  25  3 


NSTATttMUAN 
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APPENDIX  G  -  THE  USE  OF  PUSLOB 


1.  General  Information 

The  PUSLOB  computer  code  is  included  in  the  ELSHOK  suite  of  com¬ 
puter  codes  to  produce  plots  of  velocity-time  histories  for  any  point 
on  the  main  body  S  and  for  any  nodal  point  on  any  substructure  a.  As 
discussed  in  Section  III,  velocity  plots  for  the  shell  may  be  provided 
in  the  local  coordinate  system  and/or  in  the  global  coordinate  system. 

In  the  local  system,  the  inward  normal  (n) ,  meridional  (s),  and  circum¬ 
ferential  (9)  directions  of  Fig.  1  are  employed.  In  the  global  system, 
the  athwartship  (-Z) ,  fore-aft  (+X) ,  and  downward  (-Y)  directions  of 
Fig.  1  are  used.  All  plotting  for  a  substructure  is  performed  with 
reference  to  the  coordinate  system  of  the  given  substructure.  In  addition 
to  plotted  results,  coded  card-image  time  histories  may  also  be  prepared. 

All  plotting  of  velocity-time  histories  using  the  PUSLOB  code  is 
performed  in  two  stages.  The  first  stage  is  the  execution  of  PUSLOB 
to  produce  a  file,  called  the  "plot  file",  which  must  be  saved  on  a 
permanent  file  device  or  any  other  storage  unit  accessible  by  means 
of  a  TEKTRONIX  graphics  terminal.  The  second  stage  of  the  operation 
is  the  transfer  of  the  plot  commands  and  related  data  from  the  saved 
plot  file  to  a  TEKTRONIX  graphics  terminal  (usually  via  some  type  of 
COPY  or  LIST  command)  and  the  subsequent  production  of  the  desired  time 
histories  on  a  TEKTRONIX  hard  copy  unit. 


(1)  An  input  file,  a  generalized  velocity  file,  and  a  shell-substructure 
transformation  file  (for  those  problems  involving  at  least  one  sub¬ 
structure)  must  be  available  prior  to  the  execution  of  the  PUSLOB  code. 

(2)  The  input  file  must  be  supplied  to  PUSLOB  as  a  local  file  assigned 
to  logical  unit  4  and  referred  to  as  TAPE4. 

(3)  The  generalized  velocity  file  must  be  supplied  to  PUSLOB  as  the 
local  file  TAPE 10. 

(4)  If  a  shell-substructure  transformation  file  is  required  for  the  problem 
under  study,  it  must  be  supplied  to  PUSLOB  as  the  local  file  TAPES. 

(5)  The  user  must  supply  control  cards  or  control  statements  to  ensure 
that  the  plot  file  produced  by  PUSLOB  (TAPE99  on  logical  unit  99) 
is  comprised  of  system-logical  records. 

(6)  The  plot  file  (TAPE99)  must  be  saved  for  subsequent  use  with  TEKTRONIX 
graphics  equipment. 

(7)  The  user  must  decide  the  disposition  of  any  coded  card-image  time 
histories  produced  by  PUSLOB.  These  time  histories  are  written  on 
a  coded  file  assigned  to  logical  unit  7  and  referred  to  as  TAPE7. 


.  Description  of  Input  Data 


The  input  data  specifications  required  for  an  execution*  of  the  PUSLOB 
code  are  described  below.  Following  Appendix  B,  these  specifications  are 
written  in  a  style  similar  to  FORTRAN. 


Input  Set  1  FORMAT  (515) 

NTIME,  NSKIP,  NSUBS,  NTEK,  NCARD 

where  NTIME  =  number  of  time  points  to  be  used  in  preparing  velocity¬ 
time  histories.  The  value  supplied  may  not  exceed  the 
number  used  In  the  execution  of  the  USLOB  code  which 
produced  the  generalized  velocity  file  for  the  problem 
under  study. 

NSKIP  =  Integer  parameter  for  processing  generalized  velocities. 

For  the  current  version  of  PUSLOB,  supply  the  same  value 
of  NSKIP  used  in  the  execution  of  the  USLOB  code  which 
produced  the  generalized  velocity  file  for  the  problem 
under  study. 

NSUBS  =  number  of  substructures.  This  parameter  may  be  set  to 
zero  if  no  velocity-time  histories  are  desired  from 
PUSLOB  for  a  problem  involving  substructures. 

NTEK  »  integer  parameter  Indicating  whether  or  not  plotted 

velocity-time  histories  are  to  be  produced  (l=yes;  0=no) . 
If  NTEK>=0,  the  plot  file  (TAPE99)  is  not  produced. 

NCARD  =  Integer  parameter  indicating  whether  or  not  coded  card- 
image  time  histories  are  to  be  produced  (l=yes;  0=no) . 

If  NCARD=0,  the  file  used  for  the  coded  velocity-time 
histories  (TAPE?)  is  not  produced. 

Note  :  The  current  version  of  PUSLOB  employs  main  memory  management , 
but  does  not  allow  a  restricted  execution  to  determine  the 
minimum  size  of  the  dynamic  array  needed  for  a  given  problem. 
Thus,  the  user  must  not  supply  a  zero  or  negative  value  for 
the  input  parameter  NTIME  when  executing  the  PUSLOB  code. 
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FORMAT  (3E10.0) 


Input  Set  2 

DELT,  XMULT,  YMULT 

where  DELT  *  size  of  time  step  used  in  execution  of  USLOB  code  which 
produced  the  generalized  velocity  file  for  the  problem 
under  study. 

XMULT  =  multiplicative  scale  factor  for  time  or  horizontal 
axis  of  every  plotted  velocity-time  history.  This 
parameter  allows  the  units  of  the  time  scale  to  be  changed 
easily. 

YMULT  =  multiplicative  scale  factor  for  velocity  or  vertical 
axis  of  every  plotted  velocity-time  history.  This 
parameter  allows  the  units  of  the  velocity  scale  to  be 
changed  easily. 

NOTE;  The  scale  factors  do  not  apply  to  the  coded  card-image  time  histories 
written  on  local  file  TAPE?.  Velocity-time  histories  prepared 
for  coded  file  TAPE?  have  units  consistent  with  those  employed 
tor  the  execution  of  program  USLOB. 

Input  Set  3  FORMAT  (3A10) 

(LTITLE(J),  J  =  1,3) 

where  LTITLE  =  thirty-character  alphameric  main  title  or  heading 

to  be  placed  on  every  velocity-time  history.  During 
the  execution  of  PUSLOB,  a  ten-character  alphameric 
identification  label,  supplied  by  the  user  for  each  time 
history,  is  appended  to  the  main  heading  contained  in 
array  LTITLE. 

Input  Set  4  FORMAT  (4A10) 

(NOTEX(J),  J  =  1,2),  (NOTEY(J),  J  =  1,2) 

where  NOTEX  =  twenty-character  alphameric  label  for  time  or  horizontal 
axis  of  every  plotted  velocity-time  history. 

NOTEY  =  twenty-character  alphameric  label  for  velocity  or  vertical 
axis  of  every  plotted  velocity-time  history. 


INPUT  Set  5 


DO  90  K  -  1,3 
READ  NPTSHL 

IF  (NPTSHL  .LE.  0)  GO  TO  90 
DO  80  J  =  1,  NPTSHL 
80  READ  NGAGE,  LBOSEG,  LBOSPT, 

NSPHC,  NANG,  ANGDEG 

90  CONTINUE 

where  NPTSHL  =  number  of  shell  stations  at  which  velocity-time  histories 
are  to  be  prepared  using  K-th  component  of  modes  of  S. 

If  no  shell  responses  are  required  for  a  given  component 
of  velocity,  use  NPTSHL  ^  0  and  supply  no  further  data  for 
that  component. 

K  =  1  ;  Inward  normal  velocity  (w) . 

K  =  2  :  meridional  velocity  (u) . 

K  =  3  ;  circumferential  velocity  (v) . 

NGAGE  =  ten-character  alphameric  label,  to  be  appended  to  main 

title,  for  identifying  shell  station  to  which  a  velocity¬ 
time  history  pertains. 

LBOSEG  =  segment  of  shell  S  in  which  the  J-th  velocity  station,  of 
the  NPTSHL  stations  considered,  is  located. 

LBOSPT  =  mesh  point  in  segment  LBOSEG  of  shell  S  at  which  J-th 
velocity  station  is  located. 

NSPHC  =  dual  purpose  integer  flag  indicating  (a)  whether  or  not 
current  K-th  local  component  of  velocity  is  to  be  used 
for  projection  into  global  coordinate  system  of  shell  S 
and  (b)  whether  or  not  a  hard  copy  (i.e.,  plot  or  card- 
image  history)  of  resulting  K-th  projected  global  component 
of  velocity  is  to  be  produced.  When  producing  global 
components  of  velocity  by  projection,  the  hard-copy 
options  indicated  by  NSPHC  for  K  =  1,  2,  and  3  apply 
to  the  athwartshlp  ,  fore-aft,  and  downward  global  com¬ 
ponents  of  the  velocity  of  S,  respectively,  at  the  station 
of  Interest. 

NSPHC  =  0  :  do  not  use  current  K-th  local  component  of 
velocity  for  projection  into  global  coordinate  system  of  S. 
NSPHC  =  1  :  use  current  K-th  local  component  of  velocity 
for  projection  into  global  coordinate  system  of  S,  and 
provide  hard  copy  of  resulting  K-th  global  component  of 
velocity. 

NSPHC  =-l  :  use  current  K-th  local  component  of  velocity 
for  projection  into  global  coordinate  system  of  S,  but 
do  not  provide  hard  copy  of  resulting  K-th  global  component 
of  velocity. 

NANG  “  number  of  angular  locations,  at  current  J-th  shell  station, 
at  which  velocity-time  histories  are  to  be  produced.  For 
the  current  version  of  PUSLOB,  use  NANG  =  1  or  2. 

ANGDEG  ■  circumferential  coordinate  angles  (in  degrees)  of  above 
NANG  angular  locations. 


FORMAT  (15) 

FORMAT  (A10,4I5,2E10.0) 
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Note  1 ;  When  global  components  of  the  shell  velocity  are  being  produced  at 
a  given  shell  station  (NSPHC  0) ,  the  labels  supplied  as  NGAGE 
for  K  =  1,2,  and  3  apply  to  the  athwartship,  fore-aft,  and  down¬ 
ward  global  components  of  the  velocity  of  S,  respectively. 

Note  2;  The  parameter  NSPHC  overrides  the  NTEK  and  NCARD  parameters  on  a 
component-by-component  basis  if  and  only  if  global  components  of 
the  shell  velocity  are  being  produced  by  projection  of  the  local 
components . 

Note  3;  When  producing  global  components  of  shell  velocity  at  a  given 

station,  consistent  input  data  must  be  supplied  for  that  station 
for  each  of  the  local  components  (K  =  1,2,3)  of  shell  velocity. 

When  more  than  one  shell  station  is  being  considered,  the  input 
parameters  for  those  stations  must  be  arranged  in  the  same  sequence 
for  each  local  component  (K  =  1,2,3).  Input  for  which  NSPHC  =  0  may 
be  interspersed  with  input  for  which  NSPHC  =  +  1,  if  desired,  without 
upsetting  the  sequence  of  data  used  for  projection. 

Input  Set  6 

DO  110  K  =  1,  NSUBS 
READ  NPTSUB 

IF  (NPTSUB  .LE.  0)  GO  TO  110 
DO  100  J  =  1,  NPTSUB 

100  READ  NGAGE,  NSTAT,  NDOFPT,  NQUAN 

110  CONTINUE 

where  NPTSUB  =  number  of  stations  at  which  velocity-time  histories  are  to 
be  prepared  for  the  K-th  substructure.  If  no  responses  are 
required  for  the  K-th  substructure,  use  NPTSUB  <  0  and  supply 
no  further  data  for  the  K-th  substructure. 

NGAGE  =  ten-character  alphameric  label,  to  be  appended  to  main  title, 
for  identifying  nodal  point  of  K-th  substructure  to  which  a 
velocity-time  history  pertains. 

NSTAT  “  nodal  point  of  the  K-th  substructure  at  which  a  velocity¬ 
time  history  is  to  be  printed. 

NDOFPT  *  number  of  physical  degrees  of  freedom,  at  current  J-th 
substructure  station,  for  which  velocity-time  histories 
are  to  be  produced.  For  the  current  version  of  PUSLOB, 
use  NDOFPT  =  1,2,  or  3 . 

NQUAN  ■  Identification  numbers  of  above  NDOFPT  physical  degrees 
of  freedom. 

NQUAN  *  1;  local  x-component  of  translation. 

NQUAN  =  2:  local  y-coraponent  of  translation. 

NQUAN  =  3:  local  z-component  of  translation. 

NQUAN  =  A;  local  x-component  of  rotation. 

NQUAN  *  5:  local  y-component  of  rotation. 

NQUAN  ■  6:  local  z-component  of  rotation. 

Note;  The  substructures  are  processed  sequentially  in  the  order  that  the 

substructure  files  were  appended  to  the  shell-fluid  file. 


FORMAT  (15) 
FORMAT  (A10,5I5) 
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4.  Sample  Control  Cards 

A  sample  set  of  control  cards  or  job  control  language  statements  for 
an  execution  of  PUSLOB  under  a  NOS  1  operating  system  is  shown  below.  The  next 
subsection  contains  a  brief  description  of  these  control  statements.  Various 
other  arrangements  of  the  job  steps  are  possible. 


PUSLOB (T1 77)  RANLET 
USER (RANLETOf PASWORO) 

CHARGE (602201R«PROJECT) 

SETBSL(777) 

ATTACH (OLOPL*PUSLOB) 
UP0ATE<F«LmA124) 

RETURN (OLOPL) 

FTN(IsCOMPILE«L«0*OPT«l*A) 

RETURN(COMPILE) 

FILE(TAPE4.BT*C*RT*S) 

FILE(TAPE99#BT«C*RT*S) 

ATTACH (TAPE4mPR0B0NE) 

ATTACH (TAPE5*PROB105) 
ATTACH(TAPE10«PROBnO) 

ATTACH <QTEK) 

PURGE (PLTFILE/NA) 

DEFINE (TAPE99sPLTFILE/CT«S.M«R) 
RFL(135000) 

LOSET(PRESET*INOEF,MAP=SBfLlBsQTEK) 
LOSET (FILESSTAPE4/TAPE99) 


5.  Description  of  Control  Cards 


PUSLOB(T177)  RANLET 

Job  card  or  job  statement. 

USER (RANLETD , PASWORD) 

CHARGE (60220 IR, PROJECT) 

SETBSL(777) 

System  dependent  accounting  Information. 

ATTACH (0LDPL=PUSL0B) 

Assign  permanent  file  PUSLOB  to  current  job  with  local  file  name  OLDPL. 
File  PUSLOB  contains  the  UPDATE  program  library  for  the  PUSLOB  code. 
UPDATE(F,L»A124) 

RETURN (OLDPL) 

Process  UPDATE  directives  supplied  in  the  job  stream  (on  file  INPUT) 
and  prepare  a  source  file  suitable  for  compilation. 

FTN ( I=C0MP ILE , L=0 . OPT- 1 , A) 

RETURN (COMPILE) 

Produce  compiled  code  or  object  module  for  the  PUSLOB  code. 

FII E (TAPE4 , BT=C , RT=S) 

FILE (TAPE99 , BT=C ,RT=S) 

Specify  the  record  type  of  each  of  the  indicated  local  files.  This 
assures  that  appropriate  system  routines  are  loaded  for  processing. 

ATTACH (TAPE4-PR0B0NE) 

ATTACH (TAPE5-PR0B 105) 

ATTACH(TAPE10=PROB1 10) 

Assign  permanent  files  PROBONE,  PROB105,  PROBllO  to  current  job  with  local 
file  names  TAPE4,  TAPES,  TAPElO,  respectively.  These  files  are  the  input 
file,  the  shell-substructure  transformation  file,  and  the  generalized 
velocity  file,  respectively. 

ATTACH (QTEK) 

Assign  permanent  file  QTEK  to  current  job  with  local  file  name  QTEK. 

File  QTEK  contains  subroutines  used  for  plotting  on  TEKTRONIX  graphics 
terminals  and  serves  as  a  library  file  for  satisfying  externals  required 
by  PUSLOB. 

PURGE (PLTFILE/NA) 

DEFINE(TAPE99-PLTFILE/CT-S,M-R) 

Assign  space  on  permanent  file  device  for  the  plot  file  (TAPE99)  under  the 
permanent  file  name  PLTFILE.  Information  written  on  TAPE99  during  the 
execution  of  PUSLOB  is  written  directly  on  the  permanent  file.  File 
TAPE99  must  be  listed  on  a  TEKTRONIX  graphics  terminal  to  produce  plots 
of  the  desired  velocity-time  histories. 

RFL( 135000) 

Request  field  length  for  execution.  This  is  optional  on  many  systems. 
LDSET (PRESET-INDEF ,MAP-SB ,LIB=QTEK) 

LDSET(FILES»TAPE4/TAPE99) 

LGO. 

Select  load  options,  load,  and  execute  the  compiled  version  of  PUSLOB. 

The  required  set  of  input  data  is  supplied  in  the  job  stream  (on  file 
INPUT).  The  FILES  parameter  must  be  supplied  on  a  LDSET  control  statement 
for  each  local  file  name  referenced  on  a  FILE  control  statement. 

EXIT. 


Terminate  execution. 


6.  UPDATE  Directives  and  Input  Data  for  Sample  Problem  1 


The  UPDATE  directives  employed  to  produce  plots  of  velocity-time 
histories  for  Sample  Problem  1  are  listed  below. 

*IDENT, STORAGE 
*DELETE .PUSLOB . 19 .PUSLOB . 20 
COMMON  A(25000) 

NWORDA  =  25000 


The  size  of  the  dynamic  array  shown  above  is  somewhat  in  excess  of  the  size 
actually  required  for  Sample  Problem  1.  In  general,  a  good  estimate  of  the 
size  of  the  dynamic  array  A  for  PUSLOB  is  provided  by  the  size  of  the  dynamic 
array  employed  for  the  execution  of  the  USLOB  code  for  the  problem  under 
study.  A  field  length  of  almost  135000  words  (octal)  is  required  for  an 
excutlon  of  PUSLOB  when  the  dynamic  array  is  dimensioned  as  above. 

The  set  of  input  data  shown  below  directs  the  PUSLOB  code  to  produce  the 
plot  file  from  which  velocity-time  plots  may  be  obtained  for  Sample  Problem  1. 


191  2  1  1  0  NTIMEfN$ICIP»NSUBS»NTEK»NCAIID 

2.50E-05  1000.0  1.0  OELTtXIHH.TtyNULT 

SAMPLE  PPOBLEM  1  (PROBOME}  (LTITLEUI tJ«lf 3) 

TIME  (MSECI  VELOCITY  (IN/SEC)  NOTEXtNOTEV 


5  NPTSHL-M 


V-ATH 

2/01 

2 

1 

1 

2 

0.0 

IBO 

0 

NBAOEtSB 

PT • NSPHC  t  NAN6 1 ANODEO 

V-ATH 

2/12 

2 

12 

1 

2 

0.0 

IBO 

0 

NBAOEtSO 

PT  t  NSBhC  t  N ANB  t  AN60EB 

V-ATM 

2/22 

2 

22 

1 

2 

0.0 

100 

0 

NOAOEtSe 

PT  *  NSPHC  »  N ANB*  ANBDEO 

V-ATM 

2/A2 

2 

42 

1 

2 

0.0 

100 

0 

NOAOEvSO 

PT .NSPHC.NANO. ANBOEO 

V-ATH 

2/53 

2 

93 

1 

2 

0.0 

IBO 

0 

NBAOEtSe 

PT  »  NSPHC • NANS t  ANBOEO 

9 

V-FMO 

2/01 

NPTSHL-U 
2  1 

^1 

2 

0.0 

100 

0 

NBABE tSe 

PT • NSPHC  t  NANS • ANSDEB 

V-PHO 

2/12 

2 

12 

•1 

2 

0.0 

100 

0 

NBABE »S6 

PT • NSPHC  *  NANS  *  ANSDEB 

V-FMO 

2/22 

2 

22 

•  1 

2 

0.0 

100 

0 

NOAOEtSe 

PT  t  NSPHC  «  NANO  «  ANBOEO 

V-FMO 

2/42 

2 

42 

•1 

2 

0.0 

100 

0 

NSABEtSB 

PT  t  NSPHC • NANS  t  ANSOEB 

v-fmo 

2/S3 

2 

93 

*1 

2 

0.0 

100 

0 

NOABEvSO 

PT  t  NSPHC  t HANS  «  ANSOEB 

5 

V-OMN 

2/01 

NPTSHL-V 
2  1 

•1 

2 

0.0 

100 

0 

NBABE f SO 

PT • NSPHC  *  NANB  *  ANSOEB 

V-OMN 

2/12 

2 

12 

•  1 

2 

0.0 

100 

0 

NBABE t SB 

PT  t  NSPHC •HANS • ANSOEB 

V-OMN 

2/22 

2 

22 

•  1 

2 

0.0 

100 

0 

NBABE vSB 

PT • NSPHC •NANS • ANSOEB 

V-OMN 

2/42 

2 

42 

•  1 

2 

0.0 

lOO 

0 

NBABEvSB 

PT • NSPHC • NANB • ANSOEB 

V-OMN 

2/S3 

2 

93 

•1 

2 

0.0 

100 

0 

NBABE tSB 

PT tNSPHC *NAH6* ANBOEO 

3 

TIP  I0*> 

NPTSUB 

*  1 

3 

NBABE t NST AT t NOOFPT • MOUAN 

JCN  (IS) 

19 

1 

3 

NBABE 

tNSTAT«ND0FPT« 

MOUAN 

MtO  (tS) 

29 

1 

3 

NBABE 

fNSTATtNOOFPT* 

MOUAN 
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.  UPDATE  Directives  and  Input  Data  for  Sample  Problem  2 


The  UPDATE  directives  employed  to  produce  plots  of  velocity-time 
histories  for  Sample  Problem  2  are  listed  below. 


*IDENT, STORAGE 
*DELETE , PUSLOB . 19 .PUSLOB . 20 
COMMON  A( 15000) 

NWORDA  =  15000 


The  size  of  the  dynamic  array  shown  above  is  somewhat  in  excess  of  the  size 
actually  required  for  Sample  Problem  2.  In  general,  a  good  estimate  of  the 
size  of  the  dynamic  array  A  for  PUSLOB  is  provided  by  the  size  of  the  dynamic 
array  employed  for  the  execution  of  the  USLOB  code  for  the  problem  under 
study.  A  field  length  of  almost  115000  words  (octal)  is  required  for  an 
execution  of  PUSLOB  when  the  dynamic  array  is  dimensioned  as  above. 

The  set  of  input  data  shown  below  directs  the  PUSLOB  code  to  produce  the 
plot  file  from  which  velocity-time  plots  may  be  obtained  for  Sample  Problem  2 


191  2  1  1  0  NTIHE*NSKIPtNSlieS*NTCK*NCARD 

2.50e«9S  1000.0  0.0033333  DCLTtXMULT.YNULT 

SAMPLE  PROBLEM  2  (PROBTMO)  (LTITLEC J) .J>1»3> 


2 

TIME 

(MSEC) 

NPTSML-R 

velocity  (PT/SEC) 

NOTEX.NOTEY 

V'ATH 

3/01 

3 

1 

1 

2 

0.0 

100.0  MB ABE  »SB • PT  *  NSPHC • N AN6 *  ANBOEB 

V-ATH 

2 

3/19 

3  19 

NPTSHL'U 

1 

2 

0.0 

100.0  NBABE t SB tPT.NSPNCtNANBt ANBOEB 

V-PRO 

3/01 

3 

1 

«1 

2 

0.0 

IBO.O  NBABEtSB.PT.NSPNC.NANBt ANBOEB 

V-PRO 

2 

3/19 

3  19 

mptshl-v 

•I 

2 

0.0 

V'DRN 

3/01 

3 

1 

•1 

2 

0.0 

IBO.O  NBABE*S6*PT«NSPHCtNAN6«ANB0EB 

V»0RN 

3 

3/19 

3  19 

nptsub 

•1 

2 

0.0 

IBO.O  NBABEvSB.PT.NSPHC.NANB.ANBOEB 

TIP  (061 

6 

1 

3 

NBA6E 

.NSTAT.NOOPPTfNOUAN 

JCN  (IS) 

IS 

1 

3 

nbabe.nstat. 

MID  (2SI 

2S 

1 

3 

NBABE 

.NSTAT.NDOPPTtNOUAN 

DISTRIBUTION  LIST 


DEPARTMENT  OF  DEFENSE 

Defense  Advanced  Rsch  Proj  Agency 
ATTN:  TIO 

Defense  Intelligence  Agency 
ATTN;  DB-4C2 
ATTN:  DB-4C2,  C.  Wiehle 
ATTN:  DT-IC 
ATTN:  DT-2 
ATTN:  RTS-2A 
AHN:  DB-4C3 

ATTN:  DB-4C,  Rsch,  Phys  Vuln  Br 
ATTN:  DB-4C1 

Defense  Nuclear  Agency 
ATTN:  SPSS 
4  cy  ATTN:  TITL 

Defense  Technical  Information  Center 

12  cy  ATTN:  DD 

Field  Command 

Defense  Nuclear  Agency,  Det  1 

Lawrence  Livermore  Lab 
ATTN:  FC-1 

Field  Command 

Defense  Nuclear  Agency 
AHN:  FCPR 
ATTN:  FCT 
ATTN:  FCTX 
ATTN:  FCTT,  G.  Ganong 
ATTN:  FCTT,  W.  Summa 
ATTN;  FCTXE 

Interservice  Nuclear  Weapons  School 
ATTN:  TTV 

Joint  Strat  Tgt  Planning  Staff 
ATTN:  NRl-STINFO  Library 
ATTN:  JLA,  Threat  Applications  D 
ATTN:  JLTW,  Rautenberg 
ATTN:  JLTW-2 
ATTN:  OOXT 
ATTN:  XPFS 

Under  Sec  of  Def  for  Rsch  &  Engrg 

ATTN:  Strategic  &  Space  Sys  (OS) 

DEPARTMENT  OF  THE  ARMY 

BMD  Advanced  Technology  Center 
ATTN:  ICRDABH-X 
ATTN:  ATC-T 

Chief  of  Engineers 

ATTN:  DAEN-RDL 
ATTN:  DAEN-MPE-T 

Dep  Ch  of  Staff  for  Ops  &  Plans 
ATTN:  DAMO-NC,  Nuc  Chem  Dir 

Dep  Ch  of  Staff  for  Rsch  Dev  &  Acq 
ATTN:  DAMA 


DEPARTMENT  OF  THE  ARMY  (Continued) 

Engineer  Studies  Center 

ATTN:  DAEN-FES,  LTC  Hatch 

Harry  Diamond  Laboratories 
ATTN:  DELHD-NW-P 
ATTN:  DELHD-TA-L 

US  Army  Concepts  Analysis  Agency 
ATTN:  CSSA-ADL 

US  Army  Engineer  Ctr  &  Ft  Belvoir 
ATTN:  ATZA-OTE-ADM 

US  Army  Engineer  School 


ATTN 

ATZA-CDC 

Army  Engr  Waterways  Exper  Station 

ATTN 

R.  Whalin 

ATTN 

WESSE 

ATTN 

WESSD,  J.  Jackson 

ATTN 

J.  Strange 

ATTN 

J.  Zelasko 

ATTN 

F.  Brown 

ATTN 

Library 

ATTN 

WESSA,  W.  Flathau 

ATTN 

WESSS,  J.  Ballard 

US  Army  Foreign  Science  &  Tech  Ctr 
ATTN:  DRXST-SD 

US  Army  Mat  Cmd  Proj  Mngr  for  Nuc  Munitions 
ATTN:  DRCPM-Nuc 

US  Army  Material  &  Mechanics  Rsch  Ctr 
ATTN:  DRXMR,  J.  Mescal  1 
ATTN:  Technical  Library 

US  Army  Materiel  Dev  &  Readiness  Cmd 
ATTN:  DRCDE-D,  L.  Flynn 
ATTN:  DRXAM-TL 

US  Army  Nuclear  &  Chemical  Agency 
ATTN:  Library 

US  Army  War  College 
ATTN:  Library 

USA  Military  Academy 

ATTN:  Document  Library 

USA  Missile  Command 

ATTN:  Documents  Section 
ATTN:  DRSMI-RH 

DEPARTMENT  OF  THE  NAVY 

Naval  Civil  Engineering  Laboratory 
ATTN:  Code  LSI,  J.  Crawford 

Naval  Coastal  Systems  Laboratory 
ATTN;  Code  741 
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^»»CVlous  p 

•s  blank 


David  Taylor  Naval  Ship  R&O  Ctr 
ATTN:  Code  L42-3,  Library 
ATTN:  Code  1700,  W.  Hurray 
ATTN:  Code  1844 
ATTN:  Code  177,  E.  Palmer 
ATTN:  Code  172 
ATTN:  Code  1770.1 
ATTN:  Code  174 
ATTN:  Code  2740 
ATTN:  Code  1740.4 
ATTN:  Code  1740,  R.  Short 
ATTN:  Code  173 
ATTN:  Code  11 
ATTN:  Code  1740.5 
ATTN:  Code  1740.6 
ATTN:  Code  1740.1 

Naval  Electronic  Systems  Command 
AHN:  PME  117-21 

Naval  Explosive  Ord  Disposal  Fac 

ATTN:  Code  504,  J.  Petrousky 

Naval  Facilities  Engineering  Command 
ATTN:  Code  04B 

Naval  Material  Comnand 
ATTN:  MAT  08T-22 


Naval  Surface  Weapons  Center 
ATTN:  W.  Wishard 

ATTN;  Tech  Library  &  Info  Svcs  Br 

Naval  War  Col  lege 

ATTN:  Code  E-11,  Tech  Service 

Naval  Weapons  Center 

ATTN:  Code  343,  FKA6A2,  Tech  Svcs 
ATTN:  Code  266,  C.  Austin 
ATTN:  Code  3263,  J.  Bowen 

Naval  Weapons  Evaluation  Facility 
ATTN:  G.  Binns 
ATTN:  Code  10 
ATTN:  Code  210 
ATTN;  R.  Hughes 

Naval  Weapons  Support  Center 

ATTN;  Code  70553,  0.  Moore 

New  London  Laboratory 

ATTN:  Code  4494,  J.  Patel 
ATTN;  Code  4492,  J.  Kalinowski 

Newport  Laboratory 
ATTN:  Code  EM 

ATTN:  Code  363,  P.  Parantino 


Naval  Ocean  Systems  Center 

ATTN:  Code  013,  E.  Cooper 
ATTN:  Code  4471 

Naval  Postgraduate  School 


ATTN; 

Code  69NE 

ATTN: 

Code  1424,  Library 

ATTN: 

Code  69SG,  Y. 

Shin 

Research  Laboratory 

ATTN; 

Code  8403,  R. 

Bel  sham 

ATTN: 

Code  8440,  G. 

O'Hara 

AHN: 

Code  6380 

ATTN: 

Code  8100 

ATTN; 

Code  8301 

ATTN: 

Code  8406 

ATTN: 

Code  2627 

ATTN; 

Code  8445 

ATTN: 

Code  8404,  H. 

Pusey 

Naval  Sea  Systems  Command 
ATTN:  SEA-033 
ATTN:  SEA- 323 
ATTN:  SEA-06J,  R.  Lane 
ATTN:  SEA-09G53 
ATTN:  SEA- 55X1 
ATTN:  SEA-08 
ATTN:  SEA-0351 
ATTN:  SEA-9931G 

Naval  Surface  Weapons  Center 
ATTN:  Code  F34 
ATTN:  Code  R13 
ATTN;  Code  RIO 

ATTN:  Code  U401,  M.  Kleinerman 
ATTN;  Code  R14 
ATTN;  Code  F31 
ATTN;  Code  R15 


Ofc  of  the  Deputy  Chief  of  Naval  Ops 
ATTN:  OP  987 

ATTN:  NOP  982,  Tac  Air  Srf  &  Ewdev  Div 
ATTN:  NOP  981 

ATTN:  NOP  654,  Strat  Eval  S  Anal  Br 
ATTN :  OP  098T8 
ATTN:  OP  982E,  H.  Lenzini 
ATTN:  OP  957E 

ATTN:  NOP  953,  Tac  Readiness  Div 

ATTN;  OP  37 

ATTN:  OP  225 

ATTN:  OP  03EG 

ATTN:  OP  21 

ATTN:  NOP  951,  ASW  Div 

ATTN:  OP  605D5 

ATTN;  OP  98 INI 

ATTN:  OP  223 

Office  of  Naval  Research 

ATTN:  Code  474,  N.  Perrone 

Strategic  Systems  Project  Office 
ATTN:  NSP-272 
ATTN:  NSP-43 
ATTN:  NSP-273 

DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Institute  of  Technology 
ATTN:  Commander 
ATTN:  Library 

Air  Force  Systems  Command 
ATTN:  DLW 

Assistant  Chief  of  Staff 
Intelligence 
ATTN:  IN 
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DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 

DEPARTMENT 

OF  ENERGY  CONTRACTORS  (Continued) 

Air  Force  Weapons  Laboratory 

Los  Alamos 

National  Laboratory 

ATTN: 

NTES-G,  S.  Melzer 

ATTN 

R.  Whitaker 

ATTN: 

NTE,  M.  Plamondon 

ATTN 

MS530,  G.  Spillman 

AHN: 

NTES-C,  R.  Henny 

ATTN 

Reports  Library 

AHN: 

SUL 

ATTN 

MS634,  T.  Dowler 

AHN: 

NTED 

ATTN 

R.  Sanford 

ATTN 

MS  670,  J.  Hopkins 

Ballistic  Missile  Office 

ATTN: 

DEB 

Oak  Ridge  National  Laboratory 

ATTN 

Civil  Def  Res  Proj 

Deputy  Chief  of  Staff 

Research,  Development,  &  Acq 

ATTN 

Central  Rsch  Library 

ATTN: 

AFRDQI 

Sandia  National  Lab 

ATTN: 

R.  Steere 

ATTN 

Tech  Lib  3141 

ATTN 

L.  Vortman 

Deputy  Chief  of  Staff 

Logistics  &  Engineering 

Sandia  National  Labs,  Livermore 

ATTN: 

LEEE 

ATTN 

Library  &  Security  Classification 

Foreign  Technology  Division 


AHN: 

Nils  Library 

AHN: 

TQTD 

ATTN: 

SOBG 

AUN: 

SD8F,  S.  Spring 

Rome  Air  Development  Center 
ATTN:  RBES,  R.  Mair 
ATTN:  Commander 
ATTN:  TSLD 

Strategic  Air  Command 

AHN:  NRI-STINFO  Library 

DEPARTMENT  OF  ENERGY 

Albuquerque  Operations  Office 
ATTN:  CTID 

Office  of  Military  Application 
ATTN:  OMA/RO&T 

Nevada  Operations  Office 

ATTN:  Doc  Con  for  Technical  Library 

OTHER  GOVERNMENT  AGENCIES 

Central  Intelligence  Agency 
ATTN:  OSWR/NED 
ATTN:  OSR/SE/F 

Department  of  the  Interior 

US  Geological  Survey 
ATTN:  0.  Roddy 

NASA 

ATTN:  F.  Nichols 
ATTN:  R.  Jackson 

US  Nuclear  Regulatory  Commission 

ATTN:  R.  Hhipp  for  Div  Sec  L.  Shao 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

University  of  California 

Lawrence  Livermore  National  Lab 
ATTN:  S.  Erickson 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Applied  Research  Associates.  Inc 
ATTN;  D.  Piepenburg 

Applied  Research  Associates,  Inc 
ATTN:  B.  Frank 

BDM  Corp 

ATTN:  T.  Neighbors 
ATTN:  A.  Lavagnino 
ATTN:  Corporate  Library 

California  Institute  of  Technology 
ATTN:  T.  Ahrens 

California  Research  &  Technology,  Inc 
ATTN:  M.  Rosenblatt 
ATTN:  S.  Schuster 
ATTN:  Library 
ATTN:  K.  Kreyenhagen 

Columbia  University 

ATTN:  H.  Bleich 
ATTN:  F.  Dimaggio 

University  of  Denver 

ATTN:  Sec  Officer  for  J.  Wisotski 

Electric  Power  Research  Institute 
ATTN:  G.  SI  iter 

Electro-Mech  Systems,  Inc 
ATTN:  R.  Shunk 

General  Dynamics  Corp 
ATTN:  J.  Miller 
ATTN:  J.  Mador 
ATTN:  H.  Pakstys 

Kaman  AviDyne 

ATTN:  R.  Ruetenik 
ATTN:  G.  Zartarian 
ATTN:  Library 
ATTN:  N.  Hobbs 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued 


Kaman  Sciences  Corp 
ATTN:  Library 
ATTN:  F.  Shelton 

Kaman  Sciences  Corp 
ATTN:  0.  Sachs 

Kaman  Tempo 

AHN:  DASIAC 

Karagozian  and  Case 

ATTN:  J.  Karagozian 

Lockheed  Missiles  &  Space  Co,  Inc 

ATTN:  Technical  Information  Center 
ATTN:  T.  Geers 
ATTN:  B.  Almroth 

Lockheed  Missiles  1  Space  Co,  Inc 
ATTN:  TIC-Library 

M  &  T  Company 

ATTH:  0.  McNaight 

Management  Science  Associates 
ATTN:  K.  Kaplan 

McDonnell  Douglas  Corp 
ATTN:  R.  Halprin 

NKF  Engineering  Associates,  Inc 
ATTN:  R.  Belsheim 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 

Pacifica  Technology 
ATTN:  R.  Bjork 
ATTN:  G.  Kent 
ATTN:  A.  Kushner 

Physics  Applications,  Inc 
ATTN:  C.  Vincent 

Science  Applications,  Inc 

ATTN:  Technical  Library 

Southwest  Research  Institute 
ATTN:  A.  Wenzel 
ATTN:  W.  Baker 


Physics  International  Co 
ATTN:  L.  Behrmann 
ATTN:  F.  Sauer 
ATTN:  J.  Thomsen 
ATTN:  E.  Moore 
ATTN;  Technical  Library 

S-CUBED 

ATTN:  T.  Cherry 
ATTN:  R.  Sedgewick 
ATTN:  D.  Grine 
ATTN:  T.  Riney 
ATTN;  Library 
ATTN:  K.  Pyatt 
ATTN:  T.  McKinley 

SRI  International 

ATTN:  G.  Abrahamson 
ATTN:  W.  Wilkinson 
ATTN:  A.  Florence 

Teledyne  Brown  Engineering 
ATTN:  J.  Ravenscraft 

Tetra  Tech,  Inc 

ATTN:  L.  Hwang 

TRW  Electronics  &  Defense  Sector 
ATTN;  P.  Bhuta 
ATTN:  A.  Feldman 
ATTN:  N.  Lipner 

ATTN;  Technical  Information  Center 
ATTN:  D.  Jortner 
ATTN:  B.  Sussholtz 

TRW  Electronics  &  Defense  Sector 


ATTN: 

P. 

Dai 

ATTN; 

F. 

P leper 

ATTN: 

E. 

Wong 

ATTN: 

G. 

Hulcher 

Weidlinger  Assoc,  Consulting  Engrg 
ATTN:  J.  McCormick 
ATTN:  H.  Baron 
A  cy  ATTN:  R.  Vasudevan 
4  cy  ATTN:  D.  Ranlet 

Weidlinger  Associates 

ATTN:  J.  Isenberg 
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